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the symbol Tj. 

Due to the structure of the data achieved on the basis of the results of investigations of a real 
system, being indispensable to determine the semi-Markov process {x „LET is it is convenient to 
define such a process by determining the following three (p, P, F'(t)) elements: 

- p= (p, ieS ] - stochastic vector of the initial distribution process {x „LET 2 

- P=[p,:i,7¢S] - matrix of transition probabilities of the Markov chain embedded in the 
process Xt € Tt. 

- F (t) =[F; (t) ii,jeS] - matrix of distribution functions of the conditional distributions of the 
random variables Tij, i,j e S of the times remaining in the states. 

To simplify further investigations it has been assumed that F;(t) = G;(t), i, j e S. The function 
G;(t) is a distribution function of distribution of the duration time of the state ieS. This 
assumption means that the duration of the current state does not depend on the next process state. 

The state space S consists of five states S = (51, S2, S3, S4, S5} in the considered example. The 
possible transitions between the selected bus maintenance states, presented by the relation (15), 


were determined on the basis of identification of a real investigation object (a maintenance system 
of an urban bus transport system) and the maintenance processes performed therein. 
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P=| P3 Pa 0 Py 0 (15) 


LPs; 0 0 py 0 
The transition matrix of the embedded Markov chain estimated on the basis of the results from 

the performed preliminary maintenance investigations is formulated as follows: 

0 0,9998 0 0,0001 0,0001] 


0,986 0 0,013 0,001 0 


P= 011 0,84 0 0,05 0 (16) 
1 0 0 0 0 
| 0,9999 0 0 0,0001 0 


The values of the estimators of the maximum credibility, formulated as stated below, have been 
adopted as statistical estimations of the transition probabilities of the embedded Markov chain: 


Py (17) 


where: 
nij- the number of the direct chain state changes, from the state i into the state j in its finite 
process realization. 
The determined limiting distribution of the Markov chain embedded in the considered semi- 
Markov process is presented in the Table 1. 


Table 1. Limiting distribution of the embedded Markov chain 


la | e | s | m | s | 
0,00648662 


The expected values E(T;), i e S of the distributions determined by the empiric distribution 
function are the estimations II; of the expected value of the time T; of duration of the state i € S: 
1 i m 
IN, =) 6”) (18) 
n 


i m=l 


(1) (n;) 
ths 


i 


independent realizations of the random variable T; which stands for the 
duration of the state 1. 

The expected times of duration II; of the considered maintenance states (expressed in hours) 
estimated on the basis of the preliminary data are presented in the Table 2. 


Table 2. Expected values of times of duration of the analysed maintenance states 


[6 | 15.7608 | 0.0456 | 216 | 0,0336 | 
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The determined limiting distribution of the semi-Markov process constituting the model of the 
maintenance process of a single bus is presented in the Table 3. 


Table 3. Limiting distribution of the analysed process 


Boe | ee e | ee 


4. Model of maintenance process of a set of vehicles 


Assuming that n of homogenous vehicles with the codes k e N, forming the set Z of vehicles, 
are used in the maintenance system, a simplified model of the state change process of the elements 
in this set has been developed. 

It has been assumed that the independent semi-Markov processes {X,(t):t 20}, 
ke N= fl... n} determined by identical kernels Q(t) described with the relation (4) are the model 
of the bus maintenance process. The assumption that the processes are independent is just an 
approximation of the real maintenance processes. In fact, some bus states may have influence over 
the maintenance states of other buses. 

The random vector: 

X(t) = [Xi(t), X2), ..., Xa(t)] (19) 
describes the process of changes of maintenance states of the set of buses of the considered vehicle 
maintenance system. With such an assumption it is possible to determine the total distribution of 
this random vector, that is: 

PiXr(0 = i, X2(t) = i2, ..., Aa) = in} = Pi (t) P ... Pi (t), (20) 
where: i1, i2 ..., in € S. 

The random variable: 

Ni(t)=#{k e N: X(t) =i}, t20, ie S (21) 
stands for the number of buses in the state i at the moment t. 
The random vector {(Ni(t), Na(t), ... , Ng(t)):t20} has multinomial distribution, which 
means: 
| : : 
PIN 6) =i, NO =i,,..., N0) =i} = al OFP OP,- [POT (22) 
alili! 
where: i) +i2 +... +i5 =n. 
The equation [7] results from the limit theorem for the semi-Markov processes: 


lim PiN, (1) =i, Ni) =i, N;(ł)=i;j = 


n! . i i i n! ipi is 
=F lm [FOP TP, OF". LPT? cra ee Borak: = 


ii,!..i;! WRZE 
hee l B b is is 
_ nlm m"a, mą” ...5 m; 


ititi Tem | 


ieS 


(23) 


The expected value of the considered random vector is formulated as: 
E[(Ni(t), N2(t),..., Ns(t))] = (nP;(t), nP>(t),..., nPs(t)). (24) 
The coordinates of this vector stand for the expected number of buses in the analysed 
maintenance states i, i=1,2,...,5, at the moment t. Therefore, the value: E[(N;(t)] = nP;(t), 
i € S, t 2 O stands for the expected number of buses in the state i at the moment t. 
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The variance of the random variables N; takes the form: 


VIN] =nPi(HU-Pi®), ieS,t20, (25) 
while the covariance is expressed with the relation: 
Cov[(Ni(t), Ni] =- nP- Pj), iyEeS,i4j, t20. (26) 


The expected value of the random variable N; which stands for the number of buses in the state 
i in a system working for a long time may be determined using the relation: 
ELN,]=limE[(N,())]=nP, ieS. 27 
The expected value of the random vector [N;, Na, ..., Ns] in the considered model is: 
E[(Ni, N>...., Ns)] = (nP, nP>,..., nPs). (28) 
The coordinates of this vector stand for the expected number of buses in the determined states 
in a system working for a long time. The statistical estimation of the elements of this vector is: 


E(Ni) im n Pi, ieS. (29) 
The limiting variance for the state i € S is: 
V[N;]=nP;(1-P;), ieS, (30) 
while the limiting covariance of the state pairs is expressed with the relation: 
Cov[(Ni(t), (Nj(t)] = - nP;(1-P;),ijeS,iźj. (31) 


The determined expected numbers of buses in the specified states in a system working for a 
long time for the investigated bus maintenance system and the adopted data and n=203 are 
presented in the Table 4. 


Table 4. Statistical estimation of the expected number of buses in the specified states in a system working for a long 
time (n = 203) 


[35.53 | mas | oor | 004 | 000 


It is possible to simulate such an event that the scope of the transport service realization is 
reduced. Such a situation may take place when it is not profitable to realize transportation by 
means of specific bus lines or in case of a failure to win another tender for the transportation tasks 
performed before. If it is decided not to handle one of the bus lines, the number of the necessary 
vehicles would be reduced by 17 vehicles. For such a variant, the determined expected values of 
the number of buses in the specified states in a system working for a long time are presented in the 
Table 5. The calculation results may serve to preliminarily forecast the demand for the corrective 
service realization both in the system (service station) and in its environment (technical emergency 
service units), as well as the system operation costs and the number of vehicles necessary to 
provide continuous transportation. 


Tab. 5. Statistical estimation of the number of buses in the specified states in a system working for a long time 


(n = 186) 


| 50,88 | 135,09 | 001 | 003 | 0,00 | 


5. Summary 
The purpose of the investigations presented herein was, among others, to show a possibility to 


use selected stochastic processes for mathematical modelling the system and the process of vehicle 
maintenance. 
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The considered example of the model of maintenance process of urban transport buses is 
characterised by a considerable simplification (due to the nature of the work). However, the 
presented way to create models of that type and analyse them indicates that they may be used to 
preliminarily forecast the system state, and to simulate changes of the internal and external 
interactions affecting the system and to support the decision makers in the decision making 
process concerning control of the maintenance process and the system in which it is performed. 

The analysis of results of the investigations performed shows that the developed model is 
susceptible to a change of the input parameter values. 

The mathematical models of the maintenance processes, performed in complex systems, are by 
their very nature a considerable simplification of the real processes. The consequence of this is a 
necessity to carefully formulate practical conclusions resulting from investigating those models. 
However, in the opinion of the author, it does not limit the purposefulness of developing and 
analysing models of that type. The analysis of the results of investigating those models, for the 
model parameter values, determined on the basis of the results of maintenance tests performed in a 
real transport system, makes it possible to formulate conclusions and opinions both of a qualitative 
and (to the limited extent) quantitative nature. 

Appropriate formulation of the investigation goals, expected results and variants of simulation 
of behaviour of the analysed system makes it possible, due to analysing the models of that type, to 
perform, among other things, preliminary evaluation of influence of the considered decision 
variants (concerning such things as the number of vehicles required to assure effective realization 
of the transport tasks, frequency and scope of the technical services performed, etc.) on the course 
and effectiveness of the investigated maintenance process. 

The presented method to build models of changes of vehicle maintenance states and the semi- 
Markov process used to model the maintenance process of urban transport buses may be, in a 
relatively simple way, applied for other maintenance systems than the considered one. 


References 


[1] Bobrowski D.: Mathematical models and methods of reliability theory in examples and 
exercises. WNT, Warsaw 1985. 

[2] Bobrowski, D., Probabilistics in technical applications. WNT, Warsaw 1986. 

[3] Buslenko, N., Kałasznikow, W., Kowalenko, I., Theory of complex systems. PWN, Warsaw 
1979. 

[4] Koroluk, W.S., Turbin, A.F., Połmarkowskije processy i ich priłożenija. Naukowa Dumka, 
Kijów 1976. 

[5] Landowski, B., Maintenance model of a certain class of technical objects. Scientific 
Journals No. 229. Mechanics 48. University Publishing House ATR, Bydgoszcz 2000. 

[6] Mine, H., Osaki, S., Markowskije processy priniatija reszenij. Science, Moscow 1977. 

[7] Woropay, M., Grabski, F., Landowski, B., Semi-Markov model of the vehicle operation 
and maintenance processes within an urban transportation system. Scientific Publishers 
PTNM, the Archives of Automotive Engineering Vol. 7, No. 3, 2004. 

[8] Woropay, M., Knopik, L., Landowski B., Modelling the maintenance and operation 
processes in a transport system. Library of Maintenance Problems. ITE, Bydgoszcz - Radom 
2001. 

[9] Woropay, M., Landowski, B., Neubauer, A., Controlling reliability in transport systems. 
Library of Maintenance Problems. ITE, Radom 2004. 

[10]  Woropay, M., Landowski, B., Simulation analysis of the maintenance and operation 
processes in an urban transport system. Scientific Journals No. 212. Mechanics 42. University 
Publishing House ATR, Bydgoszcz 1998. 


168 


Journal of POLISH CIMAC 


Faculty of Ocean Engineering & Ship Technology 
GDANSK UNIVERSITY OF TECHNOLOGY 


EFFECT OF SUPPLY VOLTAGE ON THE DOSAGE OF FUEL 
IN INJECTION SYSTEM THE COMMON RAIL TYPE 


Kazimierz Lejda, Adam Ustrzycki 


Rzeszow University of Technology 
Al. Powstańców Warszawy 8, 35-959 Rzeszów, Poland 
tel.: +48 17 8651597, +48 17 8651531 
fax: +48 17 8543112 
e-mail: klejda(QQprz.edu.pl, austrzyc@prz.edu.pl 


Summary 


Common Rail systems are most frequently used solution to the fuel supply system of Diesel engines. This is due to 
due to their high potential to shape the characteristics of fuel injection. One of the important requirements for each 
injection system is precision dose control and its onset, that is, as a consequence of the start of the combustion 
process. The most important parameters affecting the accuracy of dosing injection systems are, of course injection 
time, pressure of fuel in the rail and the fuel temperature. Influence on the size of the dose but also other factors, 
which include the supply voltage, which may be subject to significant change during the start-up of engine. The issue 
of start-up is becoming of interest to EURO standards, as the processes occurring in the course have 
a significant impact on the level of pollutants in the exhaust of internal combustion engine. Increased in this period, 
emissions of toxic compounds results from the unstable processes in the engine cylinder, as well as dynamic changes 
in the parameters of injection and in the start-up system, in particular, changes in voltage resulting from a large 
collection of current from the battery. Thus, the article attempted to determine the impact of voltage supply on the 
dosage of fuel in different operation conditions of the injection system. In particular attention was paid to the size of 
the dose and the actual changes at the beginning of fuel injection voltage changes. The study was conducted on a test 
stand, which were the main components of the test bench Bosch EPS-815 with electronic measurement of fuel dose 
and the visualization system AVL Visioscope. Control of the injector was performed using the self-developed 
controller. 


Key words: Common rail, fuel injection, control injection, precision dosing of fuel 
1. Introduction 


Common Rail fuel supply systems have high potential in shaping the fuel injection 
characteristics. For that reason, they are now widely used systems of fuel injection for diesel 
engines [3, 7, 8]. Diesel injection system is expected to provide a high precision metering of fuel 
amount injected into the cylinder, accuracy in injection start timing control and repeatability of 
that. Many factors contribute to above capacities including, inter alia, accurate control of injection 
duration and timing, which depends on the electronic control unit design and the kind of input 
signals that it uses, the injection pressure, the fuel parameters [4, 5, 10, 11], and the power supply 
voltage. On the accuracy of injected fuel amount, the factors related to the injection process 
progressing have also a significant impact, especially when multi-split injection is implemented. 
Pressure pulsations triggered by pilot injection, the amplitude of which may reach even 60 MPa, 
change the fuel quantity of main injection and its start angle position as well [9]. The changes in 
voltage, especially during the engine start-up, can be so considerable, that in effect they may cause 
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the changes of the injector coil current value and emerged magnetic stream. Fluctuations of this 
stream will affect the magnetic force lifting the injector control valve, so the entire range of 
injector operation. 

In this article the impact of supply voltage on the injected fuel quantity and the actual start of 
fuel injection are studied. 


2. Test stand and measurement methodology 


To determine the effect of voltage variations on the dosage of fuel, the test stand was built, 
which consist of a few key devices like: the test bench Bosch EPS-815 joined to the fuel delivery 
electronic measurement device Bosch KMA-822 with heat exchanger, and the module for testing 
the Common Rail system CRS-845. The fuel was injected into the special transparent chamber 
filled with the fuel, that enables the observation of injected fuel sprays. Fuel propagation inside the 
measuring chamber was recorded with use of visualization system AVL Visioscope. It was 
captured at every 0.1° of pump shaft angle. The layout of the test stand is shown on Fig. 1. 
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Fig. 1. The test stand layout: 1-test bench Bosch EPS 815, 2-high pressure rail, 3-fuel delivery measurement device 
Bosch KMA 822, 4-control module of high-pressure pump, 5,15- piezoresistive pressure sensors, Kistler 4067A2000, 
6- rail pressure sensor, 7- rail pressure regulators, 8-fuel temperature sensor, 9-high pressure pump 10-pump shaft 
position and rotational speed optical sensor of AVL 365C, 11-light pipe, 12-digital camera PixelF ly VGA, 
13-endoscope, 14-visualization chamber, 16-tested injector, 17-measurement connector of injector control voltage, 
18-injector current sensor PA-55, 19- microprocessor controller of injector, 20, 21-charge amplifiers Kistler 4618A0, 
22-light unit, 23-computer with control software of injector,24-visualisation system AVL Visioscope, 25-computer with 
software for test bench controlling and for measuring of fuel delivery, 26-computer with data acquisition devices for 
high-speed parameters recording, 27,28-temperature sensors, 29,30-termometers EMT 101 
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Research in each measurement series were carried out at a fixed temperature of the injected 
fuel, which was measured in outlet pipe of the measuring chamber. The temperature of the fuel in 
the tank was stabilized at around 40 + 1°C. As the testing device the Bosch injector from the 
engine Fiat Multijet 1.3 (marked 445 010 083) was used. For injector operation control the 
electronic unit developed in Department of Automotive Vehicles and Internal Combustion Engines 
of Rzeszów University of Technology was used. Detailed description of the controller is presented 
in the work [1]. 

The primary purpose of the study was to determine the injected fuel amount and actual start of 
fuel injection depending on supplying voltage. The study was conducted for single injection 
strategy in order to avoid the changes in injection delivery that may arise from variations of 
pressure due to implementation of multi-split injection [2,6]. In order to determine precisely the 
actual start of injection, the progress of fuel spray was recorded by visualization system AVL 
Visioscope with an angular resolution of 0.1° of pump shaft revolution. That accuracy level was 
possible to obtain owing to use of pulse angle encoder AVL 365C mounted on the shaft of test 
bench. Pulses generated by the transmitter were also used as the input signal for the electronic 
injector controller. In this case, the angular resolution of 720 pulses per revolution was available, 
what gives a very high precision of controlling the start of fuel injection. 

For recording the progress of injection process a digital VGA camera and the endoscope of 
AVL Visioscope system were used. Figure 2 presents the photographs of fuel sprays in the 
measuring chamber that were captured in the same angle position of pump shaft but for different 
values of supply voltage. 


-51.0 deg CA Supply voltage U=9 V 51.0 deg CA Supply voltage U=11 V 


Fig. 2. Progressing frames of fuel spray flowing out the 
tip of Fiat Multijet injector recorded at the same angle 
position of pump shaft for different values of supply 
voltage (pressure in the rail p, = 125 MPa, 
injection time t,, = 1 ms, n = 1000 rpm) 


610 deg CA Supply voltage U=13 V 
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As can be noticed, at the supply voltage of 9 V the fuel sprays there are at the initial phase of 
their progress, while at the supply voltage of 13 V they are already in mid-expanded stage. 

The injector control system was powered from a battery, which was also connected to the 
charger with a free and precise adjustment of supplying voltage and current feature. It allows for 
the voltage power supply to stay highly stable and reliable. The research were conducted for the 
voltage varied from 9 to 13 V, beginning from lower value. In order to perform it, the battery was 
partially discharged and loaded by an additional resistance increasing current consumption. Since 
the study was conducted at different durations of the injection, which causes the changes in the 
injector input current, each time after changing the injection time the supply voltage was re- 
adjusted. The tests were made for injection durations changed in the range of 1 to 3 ms, for three 
values of fuel pressure in the whole system equals 75, 100, and 125 MPa and at rotational speed of 
1,000 rpm. In addition, for specific values of injection time and pressure some studies on the 
impact of supply voltage on the fuel injected quantity at different pump shaft speeds changed from 
600 to 2000 rpm have been carried out. 


3. Test results 


Fig. 3 presents the volume of fuel injected versus the time of injection for different values of 
voltage U and fixed rotational speed n = 1000 rpm and pressure in the rail 75 MPa. For the voltage 
below 11 V this relationship ceases to be linear, and further, for the voltage below 9.5 V, the 
amount of fuel injected becomes a random, and increasing the time of injection does not make it 
higher. Thus, these results are not included in the figure. Reducing supply voltage results in a 
delay in opening the electromagnetic valve and thereby the start of injection delays also (Fig. 4). 
For the change in supply voltage from 13 V to 9.5 V, the delay at this rotational speed is exceeds 
1°of pump shaft angle. 
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Fig. 3. The fuel dose quantity Q versus the injection time ti; at different values of voltage U 
(the pressure in the rail Pai = 75 MPa, n = 1000 rpm) 


Fig. 5 presents the changes of the value of fuel dose quantity Q depending on the supply 
voltage for different values of fuel pressure in the rail and for fixed injection time of 1 ms. The 
temperature of the fuel flowing out the chamber was 90°C for the injection pressure of 100 and 
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125 MPa, but 75°C for the injection pressure of 75 MPa due to the fact, that attainment 90°C of 
fuel temperature at this value of injection pressure was not possible. 

As can be seen, the variation range of injected fuel quantity is very wide; it changes about 
5 mm '/stroke per each 1 V of drop, where at the voltage less than 10 V this value still increases. 
For longer injection time of 3 ms, these changes fit in the exponential curve; in this case a 
significant reduction in the injected fuel amount is already noticed at a voltage of 11 V (Fig. 6). 
Minor changes of fuel amount at shorter injection time can be explained by the influence of 
capacitor, whose stabilizing action becomes little at increasing injection time. 
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Fig. 4. Changes in the actual start of injection a versus the injection time tin for different values 
of voltage U (the pressure in the rail P,a = 75 MPa, n = 1000 rpm) 
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Fig. 5. Effect of supply voltage U on the fuel injected quantityQ for different values of pressure in the rail Pyqij 
(injection time ti; = 1 ms, n = 1000 rpm) 
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In Fig. 7 the relationship between injected fuel amount Q and pump shaft speed n for three 
different values of supply voltage of 9, 11, and 13 V, injection time of 3 ms and the fuel pressure 
of 75 MPa is presented. At 9 V, the injected fuel amount corresponds to only about 15% of the 
amount at higher voltage values. In principle, the injected fuel quantity does not dependent on the 
speed, although for larger values of supply voltage a slight increase of injected fuel volume at 
higher rotational speeds can be observed. 
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Fig. 6. Effect of supply voltage U on the injected fuel amount Q for different values of pressure in the rail Dai 
(injection time tin; = 3 ms, n = 1000 rpm) 
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Fig. 7. The injected fuel amount Q versus the pump shaft speed n for different values of voltage U 
(injection time tin; = 3.0 ms, the pressure in the rail P,ait = 75 MPa) 


Figure 8 presents the changes in the start of injection angle position depending on the pump 
shaft speed. As can be seen, the time delay of needle lift associated with a reduction of voltage 
value is constant, what results in increasing the value of changes in injection beginning angle 
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position at increasing pump shaft speed. However, throughout the whole analysis range the change 
of voltage from 11 V to 9 V causes greater differences in the injection beginning angle position 
than as for the change from 13 to 11 V. In the same conditions of system operation, similar 
changes occur for a longer injection time. 
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Fig. 8. Changes of the injection start angle position a depending on the pump shaft speed n for different values of the 
voltage U (injection time tiny = 1 ms, n = 1000 rpm) 


4. Conclusions 


Analyzing the results obtained, it should be noted that the change in supply voltage 
significantly affects the injected fuel quantity. The lower injection pressure, the more significant 
the relative changes of that are. It may further affect the injected fuel amount, especially at low 
engine speeds. Changes in the angle position of the injection start for low speeds are negligible, 
and taking into account that at high speeds and alternator order working, the voltage fluctuations 
are smooth and should not be large, so their impact on this parameter can be neglected.. 

However, the injector may not open at all, if the voltage would be to low, especially at low fuel 
temperature. The injector operation principle at low values of voltage depends to a large extent on 
the kind of capacitor used in the controller. The boundary values of voltage, at which the injector 
will still work, also depend on the parameters of the injector coil, and so, also on the coil 
temperature, that may change due to the temperature of the fuel, as well as working time of the 
injector. 

In conclusion, the voltage fluctuations can significantly affect the value of fuel quantity 
injected by the Common Rail system into the engine cylinders, which may adversely affect the 
emissions of pollutants in the engine exhaust. The elimination of these changes can be obtained by 
applying an additional electricity storage device e.g. by placing the suitably selected capacitor in 
the controller or by the use of a separate power source for the injectors, so that the change of the 
battery voltage during engine start-up and immediately after does not affect the accuracy and 
quality of injection control in the meaning of fuel dosing in these phases of Diesel engine work. 
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Abstract 


Reliability needs of the military focus primarily on operational readiness, longevity, supportability, 
robustness. That means product performance on demand, long useful life, repair ability and satisfactory performance 
over environmental extremes. Availability requirements of modern systems and equipment demand that designing for 
maintainability must be as important as designing for other performance characteristics. 


Keywords: helicopter, exploatation, gas turbine engine 


Customer needs are always very high. The challenge to the manufacturer or service provider is 
how to assess and define true customer expectations and then how to design, manufacture and sell 
the product to best meet those expectations. Military customer is specific. The air force is very 
specific but navy aviation is extremely specific. There are numerous differences between the needs 
of a typical customer and a navy aviation customer. Reliability needs of the military focus 
primarily on operational readiness, longevity, supportability, robustness. That means product 
performance on demand, long useful life, repair ability and satisfactory performance over 
environmental extremes [1]. A thorough understanding of the environment in which the product is 
intended to operate is the major consideration in determining the ability of a product to meet the 
operational reliability needs of the military. The environment is a critical factor in determining the 
suitability of the equipment for military use, especially in the navy aviation. Each product must 
provide operational capabilities which allow military forces to maintain technological and logistics 
superiority over a potential adversary even if it is just natural or induced environment which can 
be expected over the operational and maintenance portions of life cycle. Degradation or loss of 
these operational capabilities typically can result in immediate loss of human life and equipment 
(Fig. 1, 2) [2]. 
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Fig. 1 


Availability requirements of modern systems and equipment demand that designing for 
maintainability must be as important as designing for other performance characteristics. According 
to military definition, maintainability is a measure of the ability of an item to be retained in or 
restored to a specified condition when maintenance is performed by personnel having specified 
skill levels, using prescribed procedures and resources, at each prescribed level of maintenance or 
repair. You can find some other definitions connecting effectiveness and economy and efficiency. 


Fig.2 


LAY) 


The T700-GE-401 engine is a front drive, turbo shaft engine featuring a single spool gas 
generator section consisting of a five-stage axial, single stage centrifugal flow compressor, a 
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through—flow annular combustion chamber, a two-stage axial flow gas generator turbine, and a 
free or independent two-stage axial flow power turbine. The power turbine shaft, which has a 
rated speed of 20,900 rmp, is co-axial and extends to the front of the engine. The engine 
incorporates an integral inlet particle separator, a top-mounted accessory package, an engine — 
driven fuel boost pump, a self-contained lubrication system, condition monitoring, diagnostic 
provisions, a hydro mechanical control system, and an electrical control system or digital 
electronic control system (401C version) providing power turbine speed control, dual engine load 
sharing, temperature limiting, and redundant power turbine over speed protection. The T700-GE- 
401 engine also has a contingency power rating for use under one engine inoperative conditions. 

Engine control system has all the components necessary for the proper and complete control of 
the engine. This system provides for fuel handling, fuel computation, compressor bleed and 
variable geometry control, PT over speed control protection, over temperature protection and 
torque matching. This system is self-contained and does not require external electrical power for 
any controlling functions. It also allows replacement of any weapon replaceable assembly (WRA) 
without calibration. The following units make up the control system: 

- Hydro mechanical control unit, The HMU schedules fuel for combustion. It contains a 
high-pressure pump. The HMU has an actuator arm that positions the inlet guide vanes, the 
variable compressor vanes, and the anti-icing bleed and start valve. The HMU responds to 
engine inlet air temperature (T2), to compressor discharge air pressure (P3), and to a trim 
signal from the ECU/DEC to set fuel flow and variable vane positioning. It is designed to 
be adjusted at depot only. Adjustments to the HMU are safety wired to prevent adjustment 
at intermediate maintenance level, 

- Electrical control unit (ECU) /Digital electronic control (DEC), Used for diagnostic testing 
of the engine electrical control system, interface with the test stand system and some 
sensors, 

- Fuel boost pump. The fuel boost pump is an engine-mounted suction-type pump. It is not 
self priming. 

- Fuel filter. The main parts of the fuel filter are: a disposable filter element, an impending 
bypass indicator button, a bypass relief valve, and an actual filter bypass sensor. 

- Overspeed and drain valve which has three functions. First, it sends fuel through the main 
fuel manifold to the fuel injectors for starting, accelerating, and engine operation. Second, 
it purges the fuel injectors of fuel when the engine is shut down. It does this by allowing 
compressor discharge air (P3) to pass through the injectors. Third, it controls Np overspeed 
by cutting off engine fuel flow. 

- Anti-icing bleed and start valve 

- Ignition exciter. The ignition system is a noncontiguous, ac-powered, capacitor discharge 
type. It includes two igniter plugs, two electrical ignition leads, and an ignition exciter 
assembly. Power is supplied to the ignition exciter assembly by the engine's alternator 
during starting only. 

- Np sensor: The power turbine rotor speed (Np) governing system compares the signal sent 
from the Np sensor with the Np selected by the operator. It varies fuel flow by actuating 
the torque motor in the Hydro mechanical control unit (HMU). Power turbine rotor speed 
is governed with +1% of required rotor speed. 


There are some engine control components, test cell systems and diagnostic equipment. 
Turbine gas temperature (TGT) limiting system. The TGT limiting system overrides the power 
turbine rotor speed governing system and the load-sharing system when TGT reaches 1542°F. It 
limits fuel flow to hold a maximum 1542°F TGT by actuating the torque motor in the HMU. The 
TGT limiting system is accurate to within +9°F. Np over speed protection system. The Np over 
speed protection system (fig. 12) receives a power turbine speed signal from the torque and over 
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speed sensor. When Np exceeds 25,000 + 250 rpm, output from the protection system activates a 
solenoid in the ODV. This shuts off fuel flow causing the engine to shut down. Load-sharing 
system. In twin-engine installations, the digital electronic controls (DEC) compare torque signals, 
for automatic load sharing. Hot Start Prevention (HSP). The HSP system prevents over 
temperature during engine start. The HSP system receives power turbine speed signal, gas 
generator speed signal, and turbine gas temperature (TGT). When Np and Ng are below their 
respective hot start reference and turbine gas temperature (TGT) exceeds 1652°F and output from 
the HSP system activates a solenoid in the ODV. This shuts off fuel flow and causes the engine to 
shut down. Fault Indication. The DEC contains signal validation for selected input signals within 
the electrical control system. Signals are continuously validated when the engine is operating at 
flight idle and above. If a failure has occurred, the failed component or related circuit will be 
identified by a pre-selected fault code. It is possible to have more than one fault detected. Each 
code should be treated as an individual fault. It should be noted that the signal validation does not 
recognize aircraft instrument related failure. Fault codes will be displayed on the engine torque 
indicator. 

Next the TW3-117M engine for the Mi-14 helicopters (Fig. 4) designed by S. Izotow team at 
Klimow Company. The turbine converts the potential energy of gas to mechanical work of the 
rotor that can be used by the consuming device. New materials have enforced some improvements 
to increase the gas temperature before the turbine without changing the geometry of the flow 
section, and the takeoff power. Further improvements are made to develop an advanced cooling 
system, use new materials including cermets, and reduce the weight of the turbine. 


Each type of engine has its own manual containing introductory information necessary to 
assure proper maintenance of the system. It includes all items having an approved mandatory 
component replacement interval and those items requiring Conditional inspections that shall be 
accomplished at the occurrence of an over limit situation. Conditional Inspection requirements 
include a brief description of what is to be performed and a reference to the manual or directive 
containing detailed requirements. Replacement items are indicated in unit hours, calendar time, 
cycles or event and are arranged by engine systems. All T700 engine hardware is tracked with 
Aeronautical Equipment Service Record. There are some cases like: salt water landing, hard 
landing, over—water hoisting operation, over—torque, overspeed, sudden stoppage, over 
temperature, oil pressure loss, engine operation limits exceeded, required conditional inspection in 
compliance with technical manual cards. The Chief Engineer established maintenance program for 
Naval Air Force which outlines the system of three levels of maintenance: organizational level (O- 
level), Intermediate (I-level) and depot-level (D-level) maintenance. The most interesting is the I- 
level maintenance. The objective of to enhance and to sustain the combat readiness and mission 
capability of support at the nearest location with the lowest practical resource expenditure. It 
consists of: maintenance on aeronautical components, calibration of designated equipment, 
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processing of components from non—mission capable planes, incorporation of technical directives, 
age exploration of a plane. 

Correct and very precise diagnose and troubleshooting need technicians with years of 
experience on aircraft [2, 3]. Training is given with the emphasis on real-life incidents as a result 
of a human error and structure fatigue. The key is our ability to manage the volume and diversity 
of information required to keep engines running, and easy and fast access to that for those making 
repairs. 
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Abstract 


Safe operational use of an technical system, depends on properties and characteristics thereof, external 
excitations, and the correctly performed process of operation/maintenance of this system. Unserviceability and 
failures/damages to the system are the most essential problems, which — if not identified on time — can result in an 
accident (a crash or failure). The issues of safety, unserviceability and failures/damages have been exemplified with 
those occurring to gas turbine engines. 


Keywords: exploatation, failures, gas turbine engine 
1.Introduction 


Causes of damages are pretty well recognized owing to rich maintenance experience and 
examination of engineering objects and units thereof. They are usually deeply rooted in physical 
and chemical phenomena (effect of the process of the IV type, according to description in [6.13]), 
but can also be effected by maintenance activities that may provoke the boundary limits, e.g. of 
strength/resistance, wear-and-tear limits, etc. of particular components are exceeded. Causes of 
damages/failures are grouped into four classes: direct, indirect, primary, and secondary ones. 

Effects of damages have their impact on the performance of operational tasks, reliability, 
safety, and economic effectiveness of operation/maintenance. The more important component of 
an engineering object (from the point of view of its functionality), the more serious effects brought 
about by a damage/failure to this component. Damages/failures may occur in the course of 
operational use of the object as well as the maintenance thereof. It is obvious that damages to an 
engine (its structural component or unit) while it is operated in the air are essentially different from 
those occurring on the ground. 


2. Damages/failures to engines 
The term ‘engine damage’ comprises damages to the engine structure, units as well as 
particular components (structural members and parts) that the units are composed of. Damages to 


parts and structural members may bring about damages to the whole unit/assembly in the 
composition of which they enter; also, to other units/assemblies of a given object. 
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Any damage is defined as the loss of physical properties and/or peculiar features by the 
structural component, i.e. an event that results in the transition of the engineering structure from 
the fit-for-service condition into the unfit-for-service condition. 

The most essential criterion of rating an engineering event among damages/failures is a change 
in a feature/a set of features or its/their deviation from those, on the grounds of which the 
structural component/unit is claimed to be either ‘fit-for-use’ or ’unfit-for-use’. 

From the standpoint of the engineering-object operation, damages/failures may be either 
typical, i.e. well recognized and identified, described in operational/maintenance documentation, 
or new ones, or those taking new forms. The latter compose a new set and occur in the course of 
current operation/routine maintenance. These new forms of damages/failures or of the object’s 
remaining unserviceable may prove single events or system-attributable ones. In operational and 
maintenance practice they are responded to with suitable preventive actions. 

Damages/failures occur due to the willfully or unwittingly exceeded permissible loads (e.g. by 
the aircrew, operator, servicing staff), externally originated effects and disturbances, or natural 
physical and chemical phenomena typical of the performed tasks (workload) or time. The 
following events may be numbered among them: the wear-and-tear of tribological nodes, 
corrosion, materials ageing, etc. [1, 9, 14]. 

The problem with damages to structural components/parts of an engineering object is typical 
of the whole operational phase thereof, although different causes predominate at different 
operational-phase stages. 

Throughout the whole operational phase of an engineering object the structural members 
thereof get degraded. Operating conditions (material effort, deformations, wear-and-tear, etc.), to 
some extent unpredictable — extreme, conditions of externally induced actions/influences that are 
able to change randomly all of them prove decisive in the reduction of the ultimate strength of the 
structure(s) of the engineering object and parts/structural components in the object’s functional 
units/assemblies. 

Damages to an engineering object, in particular while in operation, considerably affect safety, 
but also all rates of engineering availability that feature the operational readiness, i.e. readiness to 
perform tasks. Damages attributable to many and various causes can be classified as a) predictable, 
and b) unpredictable ones. Predictable damages develop at locations easy of access and subjected 
to the diagnosing (with different methods out of a variety of diagnostic techniques). 

Damages to turbine engines apllied in aircraft power plants comprise approximately 25% of 
the total number of disadvantageous events of this kind. 

Damages to engineering objects (such as turbine engines), units/assemblies and structural 
members thereof are classified according to various criteria. From the point of view of damage 
identification and preventive treatment, they are classified for what follows: 

e type of an object (electrical, mechanical, hydraulic, etc. system), 
e physical nature of the damage, 

e causes of damages, 

e ways of restoring the ‘fit-for-use’ state. 


Damages/failures to engineering objects and structural components thereof may also be 
divided into the following groups: 
e primary (independent), 
e secondary (dependant), 
e natural, and 
e forced damages. 


Primary damages are usually effected by: 
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- natural ageing and wear-and-tear processes, fatigue, corrosion, etc., 
- failures after permissible regimes have been exceeded, 

- structural, manufacturing, and material-inherent latent defects, 

- faulty operational and/or servicing practices. 


Primary damages give grounds for statistical inference on the components’ life and reliability 
[6, 7, 8]. Secondary damages are recorded as effects of primary damages. 
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Fig. 1. Structure of turbine engine for maritime applications 


Operational/maintenance practice proves that each and every structural member of a turbine 
engine (Fig. 1) is susceptible to damages/failures which occur for various reasons. The most 
essential causes comprise [2 - 7, 14]: 

- variable aerodynamic loads, 

- increase in loads (in particular, exceeding permissible values), 

- impact loads, 

- vibration of rotating components/parts and those exposed to cyclic changes of position, 
e.g. in the to-and-fro motion, 

- acoustic vibration, 

- heat effects, 

- out-of-standard actions by aircrew members and engineering staff responsible for 
maintenance, 

- weather effects (turbulence, corrosion) and nature-imposed ones (birds), 

- effects of destructive warfare agents, foreign matter, etc. 
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Fig. 2. Quantitative statement of damages/failures to systems and units/assemblies of turbine engines AI-25 


Fig. 2 shows exemplary quantitative statements of damages to systems and units/assemblies of 
the AI-25 turbine engines, recorded while operating a small fleet of aircraft in the five years’ 
period. 

Damages to the engine due to the exceeding of specified values or symptoms of irregular 
engine operation often compels adjustments of engine units, replacement of a damaged unit, 
building the damaged/failed component/part or unit out of the engine. Such operations can be 
carried out either on the engine test bench (for engines after overhaul or withdrawal from service) 
or built in the object where it is operated (sea vessel or aircraft). 


3. Mechanical damages 


Damages of the mechanical type due to foreign-matter collisions contribute to the greatest 
extent to failures to the structures of the engine block, air inlets, and compressor, the fact being 
confirmed by the records. In aviation, they most often occur during the take-off and landing and 
are usually effected by concrete and/or soil chips, and birds. Detection of such damages during 
routine inspections/maintenance usually produces no difficulty, whereas any failure to detect them 
can generate top hazards (Figs 3 and 4) [1 - 3, 9, 11, 13]. 


4. Corrosion damages, tribological and fatigue failures 


Primary causes of such damages/failures are physical and chemical phenomena, and exceeded 
performance parameters (service loads). Theses are as follows: corrosion (Fig. 4), erosion, 
cavitations, tribological wear, fatigue and brittle cracking. Detection of damages/failures effected 
by such phenomena often proves very difficult. They are usually found in the form of the unit’s 
failure to operate, pits that weaken the structure, etc. The engineering diagnosis [1 - 4, 10 - 

11, 13, 15] offers tools and techniques to detect and identify them. Figs 5— 8 show examples of 
such damages/failures. 
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Fig. 3. Exemplary mechanical and heat-induced damages in the form of cracks 
and material defects in the turbine engine structure 


Fig. 4. Fatigue-corrosion-attributable damages to turbine-engine structural components 


5. Description of the occurrence of damages/failures to a set of turbine engines 


The occurrence of damages/failures to engineering systems (turbine engines included) that 
belong to a homogeneous set is featured with the ‘unfit-for-use’ states of structural components or 
units/(sub)assemblies. In the case of reparable systems (engines), any damages/failures are 
removed by means of the restoration process. 
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Fig. 6. Operation-induced III” damage to the turbine of a jet engine [12] 
A — blade ring of the turbine, one blade lacking; B — simulation of destruction of the blade locking piece; 
C- simulation of vibration destructive to the blade locking piece; D — simulation of vibration 
of turbine with one blade lacking 


Fig. 7. Examples of fatigue failures to turbine engine’s compressor blade 
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Fig. 8. A damaged bearing of the turbine engine shaft 


The stochastic occurrence of damages/failures to a set of engines can be described as soon as 
probability distributions of emerging events (damages/failures) within a certain assumed time 
interval of system’s operation. 

Engines belong to the class of reparable engineering systems featured with that the reliability 
R(t) after restoration should not change. Unfortunately, there are such engines for which both the 
number of observed damages/failures and failure rates increase after the major repair/general 
overhaul. The AL aero engine is a good example. This fact can be explained in the following way: 
the design structure proves faulty, or incorrect general-overhaul procedures have been developed 
and recommended by the manufacturer. For such systems, the probability P, of the random 
variable of the number of defects/failures U; in time interval of the system’s operation [0 < qt < 4] 
can be most often described with the following expression: 


T 


po) EL expla a) 


where: A - operational failure rate, 
n - the number of observed damages/failures. 
The expected value of equation (1) is as follows: 
E {U} = 8t =A (2) 
For any structural component/part, engine unit/assembly (x), the probability Px of the random 


variable of the number of defects/failures U, in time interval of the system’s operation [0 < t < 4] 
can be described with the following expression: 


A N 
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of the expected value: 
E {U} =A px (4) 
where: px - conditional probability that a given component, assembly/unit of the x type 
gets failed ifthe engine is damaged, 


m — the number of observed damages/failures. 


For the F set of structural components, units/assemblies in engines in service the expected 
value takes the following form: 
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E {F U} =F A pk (5) 


All the above-mentioned values of parameters and conditional probabilities are arrived at by 
means of suitable statistical surveys under real operating conditions. 
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Abstract 


Gathering reliability data is quite difficult. Potential sources of such data are: reliability guides and handbooks, 
data collecting during operation of technical devices, expert judgments. Practice shows that very often the statistical 
sample we have to our disposal is small and non — homogenous. The goal of the author was to verify, if it is possible to 
draw out useful in practice conclusions about reliability of technical items from small statistical sample and how to do 
it. In the article two potential ways to deal with the problem of small reliability data sample have been presented. First 
way is based on using mathematical statistics methods, the second way uses the fuzzy sets theorem methods. 


Key words: power system, reliability data, statistics, fuzzy sets. 
1. Introduction 


The key problem in making reliability analysis of power systems, is getting credible reliability 
data, which are the entrance data to created mathematical models. The probabilities of suitable 
technical elements failures are on generality these data. Problem in this place appears - where from 
to draw these data. The sources of reliability data for technical structures can be: reliability guides 
and handbooks, data collecting during operation of technical devices, reliability tests, experts' 
opinions. Practice shows that very often the statistical sample we have to our disposal is small and 
non — homogenous. 

The goal of the author was to verify, if it is possible to draw out useful in practice conclusions 
about reliability of technical items from small statistical sample and how to do it. 

At first some assumptions have been made: non repairable items have been considered, the 
statistical distributions of time to failure of items and parameters of those distributions have been 
objects of interest as well as probability of failure in given period of time. 

The two potential ways to deal with the problem of small reliability data sample have been 
considered. First way was to use mathematical statistics methods, the second way, to solve the 
problem, was to use the fuzzy sets theorem methods. 


2. Mathematical statistics methods 


Mathematical statistics methods are based on statistical hypothesis testing. We need to 
remember that there are two possible errors made in a statistical decision process. Error of the first 
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kind known as an a error — which is the error of rejecting a hypothesis when it is actually true. 
Error of the second kind known as a p error — which is the error of failing to reject a hypothesis 
when it is in fact false. 

In practical applications it is very hard to determine the probability value B of the error of the 
second kind. When the statistical sample is small and non — homogenous it is just impossible. That 
is why our care is focused only on the error of the first kind. In advance we impose the probability 
value of the error a, which in reliability analysis applications is mostly equal 5 % and is called the 
level of significance. 

First of all we have to state the relevant hypotheses to be tested. Than we carry out the test of 
significance to reject our hypothesis, or to say that there is no reason to reject the hypothesis. The 
very important thing is to remember that we cannot say the hypothesis is true, in the best case we 
can only establish that there is no reason to reject the hypothesis. 


2.1. Overview of parametric statistics models 


Three typical models, described in handbooks [2, 3], have been studied to answer a question: if 
it is possible to draw out useful in practice conclusions about reliability of technical items from 
small statistical sample using those models. 


Model 1 


Let's make an assumption that our statistical parameter is time to failure of an item. The time to 
failure has the normal distribution N (u, o) with the mean value u and the standard deviation o. 
Moreover the standard deviation value is known for the population of such items. 

The hypothesis about the mean value (Ho: U = po) is tested against the alternative Hypothesis Ha 
(u # po) with the significance level of a test a. 

We can see at once, that such model is not useful for drowing conclusions about reliability of 
technical items if the standard deviatin value o for population is not known. At our considerations 
the o value is not known, so the model 1 is useles. 


Model 2 


Let’s make an assumption that our statistical parameter is time to failure of an item. The time to 
failure has the normal distribution N (u, o) with the mean value u and the standard deviation o. 
Moreover both: the mean value u and the standard deviation value o are not known for the 
population of such items. The hypothesis about the mean value (Ho: U = Ho) is tested against the 
alternative Hypothesis Ha (u # Uo) with the significance level of a test a. 

To verify the hypothesis Ho in such case, we can use the hypothesis testing model based on the 
Student’s t distribution. The o value is estimate than from the sample we have in our disposal. 

The model can be useful for us on one condition: we have to be convinced, that the time to 
failure of technical items has a normal distribution. 


Model 3 


Let’s make an assumption that our statistical parameter is time to failure of an item. The time to 
failure has the unfounded statistical distribution with the mean value u and the standard deviation 
o. Both parameters u and o are not known. The hypothesis about the mean value (Ho: U = po) is 
tested against the alternative Hypothesis Ha (u # u o) with the significance level of a test a. To 
verify the hypothesis Hg in such case, we can use the hypothesis testing model based on the zero — 
one standarised normal distribution N(0,1) . The u and o values are estimate from the sample we 
have in our disposal. 

The model can be used only when we have at our disposal a large statistical sample number. 
The large sample number means at least 30 observations [2] . Some authors say 100 or more 
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observations [3]. Using this model we can estimate the mean value and standard deviation of time 
to failure of item but we cannot to find out the time to failure distribution density function shape. 


Similarly, like it was done in the above models for the mean value, we can also make statistical 
significance tests, using the above models, for the standard deviation values. However there are 
some conditions to choose the model [2]: 

When we have the normal distribution N (u, o) of time to failure, i and o values are not known 
and sample number n < 50 then we choose model 1. 


When we have the normal distribution N (u, 6) of time to failure, i and o values are not known 
and sample number n > 50 than we choose model 2. 


When we have the normal distribution N (u, o) of time to failure, u and o values not known and 
sample number n < 50 than we choose model 3. 


It's worth to notice that we can use those three models for statistical significance tests only if we 
are convinced, that the time to failure of the item has normal distribution. 


2.2. Overview of nonparametric statistics models 


Nonparametric hypothesis testing models are used when we do not know the random variable 
probability density function shape as well as parameters of the function. We have face to such 
situation very often in reliability analysis. However using data set, we have in our disposal, we are 
able to build a histogram of time to failure (which is our random variable) first, and then make 
density estimation based on the data to evaluate theoretical probability density function of time to 
failure. 

Statistical hypothesis testing in this case let us only to find out how much our theoretical model 
(theoretical density function) is adequate to our data distribution in the sample we have. 
Unfortunately, we never can be sure the model will be adequate to whole population of technical 
items being the object of our interest. 

The most popular nonparametric tests in statistics are: the Pearson's chi - square test and the 
Kolmogorov — Smirnov test, in that case used as a nonparametric tests of equality of one- 
dimensional probability distributions used to compare a sample with a reference probability 
distribution. 

The most suitable, to solve the problem given in this article, seems to be the Kolmogorov — 
Smirnov test. The Pearson's chi - square test needs a large sample (at least 100 items). For the 
Kolmogorov — Smirnov test a much smaller sample will be enough. Of course the test is adequate 
for different density functions of the random variable. Moreover the Kolmogorov-Smirnov test is 
the only test usually used to test whether a given distribution function F(x) is the underlying 
probability distribution of F,„(x), the procedure may be inverted to give confidence limits on F(x) 
itself. If one chooses a critical value of the test statistic D, such that P(D, > D,) = a, then a band of 
width +D, around F,(x) will entirely contain F(x) with probability I — a [4]. 

The last feature is very important in reliability analysis, because if the random variable is time 
to failure of the item, then the density function is just the unreliability function. So we are able not 
only to estimate the unreliability function but also to find the confidence limits for estimated 
unreliability function. 

If the density function of time to failure of technical items is normal, then we can replace the 
Kolmogorov — Smirnov test by the Lilliefors test (sample size n < 500) or the Shapiro Wilk test 
(sample size n > 30) or Anderson - Darling test. The last test is more powerful [4]. 
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3. Methods based on fuzzy logic 


At the beginning let’s make an assumption, that we have no statistical data of time to failure of 
technical items. The only way to collect needed reliability data is to use experts judgments. It is 
very difficult for experts to give answers in the form of numerical values, for example: the time to 
failure of the item is 5000 hours or the failure probability per year is 0,003. Much easier is to ask 
them using linguistic values. The problem arise - how to converse the linguistic values into 
numerical values. The fuzzy logic can be helpful to deal with the problem. 

Assume that the object of our interest is the probability of failure of technical item during one year 
period of time. The idea of the author is to test the truthfulness of the sentence given in such form: 
“The item will fail in the one year period of time”. Experts give their opinions shown in Tab. 1. 


Tab. 1. An example of using expert opinions in reliability analysis of technical items 


Tested sentence: 


Possible experts opinions: 


The item will fail during one year period of time. 


(1) Sentence is absolutely truth. 
(2) Sentence is very truth. 

(3) Sentence is truth. 

(4) Sentence is rather truth. 

(5) Sentence is truth or false. 
(6) Sentence is rather false. 

(7) Sentence is false. 

(8) Sentence is very false. 

(9) Sentence is absolutely false. 


Now we can transform expert opinions into fuzzy numbers according to so called “standard 
degrees of truth” proposed by Baldwin in [1] with the membership functions u(v) given in Fig. 1. 
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u(v) u(v) u(v) 


0 0 0 
u(v) u(v) 
1 1 
(5) 


(=N E] 


0 1 Vv 0 LY 0 
Fig. 1. Standard degrees of truth by Baldwin given as the fuzzy numbers 
As the result we receive unreliability value F, for given period of time, as a fuzzy number. For 
instance: 


u(F) u(F) 
1 1 


1 F F 
Unreliability of item for given period of time Unreliability of item for given period of time is 
is equal “about 1”. “very close to 1”. 


Fig. 2. Examples of unreliability values given in the form of fuzzy numbers 
Fuzzy logic is a very good tool when we have in our disposal very imprecise data. In contrast 
with binary logic, where the truth has only two values (1 — true, 0 — false) fuzzy logic variables of 


truth may vary form 0 to 1. It gives us possibility to deal with reasoning that is approximate rather 
than precise. And such situation is very typical in reliability analysis of technical system. 
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Another interesting way to solve the problem of reliability data evaluation, in author’s opinion, 
seems to be the possibility theory. The main idea of possibility theory is to replace probability 
measure with possibility measure and necessity measure. 

If we are not able to evaluate probability value of item failure, we can try to evaluate the 
possibility measure and necessity measure — what is much easier. Let assume A — is event of item 
failure. Between probability, possibility and necessity measure there is relationship given below 
[1]: 

N(A) < P(A) < (A) (1) 
where: 
N(A) — necessity measure of event A, 
P(A) — probability measure of event A, 
II(A) — possibility measure of event A. 


Now let's try to illustrate usage of the possibility and necessity measure in reliability analysis 
on example. Given: A — event of the item failure during one year period of time and 5 statistical 
data in form of time to failure of five items (a, b, c, d, e) presented in Fig. 3. 


to ta tb te ti ta te t 
to—tı — one year period of time, 
tot, —time to failure of item a, 
to-t —time to failure of item b, 
tot. —time to failure of item c, 
to—tg —time to failure of item d, 
to—te —time to failure of item e, 


Fig. 3. Time to failure of items 


Given the above set of data we can evaluate: 
Possibility measure of event A: II(A) = 1, 
Necessity measure of event A: N(A) = 3/5 = 0,6. 


So we can give a statement 1: Probability of item failure during one year period of time can range 


between (0,6 ; 1). 


Imagine now that we have a new information about another (f) item time to failure. 


First situation: the time to failure of a new item is less then one year — then a new necessity 
measure is N’(A) = 4/6 = 0,67 so we have a new statement 2: Probability of item failure during 
one year period of time can range between (0,67 ; 1). Conclusion: statement 2 is not contradicted 
to statement 1. Statement | is still truth. 


Second situation: the time to failure of a new item is more then one year — then a new necessity 
measure is N’’(A) = 3/6 = 0,5 so we have a new statement 3: Probability of item failure during one 
year period of time can range between (0,5 ; 1). Conclusion: statement 3 is contradicted to 
statement 1. Statement 1 is not truth, but statement 1 is more pessimistic than statement 3 what is 
consistent with the rule of the worst case which is often used in reliability analysis. Statement 1 
then is still worthy for us. 
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4. Final remarks 


In the authors opinion, it is possible to draw out useful in practice conclusions about reliability 
of technical items from small statistical sample, using statistical methods, fuzzy logic and 
possibility theory. 

We should always to remember, that received results in reliability analysis are never sure. We 
have to deal with the problem of large uncertainty. Fuzzy logic and probability are the ways of 
expressing uncertainty. 

With the very small data set it is very hard to talk about probability in the classical sense, even 
from statistical point of view. We should rather to use fuzzy numbers, fuzzy probability with the 
help of experts judgments or the possibility theory, as it has been shown in the article. 
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Abstract 


In the article there are presented different operating events leading to the same consequences, i.e. permanent 
engine damage. The analysis has been carried out of diagnostic actions performed by the operators, of the usefulness 
of diagnostic systems existing on these ships as well as analysis of the operating values of the latest diagnostic systems 
in respect of possibility to prevent the described consequences of operating events 


Key words: piston engine of the main ship drive, technical state of engine, diagnostic parameters, diagnosing 
system, operating events, efficient operation of ship engine 


1. Introduction 

Task of the diagnostic system (SDG) is to provide the diagnostic data concerning a ship engine 
(diagnosed system - SDN). SDG and SDN form a diagnostic system. Efficiency of SD depends on 
a possibility to obtain accurate and reliable diagnostic information, which enables a mechanic to 
make accurate and reliable operating decisions. 

In the presented operating events such information was missing, what resulted in considerable 
damages and threat to the crew and ship safety. 

In the article there are presented different operating events leading to the same consequences, 
i.e. permanent engine damage. The analysis has been carried out of diagnostic actions performed 
by the operators, of the usefulness of SD existing on these ships as well as analysis of the 
operating values of the latest SDs in respect of possibility to prevent the described consequences of 
operating events. 


1. Examples of damage to marine engine identification 


1“ Example 

In the four-stroke engine if the ship propulsion plant the crack and tear off of the connecting- 
rod big end bottom cover took place. 

After visual inspection of the engine it has been established that a cause of the failure was 
damage to the bottom cover of a connecting-rod big end. 
> Results of the accident examining board inspection: 

“Cause of a fatigue crack of the cover was repair by pad welding of discontinuity found in 
manufacture process after forging and machining operation. 
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As it results from examination of the damaged part, a cause of failure was the faulty heat 
treatment in manufacturing process. After pad welding no actions were performed in purpose to 
eliminate thermal stresses. The fact was neglected that the above element was not among the 
consumable elements and that it was supposed to serve through a whole period of the engine 
operation, but it was damaged after 4 years despite it had certificates of two recognised 
classification societies. 

The engine elements getting out of the crankcase caused damage to the oil and fuel pipes. 
Before black-out occurred in result of stoppage of the engine coupled with a shaft generator, the 
fuel and oil had been sprinkled around the engine room in result of a pressure generated by the fuel 
and oil pumps driven by electric motors. It should be admitted that there was an unbelievable luck 
that the ignition did not occur and fire did not break out, which could have had disastrous 
consequences”. 

“Consequences. 

In result of the occurred event a crack of piston no. 4 took place, a piston and crank system was 
pushed outside (Fig.1), the crank bearing cracked, camshaft case was torn off the cylinder block 
along 3000 mm length, camshaft was bended, material was torn out in the lower part of cylinder 
sleeve guiding, the sleeve was planished and wedged up with a connecting rod, the crank-pin no. 4 
was deeply damaged, the rotating masses of the system connecting rod were damaged, the fuel 
pipe cracked and heavy fuel mixed with lube oil (loss of approx. 6.5 t of oil), crankcase base 
fastening the fuel pumps no. 2, 3, 4 and 5 was torn out, the pistons, sleeves and connecting-rods of 
a system no. 4 + 5 were crushed. 

Explosion did not cause any victims. Despite the fuel fumes the fire did not break 
out.”(Fig.2)[1] 


Fig.1 Visible connecting rod — broken and pushed outside the crankcase (1) and broken piston (2): head (2a) and 
fragments of piston skirt (2b) [1] 


Indirect results in form of a cost of spare parts, the costs of withdrawing the ship from 
operation, indirect costs related to lack of propulsion. 


> Applied diagnostic system has not raised a warning about exceeding of the limit values.. 
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Fig.2 Fuel spilt around the engine causes fire risk; starboard [1] 
&: p 8i 
znd Exam ple 


The fuel piston ring was broken and after two days of voyage piston was seized and the 
connecting rod was pushed outside of the crankcase. 
> Results of the accident examining board inspection. 

Breaking up of the fuel piston ring resulted in deterioration of the process of lube oil collection 
from the cylinder sleeve bearing surface and damaging of the sleeve bearing surface by the 
damaged ring. Excess of lube oil caused origination of oil carbon deposit on the piston skirt, in the 
ring grooves in the piston crown and, in consequence, the piston rings were stuck with carbon 
deposit (Fig.3). Lockup of the broken oil ring caused spring crack, which was pushed into a space 
between the piston skirt and cylinder sleeve. The spring damaged a hardened sleeve surface layer 
and was pressed into the sleeve. 


Fig.3 Origination of oil carbon deposit on the piston skirt, in the ring grooves in the piston crown [source: Ch.Eng. 
Marceli Stelmaszczuk] 
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Fig.4. The spring set in the material of cylinder sleeve (1), oil ring spring lock with a fragment of broken spring (2) 
[source: Ch.Eng.Marceli Stelmaszczuk] 


In effect of friction forces and released heat the spring material was melted into the piston 
material, what is shown in Fig. 4, 5. Piston seizing caused breaking off the pins fastening a bottom 
cover of a connecting rod big end with the connecting rod and pushing the connecting-rod outside 
the casing, causing the effects analogous to those caused by crack of the bottom cover of a 
connecting rod big end presented in the 1** Example. 

In the last phase of failure a noise could be heard when the connecting-rod with a piston were 
being pushed outside the casing. 

Pressure drop of the recycle stock and high temperature of cooling water at the outlet from 
seized sleeve caused the emergency engine stoppage and ‘black-out’ (the main engine was driving 
a shaft generator). 


Fig.5. Spring from the piston ring pressed into the cylinder sleeve; visible black carbon deposit originated in 
result of incomplete combustion of cylinder oil [source:Ch.Eng. Marceli Stelmaszczuk] 


> Applied diagnostic system has not raised a warning about exceeding of the limit values. 
2. Diagnostic capabilities of the modern SDG 


In recent years the possibilities of engine action monitoring by the technologically advanced 
diagnosing systems have increased. 
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> Electronic analyser of combustion process, which monitors the engine action periodically, 
allows faster detecting of the changes in combustion process parameters 

> Temperature sensors detect lower rises of combustion gases temperature at the outlet of the 
damaged system head, 


Fig.6 Four-stroke engine section with the places marked, where the sensors are mounted — 1) combustion gases 
temperature meter, 2) combustion process analyser, 3) detector of piston rings wear, 4) sensor of cylinder sleeve 
metal temperature [3] 


> The systems like SIPWA-TP allow fast detection of temperature rise of the cylinder sleeves 
metal (Fig.7) 
> Diagnosing systems like MAPEX — PR allow fact detection of piston rings wear (Fig.7) 


Demo Part - LWT Gractic Dsolay-MAPEXPR — DemoPlunt - Sal Cat Load Speed - MAPEX-PR 


Modem for Satcom 


Floppy disk 


> 


Alarm output 


Fig.7 Monitoring of the sleeve metal temperature, a temperature sensor is marked (1), MAPEX — PR system — 
monitoring of the piston ring wear (2) [2] 
> Use of advanced electronics to monitoring of engine action. Data are sent to a ship owner and 
engine manufacturer (Fig 8). 
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Wärtsilä Switzerland 


MAPEX-SM 


MAPEX-SM 


Fig.8. MAPEX — PR programme shows differences between actual measurement and the planned values, e.g. 
measurement of piston rings wear (1), monitoring of diagnostic parameters by the manufacturer and ship owner office 


(2) [2] 


Increase of a number of the diagnostic data, development of a method of their analysing and 
processing, are leading to multicriteriality of assessment of the engine technical condition [4,5,6]. 
In result, probability value of the relation between the diagnostic parameters and engine technical 


conditions represented by them increases (Fig.9) [7] 
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Fig.9 Probabilities of relation occurrence between the diagnostic parameters and engine technical conditions 
represented by them (1) [7] 


Conclusions 
During engine operation the events occur, which can not be detected in actual time and the 


1. 
consequences of their occurrence can not be limited (E. 1) 
2. The event described as E.1 should be recognised as a human error, resulting from lack of 


experience and knowledge, and lack of the basic safety principles. 
3. Electronic combustion process analyser allows faster detecting of the change in engine 


running parameters, forcing a mechanic to analyse the reasons of these changes. 
4. Use of sensors of the sleeve metal temperature allows detecting the temperature rise faster 


than with indirect method through measurement of temperature of the sleeves cooling water 
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5. Application of a system determining an extent of the piston rings wear would probably allow 
early detection of the oil ring damage (E.2) 

6. On basis of the quoted examples it is obvious how important it is to use safety procedures at 
every stage of engine manufacturing and operating 

7. Despite the development of diagnostic systems (SD), both SDG and SDN, a very important 
factor in making operating decisions is an accurate and reliable diagnosis worked out by a 
mechanic on basis of information from SD (as good as knowledge of a design engineer and 
manufacturer) and resulting from his knowledge and experience 

8. Development of both SDN and SDG (including application of the respective mathematical 
models) increases probability value of relation between the diagnostic parameters and technical 
engine conditions represented by them 
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Abstract 


For the failure moments and the time between failures of the marine power plants fuel oil systems, an attempt of 
defining the failure distributions has been made. The analysis was based on the observations of the failure of marine 
power plants fuel oil systems elements. Due to the character of the statistical data, the analysis deals only with 
relatively simple distributions, most frequently used in the reliability theory. Failures to the marine power plant 
systems of 10 ships owned by the Polish Steamship Company of Szczecin was the subject of a statistical data analysis. 
All the ships differed in respect to their place and time of construction as well as their technical parameters. Data on 
their failures refer to the fuel oil system. The data on marine power plants failures were collected in similar 
conditions, that is, they were supplied by an engine crew member working in the marine power plant. The data on the 
failures of particular marine power plant systems were obtained accordingly to the test schedule [N, W, T], which 
means that N renewable objects were the subject of the test within the time T. Since the recovery time of the damaged 
system appeared negligibly short, when compared to the time of the test, it was assumed that consecutive recoveries 


overlap the failure moments. The statistical analysis dealt with moments t; St, S...<t,, of the particular systems’ 
consecutive failures and the length of time intervals T, between the objects’ consecutive failures. In order to carry 
out the distribution estimation, the statistical package STATISTICA was used. 


Keywords: marine power plant, fuel oil system, failure distributions 
1. Introduction 


There has been carried out a statistical analysis of the data concerning marine power plants 
systems failures on ten ships owned by the Polish Steamship Company of Szczecin. The ships 
were given symbols from S1 to S10. All the ships differed in respect to their place and time of 
construction as well as their technical parameters. The data on failures cover the following six 
marine power plant systems: lubricating oil system - LOS, sea water cooling system - SWCS, 
fresh water cooling system - FWCS, fuel oil system - FOS, compressed air system - CAS and 
steam system - SS [2, 3]. 

The data on marine power plants failures were collected in similar conditions, that is, they were 
supplied by an engine crew member working in the marine power plant. 

The data on the failures of particular marine power plant systems were obtained accordingly to 
the test schedule [V, W, 7], which means that N renewable objects were the subject of the test 
within the time 7. Since the recovery time of the damaged system appeared negligibly short, when 
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compared to the time of the test, it was assumed that consecutive recoveries overlap the failure 
moments. 
The statistical analysis dealt with moments 4 <t, <..<t, of the particular systems’ 
consecutive failures and the length of time intervals 7, between the objects’ consecutive failures. 
The following has been assumed [2, 3]: 
1) on ships S1 - S5 time is measured from the moment of the first failure repair; 
2) for ships S6 - S10 time intervals 7, between the consecutive failures do not include the time 
between the beginning of sea voyage and the first failure; 
3) time between the last failure and the end of the observation, that is after 180 days, was taken 
into consideration. 
Total numbers of failures in particular systems of ten tested ships’ marine power plants have 
been presented in table 1. 


Table 1. Total number of failures to particular marine power plant systems of all the tested ships 


Name of the system Total number of failures 


Lubricating oil system - LOS 
Sea water cooling system - SWCS 
Fresh water cooling system - FWCS 
Fuel oil system - FOS 


Steam system - SS 
Sum 


30 25 25 29 


On the basis of the table, most failures, that is 30% of all of them, can be pointed out in the fuel 
oil system. It is to be noticed, that among all the tested ships, the prevailing number of failures 
occurred in case of three of them, that is, S6, S7 and S8. 

On the basis of the obtained field data on the system failures, for particular systems, there was 
built a mathematical model of distribution of time to failure and time between failures. 

Due to the type of the statistical data, the most favored distributions appeared relatively simple 
ones, most frequently applied to the reliability theory. Thus, the possibility of building a model 
based on the following distributions was sequentially examined: 

- the exponential distribution, with probability density function 


f(A =e” dla t20, (1) 
- the Weibull distribution , with probability density function 
f(t)= ope wl dla t>0, (2) 
- the log-normal distribution, with probability density function 
1 {incu Pp 
(t) =e * dla t>0, (3) 
JO tov27 
- the gamma distribution, with probability density function 
p" 1e! e 
t) = —— e dla ¢>0. 4 
85 a) (4) 


The results of the Kruskal-Wallis test, aided by the program STATISTICA 6.0, for the statistical 
data of fuel oil system failure moments ¢, concerning all the 10 tested ships, at the assumed 
significance level a=0,05, turn out to come from one general population (for space reasons the 
analysis has not been enclosed). However, when considering the statistical data concerning the 
time between the fuel oil system failures (z„) of the 10 analyzed transport ships at the significance 
level a=0,05, the hypothesis that they come from one general population needs to be rejected. 
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Thus, multiple comparisons of rank means (post-hoc tests) for the data of the time between failures 
(tn) in reference to each pair of the 10 ships, using the program STATISTICA 6.0 have been 
performed. For each comparison, there have been computed ‘z’ values and ‘p’ values. On the basis 
of the obtained results it has been concluded that data from ships S3, S5 and S7 are significantly 
different from each other. Therefore, it was assumed that the data could belong to various 
populations. However, to the data from S1, S2, S4, S6, S8, S9 and S10 there was applied the 
Kruskal-Wallis test, whose results point out that there are no bases for rejecting the hypothesis that 
they come from one general population at the significance level a=0,05. A broader presentation of 
the above analyses goes far beyond the frames of the paper and would require a separate 
publication, what has initially been done in [1]. 

On the basis of the obtained results, data on fuel oil system failures of S1, S2, S4, S6, S8, S9, 
S10 concerning both consecutive failure moments (£,) and time between failures time (z, ) need to 
be treated as coming from one general population. Therefore, an attempt of estimation of the 
obtained empirical distributions by means of the relatively simple theoretical exponential, gamma, 
log-normal and Weibull distributions, most frequently applied to the reliability theory. In case the 
exponential, gamma, or log-normal distributions have been accepted as the models, the program 
STATISTICA 6.0 enables adjustment and application of tests of goodness of fit. But in reference to 
the Weibull distribution, it only allows for the estimation of the distribution parameters by the 
Maximum Likelihood test. Such an analysis has been carried out by means of the statistical 
package STATGRAPHICS, enclosed in [2]. 


2. Tests of goodness of fit 


Tests of goodness of fit allow for the verification of the null hypothesis that an analyzed 
random variable has a distribution which belongs to the appropriate distribution group. The oldest 
test of goodness of fit appears Pearson 's chi-square test, first presented in 1900 [4]. The test can 
not be applied to small size samples (n<50) [4, 5]. In case of small size samples the Kolmogorov- 
Smirnov test is used. 

Program STATISTICA [4, 5] performs the chi-square test by default on the basis of the number 
of the observed and expected frequencies. Categories containing expected frequencies lower than 5 
are grouped to make classes of higher frequency. The number of degrees of freedom for chi-square 
test statistics is computed in the following way: 
df = number of categories - number of parameters - 1, 
where: 

- number of categories refers to the number of classes in the frequencies table, where the expected 
frequencies are higher than 5; 
- number of parameters refers to the parameters of the appropriate theoretical distribution. 

If the chi-square test outcome is marked with *df adjusted”, it means that to compute chi- 
square test statistic, the program joined the categories with expected frequencies lower than 5. In 
particular, the neighboring categories are joined until their frequency becomes higher than 5. 

If the test demonstrates statistical significance (that is, p-value is lower than the assumed 
significance level, usually 0,05), the hypothesis, that the observed data are liable to a specific 
distribution, is rejected [4, 5]. 

The Kolmogorov-Smirnov test can be computed, or not, depending upon the settings on the 
chart Options. In case of two samples, the Kolmogorov-Smirnov test is the test of statistical 
significance differences between them. Just like the test for one sample, the test statistic compares 
cumulative distribution functions, in this case it deals with cumulative distribution functions of 
two samples (e.g. the observed and the simulated values). Big difference between the cumulative 
distribution functions shows that the data come from two different populations. The test may turn 
out useful when building a model for assessment whether the expected results (based on the 
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simulation input data) differ from the observed ones. A significant difference between the 
expected and observed outcomes usually reveals the model’s insufficiency for demonstrating 
dependencies between the input and output data [4, 5]. 

In the Kolmogorov-Smirnov test group there may be set options for Kolmogorov-Smirnov test. 
Calculations can be done on the basis of categorized data by clicking the button Yes (categorized) - 
(computing is faster) or using raw data (the procedure is slower; button Yes (continuous) must be 
pressed to make use of this method). In the first case the program calculates the value of D-max 
statistic on the basis of grouped data, whereas in the second case, at every point it performs 
calculations of cumulative expected frequencies based on the ordered data. Kolmogorov D-max 
Statistic is the biggest of the absolute differences between the observed and expected cumulative 
values of the distribution. If the Kolmogorov-Smirnov test shows statistical significance (that is, 
the p —value appears lower than the assumed significance level), we may reject the hypothesis 
assuming that the observed data are subject to the hypothetical distribution. 

Using the STATISTICA 8.0 program (chart Distribution fitting), the verifications of null 
hypothesis of fitting the distribution of the failure moments (¢,,) and the time between failures (t,) 
in marine power plants fuel oil systems with the following theoretical distributions: exponential, 
gamma and log-normal were carried out; for this purpose Kolmogorov-Smirnov (continuous) test 
was applied. For reasons of space the STATISTICA 8.0 spreadsheets with the observed and 
expected frequencies and the outcome of the tests of goodness of fit data with theoretical 
distributions fit have not been included. The results of the tests of goodness of fit have been shown 
on histograms and cumulative histograms of marine power plants fuel oil systems failure moments 
(ta) and the time between failures (t„) with probability density functions and cumulative 
distribution functions of the discussed theoretical distributions. 


3. Results of the tests of goodness of fit of the consecutive failure moments 


The performance of the Kolmogorov-Smirnov test for the exponential distribution, has resulted 
in its following parameters: number of categories: 17; lower limit: 0; upper limit: 170; lambda: 
0,0139; observed mean: 71,7130; observed variance: 1907,3664. 

Histogram and cumulative histogram of the consecutive marine power plants fuel oil systems 
failure moments with the probability density function of the exponential distribution and the 
cumulative distribution function of the exponential distribution and with the Kolmogorov-Smirnov 
test results inserted have been shown in fig. 1. 


variable:t, ; the exponential distribution variable:t, ; the exponential distributior 
d Kolmogorov-Smimov 0,18146, p < 0,01 d Kolmogorov-Smimov 0,18146, p < 0,01 
chi-square test = 39,71150, df = 9 (adjusted), p = 0,00001 chi-kwadrat test = 39,71150, df = 9 (adjusted), p = 0,00001 
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Fig. 1. Histogram (a) and cumulative histogram (b) of the consecutive failure moments of marine power plants 
fuel oil systems with probability density function of the exponential distribution (full line) (a) and cumulative 
distribution function of the exponential distribution (full line) (b) and with the results of the Kolmogorov-Smirnov 
test inserted 
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Relative frequency (%) 


There are no bases for accepting the null hypothesis that the empirical distribution of the 
marine power plants fuel systems consecutive failure moments £, complies with the exponential 
distribution. 

The Kolmogorov-Smirnov test applied to the gamma distribution has resulted in its following 
parameters: number of categories: 17; lower limit: 0; upper limit: 170; scale parameter: 64,6862, 
shape parameter: 3,2067; observed mean: 71,7130; observed variance — 1907,3664. 

Histogram and cumulative histogram of the consecutive marine power plants fuel oil systems 
failure moments with the probability density function of the exponential distribution and the 
cumulative distribution function of the exponential distribution and with the Kolmogorov-Smirnov 
test results, inserted, have been shown in fig. 2. 

variable: t, ; the gamma distribution variable: t, ; the gamma distribution 


d Kolmogorov-Smirnov 0,17071, p < 0,05 d Kolmogorov-Śmirnov 0,17071, p < 0,05 
chi-kwadrat test = 24,38039, df = 8 (adjusted), chi-kwadrat test = 24,38039, df = 8 (adjusted), p 


Relative frequency (%) 


Time [24 h] Time! [24h] 


a) b) 
Fig. 2. Histogram (a) and cumulative histogram (b) of the consecutive failure moments of marine power plants fuel oil 
systems with probability density function of the gamma distribution (full line) (a) and cumulative distribution function 
of the gamma distribution (full line) (b) and with the results of the Kolmogorov-Smirnov test inserted 


There are no bases for accepting the null hypothesis assuming that the empirical distribution of 
the marine power plants fuel oil systems consecutive failure moments ¢,, complies with the gamma 
distribution. 

The Kolmogorov-Smirnov test applied to the log-normal distribution has resulted in its 
following parameters: number of categories:17; lower limit: 0; upper limit: 170; mean (M): 
3,7573; variance: 3,20669; observed mean: 71,7130, observed variance: 1907,3664. 

Histogram and cumulative histogram of the consecutive marine power plants fuel oil systems 
failure moments with the probability density function of the log-normal distribution and the 
cumulative distribution function of the log-normal distribution and with the Kolmogorov-Smirnov 
test results, inserted, have been shown in fig. 3. 
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variable: t, ; the log-normal distribution 
d Kolmogorov-Śmirnov 0,24083, p < 0,01 
chi-kwadrat test = 102,39366, df = 6 (adjusted), p = 0,00000 
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variable: t, ; the log-normal distribution 
d Kolmogorov-Smimov 0,24083, p < 0,01 
chi-kwadrat test = 102,39366, df = 6 (adjusted), p = 0,00000 
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Fig. 3. Histogram (a) and cumulative histogram (b) of the consecutive failure moments of marine power plants fuel oil 
systems with probability density function of the log-normal distribution (full line) (a) and cumulative distribution 
function of the log-normal distribution (full line) (b) and with the results of the Kolmogorov-Smirnov test inserted 


There are no bases for accepting the null hypothesis that the empirical distribution of the 
marine power plants fuel systems consecutive failure moments ¢, complies with the log-normal 
distribution. 

To define the goodness of fit of the obtained empirical distribution of the marine power plants 
fuel oil systems’ consecutive failure moments ¢, with the Weibull distribution, a graphical method 
(the window Graphs 2D>Histograms 2D of the program STATISTICA 8.0) have been used. 
Parameters of the Weibull distribution a, p have been estimated by the Maximum Likelihood test. 
The obtained parameters for the distribution are as follows: number of categories: 17; lower limit: 
0; upper limit: 170; a=1,2665; 6=75,5095. 

Histogram and cumulative histogram of the consecutive marine power plants fuel oil systems 
failure moments with the probability density function of the Weibull distribution and the 
cumulative distribution function of the Weibull distribution have been shown in fig. 4. 
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Fig. 4. Histogram (a) and cumulative histogram (b) of the consecutive failure moments of marine power plants fuel oil 
systems with probability density function of the Weibull distribution (full line) (a) and cumulative distribution function 


of the Weibull distribution (full line) (b) with the distribution parameters a=1,2665, B=75,5095 


4. Results of the tests of goodness of fit of the time between failures 


The Kolmogorov-Smirnov test applied to the exponential distribution has resulted in its 
following parameters: number of categories: 13; lower limit: 0; upper limit: 130; lambda: 0,0704; 


observed mean: 14,1979, observed variance: 247,5795. 
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Histogram and cumulative histogram of the time between the failures of marine power plants 
fuel oil systems with the probability density function of the exponential distribution and the 
cumulative distribution function of the exponential distribution and with the Kolmogorov-Smirnov 
test results, inserted, have been shown in fig. 5. 


variable: 7 ; the exponential distribution variable: r, ; the exponential distribution 
d Kolmogorov-Smirnov 0,07620, p = n.i d Kolmogorov-Smimov 0,07620, p = n.i 
chi-kwadrat test = 1,72240, df = 3 (adjusted), p = 0,63197 chi-kwadrat test = 1,72240, df = 3 (adjusted), p = 0,63197 
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Fig. 5. Histogram (a) and cumulative histogram (b) of the time between failures of marine power plants fuel oil 
systems with probability density function of the exponential distribution (full line) (a) and cumulative distribution 
function of the exponential distribution (full line) (b) and with the results of the Kolmogorov-Smirnov test inserted 


There are no bases for rejecting the null hypothesis that the empirical distribution of the time 
between failures t, of the marine power plants fuel oil systems complies with the exponential 
distribution. 

The Kolmogorov-Smirnov test applied to the gamma distribution has resulted in its following 
parameters: number of categories: 13; lower limit: 0; upper limit: 130; scale parameter: 13,8525; 
shape parameter: 1,6149; observed mean: 14,1979; observed variance: 247,5795. 

Histogram and cumulative histogram of the time between the failures of marine power plants 
fuel oil systems with the probability density function of the gamma distribution and the cumulative 
distribution function of the gamma distribution and with the Kolmogorov-Smirnov test results, 
inserted, have been shown in fig. 6. 


variable: r, ; the gamma distribution variable: r, ; the gamma distribution 
d Kolmogorov-Smirnov 0,07907, p = n.i. d Kolmogorov-Smirnov 0,07907, p = n.i. 
chi-kwadrat test = 1,82516, df = 2 (adjusted), p = 0,40149 chi-kwadrat test = 1,82516, df = 2 (adjusted), p = 0,40149 
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Fig. 6. Histogram (a) and cumulative histogram (b) of the time between failures of marine power plants fuel oil 
systems with probability density function of the gamma distribution (full line) (a) and cumulative distribution function 
of the gamma distribution (full line) (b) and with the results of the Kolmogorov-Smirnov test inserted 


There are no bases for rejecting the null hypothesis that the empirical distribution of the time 
between failures z, of the marine power plant fuel system complies with the gamma distribution. 

The Kolmogorov-Smirnov test applied to the log-normal distribution has resulted in its 
following parameters: number of categories: 13; lower limit: 0; upper limit: 130; mean (M): 
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2,0908; variance: 1,6149; observed mean: 14,1979; observed variance: 247,5795. 

Histogram and cumulative histogram of the time between the failures of marine power plants 
fuel oil systems with the probability density function of the log-normal distribution and the 
cumulative distribution function of the log-normal distribution and with the Kolmogorov-Smirnov 
test results, inserted, have been shown in fig. 7. 
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chi-kwadrat test = 7,82505, df = 2 (adjusted), p = 0,01999 


variable: r, ; the log-normal distribution 
d Kolmogorov-Smimov 0,09969, p = n.i 
chi-kwadrat test = 7,82505, df = 2 (adjusted), p = 0,01999 
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Fig. 7. Histogram (a) and cumulative histogram (b) of the time between failures of marine power plants fuel oil 
systems with probability density function of the log-normal distribution (full line) (a) and cumulative distribution 
Junction of the log-normal distribution (full line) (b) and with the results of the Kolmogorov-Smirnov test inserted 


There are no bases for rejecting the null hypothesis that the empirical distribution of the time 
between failures (t„) of the marine power plants fuel oil systems complies with the log-normal 
distribution. 

To define the goodness of fit of the obtained empirical distribution of the time between failures 
in the marine power plants fuel oil systems (z„) with the Weibull distribution, a graphical method 
(the window Graphs 2D>Histograms 2D of the program STATISTICA 8.0) have been used. 
Parameters of the Weibull distribution a, p have been estimated by the Maximum Likelihood test. 
The obtained parameters for the distribution are as follows: number of categories: 13; lower limit: 
0; upper limit: 130; a=0,9928; 6=14,1522. 

Histogram and cumulative histogram of the time between failures (z„) of marine power plants 
fuel oil systems and the probability density function of the Weibull distribution and the cumulative 
distribution function of the Weibull distribution have been shown in fig. 8. 
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Fig. 8. Histogram (a) and cumulative histogram (b) of the time between failures of marine power plants fuel oil 
systems with probability density function of the Weibull distribution (full line) (a) and cumulative distribution function 
of the Weibull distribution (full line) (b) with the distribution parameters: a=0,9928, B=14,1522 
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5. Final remarks 


In case of the consecutive failure moments of the marine power plants fuel oil systems, non of 
the discussed theoretical distributions turns out a sufficiently appropriate model for the obtained 
empirical distribution. 

However, for the time between failures of the marine power plants fuel oil systems, there are 
no bases for rejecting the null hypothesis assuming the obtained empirical distribution comply 
with each of the considered theoretical distributions, that is, the exponential, gamma, log-normal 
and Weibull. 

Program STATISTICA 8.0 enabled the fit of the theoretical distribution and carrying out the 
tests of goodness of fit in case of the assumption that exponential, gamma and log-normal 
distributions become the models. But in reference to the Weibull distribution, due to the Maximum 
Likelihood test, only the distribution parameters could be estimated. The analysis carried out by 
means of a statistical package STATGRAPHICS, presented in [2], in reference to the Weibull 
distribution, brought slightly different, more unique results; but it was made by a much older 
statistical package than STATISTICA 8.0, and for the space reason it was not possible to discuss in 
details the comparison of the test results obtained by means of both statistical packages. 
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Abstract 


The paper present the results of emission tests of city buses (Euro 5 hybrid and Euro 4 conventional diesel), under 
real traffic conditions with the use of a portable emission testing system. The tests have been carried out in city traffic 
conditions on road portions of several kilometers each; The test results contain information on the vehicle emission 
level in operation and pertain to real road conditions. 


Keywords: emissions, road tests, diesel engines, hybrid drivetrains 
1. Introduction 


Contemporary vehicle manufacturers focus not only on the comfort and safety of their products 
but also on the issues related to the fuel economy [3, 18] and the emission level [1]. Owing to a high 
commitment of the vehicle manufacturers recent years saw a reduction of vehicle emission level 
from Euro 1 to Euro 4, which reduces the pollution by more than 90% [1]. 

European emission regulations for heavy duty vehicles set forth in Directive 1999/96/WE are 
commonly known as Euro 1...5 standards. Since October 2005 all newly homologated vehicles and 
since October 2006 all newly registered commercial vehicles (Heavy Duty Diesel, including buses) 
have had to comply with the Euro 4 standard. More stringent standard - Euro 5 will come into 
force in October 2008 (homologation) and October 2009 (registration) respectively. Additionally, 
the Directive 2005/55/WE adopted in 2005 introduced the EEV standard (Enhanced Environmen- 
tally Friendly Vehicle) for vehicles of particularly low emission level. The purpose of this directive 
is to replace the previous directives 88/77/EWG, 96/1/WE, 1999/96/WE and 2001/27/WE through 
a unification and consolidation of the regulations in a single act. Directive 2001/27/WE came into 
force on 9.11.2006. As the Euro 4 standard was introduced, the manufacturers of all powertrains of 
heavy duty vehicles including city buses had to fit them with OBD (On-Board Diagnostic) emis- 
sion monitoring system [2, 4, 10, 14]. The implementation of the Euro 4 and Euro 5 standards 
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forced the manufacturers to come up with new solutions such as EGR (exhaust gas recirculation) 
or SCR (selective catalytic reduction), the latter requiring an on-board supply of a solution of car- 
bamide (AdBlue). 

Currently, more attention is drawn to the measurement of the emissions under variable operat- 
ing conditions, particularly regarding heavy duty vehicles. Emission testing in road conditions is in 
higher demand than stationary driving cycles. On-road emission testing became possible owing to 
a rapid advancement of the measuring techniques that came in recent years [5, 6, 8, 13, 17]. The 
said advancement also aimed at the measurement of extremely small pollutant concentration in the 
exhaust gas [7, 9, 15]. The Institute of Combustion Engines and Transport with a portable emis- 
sion testing system carried out a series of on-road emission tests of two city buses with two differ- 
ent powertrains: hybrid and conventional diesel. 


2. Emission testing methodology 


The on-road emission tests were carried out in city traffic in Poznan (Fig. 1). The tests were 
performed on the city main streets in the afternoon, when the traffic was moderate. Such condi- 
tions were selected in order to ensure the highest possible use of the engine work field (regular bus 
operation was simulated including 10 second stops at the bus stops). The object of the tests were 
two buses manufactured by Solaris: one unit was fitted with a hybrid engine (Hybrid H18) and the 
other one — with a conventional diesel engine (U18 Aarhus) — characteristics — Tab. 1. The buses 
were selected based on the similarity of service routes (number of passengers) and at the same 
time in such a way that there was a possibility of comparing their functionality and ecology under 
real conditions. 
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Fig. 1. A road portion (routes marked) used for the measurement of the buses emission level 


Tab. 1. Technical data of the tested buses 


Parameter Conventional diesel bus Hybrid Bus (diesel + electric) 
1 2 3 

Engine type DAF PR 265 Cummins ISB 250 
Displacement [em'] 9200 6700 
Emission standard Euro 4 Euro 5 
Transmission VOITH DIWA 86 4.5 ALLISON Ep50 
Vehicle weight [kg] 16,700 17,800 
Vehicle weight full load [kg] 21,200 22,000 
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For the measurement of the exhaust emissions the authors used a portable SEMTECH DS ana- 
lyzer [9, 16] by SENSORS LTD (Fig. 2). The analyzer measured the fuel consumption and the 
exhaust emissions (Tab. 2) at the same time recording the exhaust mass flow. The exhaust gases 
entering the analyzer through a sensor maintaining the temperature of 191°C (Fig. 3) were filtered 
for particulate matter (diesel) and then the concentration of hydrocarbons was measured in a FID 
(flame ionization detector). Next, the exhaust gases were cooled down to a temperature of 4°C and 
the measurement of the concentration of NOx took place (non-dispersive method with the use of 
ultraviolet radiation that enabled the measurement of both nitrogen oxide and nitrogen dioxide), 
carbon monoxide, carbon dioxide (non-dispersive method with the use of infrared radiation) and 
oxygen (electrochemical analyzer). The data could be directly transferred to the analyzer central 
unit from the vehicle diagnostic system and the GPS. 


Fig. 2. View of the testing device (Semtech DS for exhaust emissions tests) fitted in the bus 


Tab. 2. Characteristics of the portable exhaust gas analyzer SEMTECH DS 


Parameter Method of measurement Accuracy 
1 2 3 
1. Concentration of 
CO NDIR — non-dispersive (infrared), range 0—10% 43% 
HC FID — flame ionization, range 0—10,000 ppm +2.5% 
NOx = (NO + NO») NDUV — non-dispersive (ultraviolet), range 0-3000 ppm +3% 
CO» NDIR — non-dispersive (infrared), range 0—20% +3% 
O2 electrochemical, range 0-20% +1% 
Sampling frequency 1—4 Hz 
2. Exhaust gas flow Mass flow intensity +2,5% 
Tmax Up to 700°C +1% 
of the range 
3. Pre-heating time 15 min 
4. Response time Too < 1s 
5. Supported diagnostic systems | SAE J1850/SAE J1979 (LDV) 
SAE J1708/SAE J1587 (HDV) 
CAN SAE J1939/J2284 (HDV) 
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Fig. 3. Schematics of the portable analyzer SEMTECH DS — exhaust gas flow (==) 
and electrical connections (- - -) shown 


3. Tests results and analysis 


With the use of the portable system the emissions of CO, HC, NOx, CO» with 1 second interval 
were measured as well as the changes in the engine speed and torque— the parameters taken from 
the vehicle OBD (CAN SAE J1939) and used to calculate the emissions related to the unit energy 
of the engine. Example data recorded during the cruise are shown in Fig. 4. The ECU of the hybrid 
vehicle limited the torque of the diesel engine when the vehicle began to move — while accelerat- 
ing up to the speed of 5 km/h the electric motor was used and the diesel engine remained idle 
(lower toxic emissions). 


Fig. 4. Example measurement recording (CO, emission overlaid on the vehicle route during the emission tests) 


The calculations of the time density characteristics for the engine parameters in the urban traf- 
fic revealed certain dependencies that characterize the share of given parameters of engine opera- 
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tion in the total cruise time (Fig. 5). For the hybrid bus a large share of the operating time falls 
within the engine operation range of 700-900 rpm and engine load of approximately 10%. For the 
bus with the conventional powertrain, over 2000 s (approx. 20% of the engine operating time, total 
engine operating time amounts to 10,200 s) of the engine operating time was 600 rpm and the en- 
gine load approximately 10%. The ranges of engine parameters used in the city traffic were: idle 
speed (for both of the tested buses) and the external characteristics operating range (maximum 


load for a given engine speed). At the same time, for the hybrid bus a higher share of part loads is 
more visible. 
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Fig. 5. Characteristics of the operating time density: a) hybrid, b) conventional diesel 


The highest intensity of the CO, emission in the tested city traffic conditions we can observe 
for medium engine speeds and for the maximum engine torque (hybrid vehicle, Fig. 6a) and for the 
conventional diesel engine the CO» emission (fuel consumption) increases proportionally to the 
engine speed and engine load (Fig. 6b). The CO» emissions for the hybrid vehicle are approximate- 
ly 20-30% lower than the analogical values for a conventional diesel vehicle. This could indicate a 
lower fuel consumption; the total fuel consumption is calculated with the use of this data as well as 
the share of the engine operation in given engine speed ranges and loads. 
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Fig. 6. Time density characteristics of the CO» emission: a) hybrid, b) conventional diesel 


The maximum intensity of CO emission, given in milligrams per second, falls within the range 
of maximum engine loads and medium engine speeds (hybrid vehicle, Fig. 7a) and for the conven- 
tional vehicle this range is approximately 4—5 times wider (maximum for high loads and medium 
engine speeds, Fig. 7b). This mainly results from the differences in the engine capacity and the fact 
that for the hybrid vehicle medium loads are more frequently the case as compared to the conven- 
tional diesel engine, which operates at high loads for most of the time. 
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Fig. 7. Time density characteristics of the CO emission: a) hybrid, b) conventional diesel 


The maximum values of the emission of hydrocarbons for both of the powertrains is similar 
(approximately 25—35 mg/s), yet, for the hybrid vehicle only the maximum values reach 25 mg/s 
(Fig. 8a); in the range of medium engine speeds and engine loads the values do not exceed 20 mg/s 
(Fig. 8b). For the conventional diesel vehicle there is a nearly linear dependence between the HC 
emission and the engine speed and load (the maximum is observed for the maximum values of the 


engine speed and load). 
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Fig. 8. Time density characteristics of the hydrocarbons emission: a) hybrid, b) conventional diesel 


The area of elevated NOx emission in a hybrid vehicle (Fig. 9a) falls within the range of low 
engine speeds in the whole range of engine loads and the maximum engine speed and high engine 
load. This could be the effect of the application of selective catalytic reduction in this vehicle, 
where, at high exhaust gas temperatures the reduction of NOx occurs. For the vehicle with the 
conventional powertrain (no SCR) a linear growth of NOx emission is observed as the engine 


speed and engine load increase (Fig. 9b). 
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Fig. 9. Time density characteristics of the NOx emission: a) hybrid, b) conventional diesel 
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4. Comparison of the emission level of the tested engines 


The emission intensity under real conditions was calculated with the help of the time density 
characteristics of the engine operation (u, wmz) and the characteristics of the emission intensity for 
aj toxic compound e;(n, M/M.): 


E seal => > Waga -e,(n,M/M.)}, (1) 


n M/M, 


where: 


Erea, j — emission of a j toxic compound under real conditions, 
n — engine speed, 
M/M, — relative engine load. 


Knowing the effective power in each range (determined by the engine speed and load), the unit 
emissions of the toxic compounds were compared, related to the engine power during the whole 
test cycle. The following values were obtained for the hybrid bus: HC — 0.193 g/(kW-h), CO — 
3.981 g/(kW*h), NOx — 2.711 g/(kW*h) and CO) — 570.2 g/(kW-h). For the conventional diesel 
bus: HC — 0.282 g/(kW-h), CO — 9.130 g/(kW-h), NOx — 3.128 g/(kW-h) and CO2 — 763.3 
g/(kW-h). The above, when compared, discloses lower values for the hybrid vehicle: CO by 
56.4%, HC by 31.8%, NOx by 13.3% and CO» by 25.3% as compared to the conventional power- 
train (Fig. 10). 
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Fig. 10. Relative values of the unit emissions of toxic compounds during the on-road tests 
(hybrid and conventional diesel) 


5. Engine emission level indexes 


From the information presented in the paper such as the characteristics of the share of engine 
operation at given engine speeds and loads as well as the characteristics of emission intensity we 
can obtain the factor of multiplication of the increase (decrease) of the emissions under real traffic 
conditions as opposed to the values obtained during the homologation test. The emission index 
(for a given toxic compound) has been defined as follows [12]: 

k = E peal. j i (2) 
* E 


cycle(norm), j 
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where: 


j — toxic compound for which the emission index has been determined, 

Eyeai,j — emission intensity obtained under real conditions ([g/(kW-h)]), 

Ecycle(norm),j — emission intensity obtained in the ESC or ETC test ([g/(kW-h)]) or the boundary val- 
ues adopted as permissible for a given emission standard. 


The knowledge of the actual emission and the test emission (or the one compliant with the 
standard) may serve to determine the emission indexes of the toxic compounds of a given vehicle. 
If there is no information on the engine toxic emissions in the ESC or ETC test, we can adopt the 
permissible values according to the Euro emission standard which is binding for a given vehicle. 

From the data included in Tab. 3 the emission indexes were obtained for each vehicle as they 
complied with different emission standards, hence the actual emission of a hybrid bus was com- 
pared with the emission values set out in Euro 5 (ETC test as the vehicle was fitted with an after- 
treatment system) and the actual emission of the conventional diesel bus was compared with the 
emission values set out in Euro 4 and the ESC test as the engine was not fitted with aftertreatment 
systems other than Oxicat. Because the standard does not provide for the CO» emission level the 
unit CO» emission was converted into unit fuel consumption and then the mean overall efficiency 
of the engine was estimated in the whole road test 


(Mo Hybrid = 0.48, No Diesel = 0.36). 


Tab. 3. Emission level under real traffic conditions and emission indexes obtained in the tests for the hybrid and con- 
ventional vehicle 


Parameter co HC NOx CO, 
1 2 3 + 5 
Actual emission (hybrid) [g/(kW-h)] 3.981 0.193 2.711 570.2 
Actual emission (conventional) [g/(kW-h)] 9.130 0.282 3.128 763.3 
ESC test emission 
Euro 4 [g/(kW-h)] 15 0.46 3.5 — 
Euro 5 [g/(kW-h)] 1.5 0.46 2.0 z 
ETC test emission 
Euro 4 [g/(kW-h)] 4.0 0.55 2.0 — 
Euro 5 [g/(kW-h)] 4.0 0.55 2.0 — 
Emission index [-] Hybrid bus 
ESC included 2.645 0.419 1.350 - 
ETC included 0.993 0.351 1.350 = 
Emission index [-] Conventional bus (Diesel) 
ESC included 6.086 0.613 0.894 = 
ETC included 2.282 0.513 1.564 = 


The analysis of the emission indexes obtained for the two different buses (Fig. 11) shows that 
the hybrid bus has advantages in terms of both — fuel economy such as lower CO2 emission and 
ecology such as lower emission indexes for CO and HC (lower than 1). This means that during its 
operation the bus releases fewer toxic compounds to the atmosphere than in the dynamic ETC test 
(in reference to a unit of energy). The emission index obtained for NOx amounts to 1.35 — this 
means that the emission of this compound during the road tests is 35% higher than in the ETC test. 
The emission indexes obtained for the conventional bus show that the emission of CO has been 
exceeded (more than 6 times) in the road test. The emission of the other toxic compounds does not 
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exceed the boundary values set out in the Euro 4 standard. It should be noted that for the conven- 
tional bus the emission indexes were obtained through the data related to the Euro 4 standard and 
the ESC test. 


E Hybrid Li Diesel 


Emission index [-] 


co HC NOx 


Fig. 11. Comparison of the emission index of the two buses in terms of real road emissions as opposed to the limits set 
out in the emission standards 


6. Conclusions 


A hybrid bus is more environment friendly because of a lower emission of CO» (fuel consump- 
tion) by 25% and a lower emission of other toxic compounds such as CO by 60%, HC by 30%, 
NOx by 15% as opposed to a conventional bus. 

The analysis of the emission indexes shows that the emission values in the homologation test 
(ESC or ETC for a Euro 4 or Euro 5 vehicle) and the values in the real operation vary. For a con- 
ventional bus the differences of certain toxic compounds (as opposed to the emission set out in 
Euro 4 in the ESC test) are high and amount to: CO 6 times higher, HC 1.5 times lower and NOx 
emission similar. For the hybrid bus (opposed to the emission set out in Euro 5 in the ETC test) the 
following emission values were obtained under real traffic conditions: NOx emission 35% higher, 
HC emission 65% lower and CO emission similar. The elevated NO, emission during the tests can 
result from an inappropriate selection of the engine for this particular vehicle (an engine of higher 
capacity and power would be required so that it could more frequently operate on part loads). 

The defined emission indexes could come in handy in the classification of vehicle fleets in 
terms of toxic emissions according to their year of manufacture, emission limit compliance, mi- 
leage or conditions of operation. 
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Abstract 


The cost of ship operation is dominated by the cost of fuel. The efficiency of ship propulsion depends on number of 
parameters. After selection of the most important we can further analyze them for taking relevant action to increase 
the ship propulsion efficiency. Biofouling of underwater part of ship’s hull has direct influence on vessel resistance 
and engines power required to maintain the speed or keep optimal fuel consumption. HISMAR is a name of the 
project directed on development of robotic system that can be used to introduce new approach to keep ship hull 
clean and almost free of bio-fouling. Frequent cleaning using HISMAR robotic system could be a new tool to 
increase ship performance as well as conduct specific types of ship hull surveys underwater. New approach for 
inspection and cleaning needs field research to validate the financial and technical benefits. 

The HISMAR Project no 012585 was founded by EU-FP6 Framework Program. 


Keywords : underwater cleaning, robotic system, underwater inspection, ship propulsion, efficiency 


1. Introduction 


We observe continuous growth both the global and EU shipping fleet. When speaking about 
the measures that can be taken to improve the shipping efficiency one of major issues is hull 
performance increase by the reduction of the wave and frictional resistance and increase of 
propulsion system efficiency. Whilst shipping is one of the cleanest forms of transportation, the 
fouling of ships reduces efficiency, increases fuel consumption and in extreme cases - corrosion. 
To minimize wave drag, special attention is given the smooth ship hull design and construction 
quality. The numerous efforts focus on frictional drag reduction using not only computational 
techniques. Modern anti-fouling paints and special low friction coatings are another large area 
of investment into research and achieved spectacular performance. When speaking about the 
increasing efficiency of the propulsor system there is a number of another opportunities like 
new highly efficient diesel engines, advanced propeller designs as well as promoted by some 
companies propeller coating techniques that prevent biofouling. Some guidelines and draft data 
regardding potential savings due to lowering of fuel consumption can be found in technical 
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literature. An example taken from data collected by Wartsila is presented in Table 1. The 
important comment to that table is that, there is very limited amount of data available that could 
be referenced as being confirmed by measurements at real sea conditions. Bearing that in mind, 
when using such data for calculation of potential savings we have to treat that figures as 
guidance. More detailed data can be collected during research on ship performance using real 
object data gathered during ship operation. This is quite complex task and will need certain 
financial resources and positive approach from ship operators to collect and validate the data 
with assistance of ship crew. 


Table 1. Estimated savings due to selected measures applied to various types of ships 
Ship Type | Tanker | Container Ro-Ro Ferry 
Efficiency measures 
Hull clearing < 3% <2% <2% <2% 
Propeller clearing and polishing <10% <10% <10% <10% 
Modern hull coatings < 9% < 9% <5% <3% 
Propeller efficiency measurement < 2% <2% <2% <2% 
Constant versus variable speed operation < 5% <5% <5% <5% 


Source: Wartsila [1] 


It is obvious, that machinery that operates for longer time and at higher power and loading is 
more susceptible to failures. Lower engine loading in general means lower wearing of certain 
components of ship drive system that could lead to longer operational life of the ship propulsion. 
Furthermore, it results in better reliability of the main propulsion and resulting ship safety 
improvement - particularly important for tankers, bulkcarriers, ro-ro and passenger ships. 
Monitoring of the ship hull condition for estimation of fouling or coating roughness 
measurements when coupled with means for cost effective ship hull cleaning (an example could 
be HISMAR system), could be valuable operational tool that could be used for ships in 
operation. Such procedure could be a part of more complex propulsion efficiency audits being 
a part of modern ship hull survey procedures. 


2. Environmental and legal aspects 


There is ongoing discussion at International Marine Organisation (IMO) with regard to Marpol 
Annex V (Regulations for the Prevention of Pollution by Garbage from Ships). IMO has already 
defined the residue from hull scrubbing as ships waste and it is therefore covered by the 
convention. As a consequence of this, it is likely that restrictions will be introduced as to where 
and how ship hulls may be scrubbed. The scrubbing of hulls by divers is forbidden within many 
European harbours. Fouling on a hull increases a ship’s drag through the water, thus increasing 
the amount of power required to maintain the same speed or reducing the speed of the vessel for 
a given power. Fundamental data on that subject can be found in reference [2]. It is estimated 
that serious fouling can increase the drag on a ship by up to 40%, reducing the speed by up to 2 
knots and increasing the fuel consumption by 10+20% [3]. For the shipping industry, this 
means increased costs and time delays. The extra fuel consumption also increases greenhouse 
gas emissions, being composed of NOx, SOx and CO, from ships by up to 20 million tonnes 
per annum. Although shipping produces the least amount of greenhouse gases annually when 
compared to other modes of transportation, under the Kyoto agreement the EU is committed to a 
5% reduction in emissions by 2012. 
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3. Underwater ship hull cleaning -- advantages and tools 


For over 130 years people have thought about machines or robots for cleaning ship hulls. 
Various devices and technology has been developed and tested. Only few developments have 
found wider application, an example could be hydraulic brush system operated by diver 
(BrushCart). In recent years, we are observing growing interest in robotic system that are 
developed to cope with the bigger ships and eliminate diver’s presence underwater. Important 
issue is stress on marine ecology requirements. Robots are able to cope with that case using 
suitable tools. However, the investment in robotic system may be substantially high, such 
system will be able to work continuously underwater and possibly above the water where divers 
do not have access. Plug in sensor modules could allow to conduct detailed inspection of the 
almost whole hull during the same mission. 


3.1 Brush technology 


Brushes are used in cleaning carts, handheld polishers and some robotic systems and are 
able to cope with almost all types of fouling. Most systems consist of one or more rotating 
brushes pneumatically, or hydraulically driven. This requires the minimum of equipment 
beyond the cleaning device itself thus reducing the cost of the system. Before the ban on 
tributyltin (TBT) came in, brush technology was preferred to underwater jetting systems, as it 
was easier and more economical to use. However, the increase in use of environmentally 
friendly low friction coatings can cause a problem, as these coatings are less durable and more 
easily damaged by the abrasive action of the brushes. Research has shown that bristle density, 
angle and gauge have a greater effect on shear and normal forces produced by brushes, while the 
brush speed and stand-off distance has little or no effect. The main point demonstrated by the 
research was the selection of the brush cleaning system and forces involved is dependant on a 
number of factors and their relationship is very complex. The major problems in cleaning using 
heavy duty brushes could have place when dealing with calcerous forms of fouling. 


3.2 Water jetting technology 


The use of high pressure water jets has become an accepted alternative to brush cleaning 
systems. Unlike the brush-based systems, water jets can be easily controlled by reducing or 
increasing the pressure from the pump. A water jet's effectiveness is dependant on the surface, 
pressure of water, jetting angle and distance from the cleaning surface. Jet nozzles, such as 
CaviJet or SwirlJet have been developed to enable effective cleaning of the hull underwater. 
Tests using cavitating water jet nozzles showed that the cleaning process can remove various 
types of fouling from hull coatings, while at the same time, minimizing the damage to the 
coating. Although jet washing provides increased control of the cleaning process, the perceived 
increase in the cost of the equipment is still thought to be prohibitive. In hase of the HISMAR 
system low pressure jetting will be sufficient to remove effectively and safely the layer of slime 
from the hull surface. 


4. Cleaning problems 


Diver(s) presence enhances the risk of diving accidents at work. Some data revealed by HSE 
show that probability of diver accidents during professional diving operation is higher than in 
farming or civil engineering. The number of accidents in offshore and inshore diving is in a 
range 20-40 for 100000 dives while fatal accidents occur in number of 6-7 for 100000 dives. 
Additionally to the risky job, diver has limited time to be able to work underwater so a team of 
divers is required to perform cleaning task. An open issue is the cost of diver safety measures 
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and environment protection requirements that must be taken into consideration discussing the 
cost benefits. Remotely controlled machines are not able to clean the whole wetted ship surface 
due to some hull features that restrict the robot operation. Ship bow and stern, due to hull shape 
are very difficult for automatic cleaning and that area will probably be cleaned by divers. 
Generally it is assumed that cleaning of 80% of the ship underwater area from slimes could 
provide suitable effect of drag reduction. 


5. HISMAR idea 


The abbreviation HISMAR stands for Hull Identification System for Maritime Autonomous 
Robotics. HISMAR is intended to be a multifunctional robotic platform which will be able to 
perform specific inspection or maintenance tasks such as structural integrity monitoring of the 
ship’s hull or cleaning operations. Apart from Project leader - Newcastle University, the other 
partners in the project were Graal Tech of Genoa; the UK's Shipbuilders and Shiprepairers 
Association; TecnoVeritas of Portugal; Polski Rejestr Statków; Robosoft of France; Carnival; 
Moscow State Technological University; Royal Thai Navy and TEPAC Technology & Patent- 
Consulting of Germany. 

Robotic platform is to be deployed for the board of the ship, harbour service craft or from 
the pier using simple care or special launching and retrieval device. The control over the vehicle 
is provided via the special umbilical with power, control lines and hoses used for removal of 
cleaning wastes to the surface. Intervention with the use of HISMAR robot needs some 
preparation before the vehicle is placed on the surface of the ship hull. The desired situation is 
the case when digital data of the ship hull construction are provided to the robot control 
computer. On the other hand, at start of the job a map of the hull is automatically charted, 
recording the location of every weld, thickness change, rivet and indentation on the ship’s 
surface. Adjustable jets of pressurised sea water blast the marine growth off the surface of the 
ship which is then sucked up into the main chamber. Here, ca 150 litres of water a minute is 
filtered and the bio-fouling removed and rendered harmless to the local environment. In this 
way, the ship’s robotic ‘vacuum’ can continuously roam the ship’s hull, preventing the build up 
the layer of slime. 

Hull surveying is an important part of any vessel's life span and a number of periodic 
inspections of the hull are required during the vessel's life. Currently, the minimum requirement 
is for a visual inspection of the hull and with some thickness measurements being taken in 
specific areas of the hull, or where a probable defect might have occurred. These are usually 
performed using an ultrasonic sensory system that is placed on the plate's surface. Due to the 
size of the detection head and the skill required to operate the equipment, only a small 
proportion of a vessel's hull can be accurately measured. A full hull inspection is required to be 
performed in dry dock every five years, but up to 20% of the hull may not be inspected due to 
the vessel's dock supports. Between these class renewal inspections, intermediate hull 
inspections are required. These are general visual underwater inspections performed by divers. 
However, in recent years a number of robots have been developed to improve the accuracy, 
coverage and reliability, while reducing the time and cost of the inspection. The current robotic 
systems available are limited to visual inspection and ultrasonic plate thickness measurements. 


6. The critical components of the HISMAR Robotic Platform 


The HISMAR robot design is presented in Fig 1 and 2. The key systems of the robot 
include: 
Drive Systems & Central Robotic Platform 


Key to the HISMAR robotic platform is the versatile central drive module. This incorporates the 
robot’s drive systems, navigational sensor systems and control electronics. Robot is linked with 
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the surface by a mean of specialised umbilical that contain power and data conductors and 
possibly hydraulic hoses for removal of cleaning debris. 


Cleaning & Debris Extraction Systems 


The cleaning system will utilize a water jetting to clean the hull. The jet spray system is 
intended to provide the customer with sufficient control for their cleaning needs while 
preventing damage to the hull coatings. A complete extraction and waste handling system is to 
be developed to comply with current and future marine environmental legislation. 


Mapping & Navigation Systems 


HISMAR’s unique mapping and navigation system will allow for full autonomy of the robot 
after the initial process of mapping the hull. This will allow for complete or partial cleaning of a 
ship’s hull, the latter being complete at the crew’s convenience with no loss of operating time. 
The map will be a permanent record of the condition of the ship’s hull, being updated by the 
navigational system whenever a cleaning operation is performed. 


7. HISMAR robotic system advantages 


Marine growth on ships is a huge environmental and financial problem for the marine 
industry and HISMAR offers a unique solution to both of these — and more. Created is a 
system that works totally independently — in or out of the water — and not only keeps the ship 
clean but also feeds back vital information about the hull’s condition. Because the map it 
follows is so detailed, if there is a change to its path caused by corrosion or a crack in the steel 
then it feeds this information back. This means it can be used as an additional check the 
condition of ship’s hull and provide important reference data for classification society 
surveyors. All other developed cleaning or inspection systems currently available are remotely 
controlled during their operation, requiring highly skilled and experienced operators to 
effectively clean the hull, while the ship is out of operation and usually in dry dock. The 
advantage of the HISMAR robot is that it is an autonomous system so it can continue cleaning 
with the ship remaining in service. 

The platform can be launched whenever the vessel is in port or at anchor. The device will be 
able to complete its tasks partially whilst in one port and be re-launched at successive points to 
complete the task. The generic platform will offer the option of using targeted plug-in modules 
to perform specific inspection or maintenance tasks such as structural integrity monitoring of the 
ships hull or carrying out cleaning and waste recovery operations. This project offers a means to 
effectively and efficiently undertake hull inspection and maintenance thereby extending the safe 
working life of the vessel. Cleaning of the hull ensures the vessel maintains the lowest possible 
hydrodynamic resistance and consequently reduces amount of fuel oil consumed. Therefore 
ensuring a clean and smooth vessel underwater hull surface reduces vessel emissions and 
reduces operating costs. 


8. Validation of fuel consumption data and effectiveness of cleaning 


There is noticeable amount of published data providing general information about cleaning 
benefits but only few data available that provide reliable figures of possible savings due to 
periodical cleaning of the ship hulls. A very optimistic diagram that presents possible 
advantages from cleaning ship hull is given in Fig. 3. Some companies are offering advisory 
optimization services of vessel performance and in some cases also hull cleaning underwater 
using hydraulic tools. Up to our knowledge, only one company is using specialized ROV system 
and no reliable commercial or scientific data regarding its performance are available. 
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Taking decision about cleaning requires certain data to be analyzed. US Navy underwater 
husbandry practice requires the decision to be taken after underwater inspection and precise 
estimation of the area covered by fouling after comparison with reference data. Other 
commercial practice requires analysis of periodically collected voyage data as well as basic data 
about the prime mover operation like — rpm, specific heating value of fuel oil, fuel oil 
consumption, power measured by special device or torsiometer, turbocharger revolutions, 
exhaust gas temperature — just as example. In that case, the collected data are subject of 
evaluation using dedicated software by commercial company on a subscription basis. In case of 
some ships and sea routs potential savings are big enough to pay for kind of ship services.. 


9. Conclusions 


The use of remotely controlled technology for underwater cleaning is subject of numerous 
R&D activities. Commercially available devices are not cheap and the cleaning in dry dock is 
much more popular and usually coupled with class surveys. Prime target for robotic tool is to 
achieve performance and cost of operation at the competitive level. Hull cleaning by cost- 
effective robotic devices that are able to cope with slime will limit fouling and lower the fuel 
consumption. The use of intelligent crawling underwater robots will expand as there will be 
more comparable data available describing ship propulsion performance increase as a result of 
more frequent robotic cleaning. Collection of reliable data about the financial and technical 
benefits of the robotic cleaning technology will be important for independent assessment of 
ship propulsion system performance. Plate thickness measurements and other data describing 
ship hull condition that will be collected by robot during cleaning missions could be further 
used by classification societies for class renewal surveys. The hull surveying using innovative 
robotic technology is in development phase and needs additional founds to be recognized as 
mature technique recommended for wider use. 
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Fig, 1 General view of HISMAR robot prototype 


Fig, 2 General view of HISMAR robotic platform 
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Abstract 

Marine transport development and related to it, an increase in emission of harmful compounds forming when 
burning marine fuels, enforces conducting research on that emission in sea regions located in the vicinity of 
agglomeration areas. The basis for formulating the models of emission of harmful compounds in marine engines 
exhausts is formulation and study of the vessels movement models. The models ought to take into consideration the 
specificity of the processes taking place while a vessel is moving, which influence the formation of harmful compounds 
and their emission intensity. Emission models can be later used in formulating models of pollutants dispersion and 
their immission in the areas distant from the source of pollution, e.g. in city agglomeration regions. The paper 
indicates the possibility of using the Automatic Identification System as a source of statistical data, which can be used 
to define parameters and characteristics of the vessels movement. Identification of such a model was illustrated by the 
examples of the Gdansk Bay region. 


Key words: movement models, emission, toxic compounds in exhausts, marine engines 


1. Introduction 


Emission of harmful compounds in engines of different applications constitutes a significant 
ecological problem of the 21* century. For a few decades there have been made many attempts to 
describe emission and its intensity with some specialised mathematical models [1-5]. When 
mobile sources of emissions are considered, majority of such models are formulated to conduct 
research on emission of motorization origin [1-3], both for road vehicles and off-road vehicles 
(mainly construction equipment and agricultural machinery). There exist also models of emission 
for aircraft vehicles (planes, helicopters). However, relatively few attempts are made to perform 
spatial modelling of pollutants emission distribution for the area where the source of these 
pollutants is the vessels movement. The problem of emission of pollutants in marine engines 
exhausts is all the more essential, as contrary to observed global land-based emissions decrease, it 
is projected that international shipping emissions will continuously be growing [6]. Such a 
situation results from dynamic development of marine transport — it is assumed that over 90 % of 
the world’s total trade is carried by sea [7]. 

Movement models to estimate emission of pollutants from engines of motorization and aircraft 
origins, which have been formulated for many years, constitute a significant input into the 
development of modelling the pollutants emission from combustion engines. However, due to 
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significant differences in topographic, hydrometeorological conditions and the specificity of 
vessels maintenance, they cannot be directly applied in estimation of emission of these pollutants 
in marine engines exhausts. 

Existing inventories of distribution of a mean number of vessels in a given sea area in given 
period of time, although they give a general idea on the vessel traffic intensity in a given area, they 
are often too general and imprecise, especially in the cases when the research do not aim at 
determining the mean number of vessels in a given sea region, but aim at describing the trajectory 
of individual vessels operating in a given sea area at given time interval. Modelling should also 
take into consideration technical parameters of an individual vessel (type, load, power etc.), which 
would be helpful in modelling the emission and dispersion of toxic compounds in the vessel 
exhausts. 


2. Deterministic model of the marine vessels movement 


The length of the segment joining any two points P,P, i=1,...,N , of coordinates (x,y;), 


between which a vessel moves: 
[PoP E) 20, )F HDE -IE > (1) 
and after transformations (applying the Lagrange’s theorem about average value): 


|P P, l= yix (0, JF +b) At, (2) 


where 6, , % are points from interval [7,_,,¢,), and At; =t; —t, ,. 


Analysis of emission and its intensity of the compounds in exhausts, should begin with 
determination of the randomly changing vessel speed at a given time /, since the emission intensity 
depends on this value. Vector stochastic process {V(t)=(V,(t),V,()): teT; describes the 
randomly changing vessel speed. At a given time ¢, this value is a two-dimensional random 
variable indicating a momentary speed vector. It should be noticed that between the random 
processes of vessel trajectory: +S(t) = (X(1),Y(1)): t20} and speed {V (£) = Vx (A), Vy): t20 
there exist the obvious relations: 


BOLYA, t20, S(t) = f V(x)dx, 120 G3) 
dt 0 
From equation (3) there come the following formulas: 
dX(t dY(t 
v,(1) = eh V,(1) = — (4) 
Realisation of the process V(t) can be thus described by the following formulas: 
pgn dy(t) 
7 dt > v(t = e (5) 
After transformations, the length of the speed vector is determined: 
Iv Flv OF +, OF (6) 


In modelling the marine vessels movement, there can be assumed that in time intervals 
[t,,,4,), i=1...,N the vessel speed is constant: | v(t)|=v,,i=1,...,N, and therefore the 
intensity of their exhausts emission will be constant: f(x(),v)=4,, teltat) i=l,..N. 

In reality, due to the yawing, the vessel route from point P,, to point P along a given 


trajectory, differs from the trajectory determined theoretically by drawing a line from point P_, to 


point 7;. The yawing can be written as: 
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where: y,,, — yawing amplitude; V — vessel speed; T — period of dominant harmonic of 
yawing. 
Therefore, due to yawing, the vessel route from point P, to point P is longer than 


As, =v,At, by a random value D,, i=1,...,N , and therefore the mass of exhausts emitted at this 
route interval will increase by a random value Am, =4,D,, i=l..,N. There can be accepted a 
preliminary hypothesis that the random variables D, i=l,..,N are mutually independent and 
have Gaussian distributions N(d,o,), i=1,....N, where d;=éAs;, 0;=pAs;, i=1,..,N. 
Factors €, are the coefficients of the average route elongation and of the standard error of 
route elongation. 


3. Stochastic model of vessels movement 


Number of vessels operating in a given sea area varies in time, which constitutes therefore a 
stochastic process. Such a process is determined with a symbol {X(t):t¢7}. The value of this 
process at time instant ¢ means the number of vessels operating in a given sea area at instant t. 
Realisations of this process are the functions of constant intervals. Any change of states appears at 
a moment of entering or leaving an analysed sea area by a vessel. The first approximated model of 
this process can be assumed the non-homogenous Markov process of discrete set of states 
S = {0,1,2,...}. 

Stochastic process {X(t):t eT) of finite or denumerable set of states is called the Markov 
process, if for any £,,f,,...,£,,£f,,, ET suchas ft) <4, <...<f, <¢,,, and any i, 7,i),/,,...,i,., ES, 

PUX (ta) =F | X(t) =i, X Ga) Shy X (to) =i} = PAX Gad) =J XG) =H (7) 


obot, are interpreted as the past instants, £, as the present instant, and £,,, as the 


orf ny 


If t 


incoming instants (future), then the equation defining the Markov process means that conditional 
distribution of ‘future’ states does not depend on the ‘past’ if the present state of the process is 
known. 

If T =N, = {0,1,2,...}, then the Markov process is called the Markov chain, and if 
T = R, =[0,0), then it is the Markov process with ‘continuous time’. 
Conditional probabilities 


Pitar = PX (yas) = JX) =I EJ ES (8) 
are called the probabilities of passing from state i atinstant 7, tostate j atinstant ¢,,,. Let 
t, =S, t, =t, where 0 <s<ż. From the definition of Markov process there result the following 
corollaries: 


1) P,(s,t)20 forall i,jeSands,teT, s<t. 
2) $ P,(s,t) =1 for all ie Sand s,reT, S<t. 


jes 
3) P,(s,t) = XP (s,u)P,(u,t) for all ijeS and s,u,teT, s<u<t, (Chapman- 
keS 
Kolmogorov equation). 
Markov process is called homogenous, if probabilities of passing depend only on the 
difference of instants (t-s): P, (s,t) = P,(t-s).For s=0, P,(0,t)= F,(t). 
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Let {x (t):¢ e[0, oo) } be Markov process of discrete set of states S. Functional matrix, whose 
elements are the probabilities of transition: 


P,(s,t)= P{X(1) = j|X(s)=i},i,7€S, OSs <1 (9) 
will be designated as: 
P(s,t) =[P,(s,1):i,j € S]. (10) 
The Chapman-Kolmogorov equation in matrix notation will be as follows: 
P(s,t) = P(s,u)P(u,t) (11) 


It can be assumed that there exist uniform (relative limits 


1-P,(t,t+h 
4,(0) = lim ult ) ieS, (12) 
h>0 h 
_ P,(t+h) 
%0) = lim jes (13) 


and they are continuous functions of parameter t. 
For homogenous Markov process, the intensities of transitions are constant A, (t) = A, . 


The Chapman-Kolmogorov equation can be written in the form of: 
P (s,t+h)= X} P,(s,)P (t,t+h), where ijeS and 0<s<t<t+h. (14) 
keS 
Subtracting from both sides of the equation P,(s,£) and dividing by A, there is obtained: 
P,(s,t + h) — P,(s,t P,(t,t+h P.(t,t+h)-1 
AG ? AG |= TB (0, p fg OO (15) 


kzj 


Passing to the limit at A > 0, assuming that in case of Markov process of countable set of states, 
the convergence defining the intensities 4 (t), k,jeS at every determined j is uniform in 


relation to k, the following system of differential equations is obtained: 


OP..(s,t 
EN Se i: ijes (16) 
ot keS i 
The initial condition constitutes the natural equality: 
1 dla j=i 
P.(s,s)=6, = 17 
pee k da gai ua 


In a similar way, the can be derived some differential equations allowing to calculate one- 
dimensional distribution of the process: 
PO=PixM=j, jes. 
From the formula for total probability, it is obtained: 
P(t+h=D POP (Gt+h), jes (18) 
ieS 
Subtracting from both sides P, (ż) , dividing by h and passing to the limit at h — 0, the following 


system of equations is obtained: 

P= POLO, jes (19) 

ieS 
The initial condition in this case constitutes the initial distribution of the process 
P(0)=p), ieS. 
For homogenous Markov process, Kolmogorov system of equations (18) is brought to the system 
of linear differential equations of constant coefficients: 
PD=XP(Qa,, jes (20) 


ieS 
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Such a system of differential equations is the most often solved using operators, applying 
Laplace’s transformation: 


UFO: 5G)=|POe"d. (21) 


Making use of the property: L[P/(1)]= s P (s) — P,(0), the Laplace system of differential equations 
undergoes transformation. As a result, there is obtained the following system of linear equations, 
where unknown are the Laplace transforms: P(s), ieS: sP (s)- p? = Y P(s)A,. jes 
ieS 
[8]. 
One of limiting theorems for homogenous Markov processes with ‘continuous time’ at time 


t—> œ goes as follows: 
Let {x (t):t e[0, oo) } be a homogenous Markov process of finite set of states S and of the intensity 


matrix A=[A, :i, je S]. If there is such a real positive number t,, that at least one column of the 
matrix W=[p,(t,):i,j €S] contains one column consisting of all the positive elements, then 


there exist limiting probabilities: 


lim P,(t)=limP,()=P,, jes (22) 
and they meet the system of linear equations 

> BA, =0, jes, ÈP =1. (23) 

ieS jes 


Therefore, the limiting distribution is determined by solving the system of linear equations, whose 
coefficients depend only on constants of the transition intensity. 

Using non-homogenous Markov process as a model of the movement of objects navigating in 
a given sea area is possible it there is assumed that the Markov process {X(t):t¢ 7}, the value of 
which at instant ¢ constitutes the number of vessels operating in a given sea area at instant ¢, isa 
single process, which means that in a short time interval the state can change by one. The matrix of 
the intensity of transitions [8]: 


A(t) =[2,(0) :i, j € S =40,1,2,...)] (24) 
of this Markov process is as follows [8]: i 
-a (t) a(t) 0 0 0 0 
at) —(a,(1)+ 8,(1)) B,(t) 0 0 0 
0 a, (t) -(2,(1) + 8, (1)) B,(t) 0 0 (25) 
0 


A(0)=| 0 0 a(t) —(@,0+B,0) BO -- 


form is practically impossible. However, it is possible to find an approximated solution, assuming 
that the elements of the matrix of the intensity of passages are constant in given time intervals and 
that the set of states of the process is finite S = {0,1,2,...,7}. Asymptotic distribution is received 


when solving the system of equations (23), in which given coefficients A, are the elements of the 


matrix of the intensity of transition: 
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-a Ay 0 

a —(a@ +f) B, 

0 a, — (0, + B,) B, 0 
A(t)=| 0 0 a, -(0,+B,) p 


(26) 


o© © © © 


A, Ana +r) Ba 
B, =, 


The solution of this system of equations, having replaced the constants with the functional 
coefficients is of the following form [8]: 


1 
P(t) 


“ay (£) a, (1)... a(t) 
a AOL (1)... AO (27) 
_ A (tat)... & iO p 
LOLO- BO | 
The number n is the maximal number of vessels in the sea area. 
If in a given time interval [æ, 2], the distribution of the number of vessels in the sea area is 


P(t) =1,2,...,7 


constant, then the number K can constitute the integer part of the expected value: 
E[X(0]=Żk Py (28) 
k=0 


In order to realize the described above stochastic model of movement of the vessels operating 
in Gdansk Bay region, it is necessary to know the number of vessels operating in the analysed 
area, their distributions in relation to the vessels types, their size, speed, power and kind of main 
engines etc. 


3. Estimation of some parameters of a model of emission of harmful compounds in marine 
vessels exhausts 


Estimation of the intensity of movement of the vessels operating in a given sea area at 
particular time interval, and on its basis — determination of emission of harmful compounds in 
exhausts of the vessels engines, is possible with the use of data transmitted in the Automatic 
Identification System (AIS), registered by the shore installation, which traces the vessels 
movement in a given sea area. For the marine vessels passing through certain test segments (fig. 
1), situated perpendicularly to shipping lanes approaching the ports of Gdynia and Gdansk, as well 
as the shipping lane splitting into those two shipping lanes (in the proximity of Hel port), there 
were analysed the movement parameters of the vessels operating in a given sea area at given time 
interval. 
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Fig. 1. Map of Gdansk Bay area where there are marked the test intervals [9] of geographic coordinates: 
AB = (54,58325°N, 18,85362°E) (54,45917°N; 18,88191°E) — gate GD; CD = (54,54200°N, 18,69511°E) 
(54,52508°N; 18,69326°E) — gate GA; EF = (54,48666°N; 18,70878°E) (54,47326°N; 18,73178°E) — gate HE. 


There were examined the vessels passing in the period of one month (from 1* to 30" June 
2006) through the test intervals (called later: ‘gates’) shown in fig.1. In the analysed Gdansk Bay 
region (fig.l), the total number of vessels operating in June 2006 amounted to 627 (fig.2). On 
average, during 24 hours, to the analysed area through the gates there were entering from 3 to 8 
vessels (in case of vessels entering the analysed area through the GD gate and GA gate), and from 
9 to 14 vessels (in case of vessels entering the analysed area through the HE gate) (rys.3). Majority 
of these vessels were entering the analysed area in the afternoon (through gates GD and GA) and 
at time interval 00.00-12.00 (through gate HE) (fig.4). 
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Fig. 2 A 24-hour distribution of the number of vessels: a) outward bound passing through gate GD . b) outward 
bound passing through ‘gate GA’; and c) passing through ‘gate HE’ and going towards of Gdansk port or Gdynia 
port (in the period of 01-30.06.2006) 
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Fig. 3. A 24-hour distribution of the number of vessels: a) outward bound passing through ‘gate GD’; b) outward 
bound passing through ‘gate GA’; and c) passing through ‘gate HE’ and going towards of Gdansk port or Gdynia 
port (in the period of 01-30.06.2006) 
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Fig. 4. Distribution of the number of vessels: a) outward bound passing through gate GD'; b) outward bound 
passing through ‘gate GA’; and c) passing through ‘gate HE’ and going towards of Gdansk port or Gdynia port, 
depending on time intervals during 24 hours (in the period of 01-30.06.2006) 


Analysis of the available statistical data allowed to determine values of momentary speed v* of 
the vessels entering the analysed area in analysed period of time. As it can be noticed, the most 
frequent speed of the vessels entering the area through individual gates, was from 8 to 18 knots 


(fig.5). 
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Fig. 5. Distribution of the momentary speed v* (velocity over ground) [in knots] for vessels a) outward bound passing 
through ‘gate GD’; b) outward bound passing through ‘gate GA’; and c) passing through ‘gate HE’ and going 
towards of Gdansk port or Gdynia port (in the period of 01-30.06.2006) 


The volume of the underwater part of a hull V was calculated for the analysed group of 
vessels on the basis of assumed simplified shape of a hull [10]. On the basis of the obtained 
results, it can be stated that most of the analysed vessels were the ones relatively small — of 
underwater hull volume up to 25000 m” (fig. 6). 
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Fig. 6. Distribution of the underwater hull volume of the vessels a) outward bound passing through ‘gate GD’; b) 
outward bound passing through ‘gate GA’; and c) passing through ‘gate HE’ and going towards of Gdansk port or 
Gdynia port (in the period of 01-30.06.2006) 


From the analysed group of vessels there were chosen the ones, which in the period of 01- 
30.06.2006 were passing through gates GD and GA and going towards gate HE (in total there were 
289 vessels). On the basis of analysed change in momentary speed v” of individual vessels, read 
when they were passing consecutive two gates, there was calculated the probability of the 
occurrence of the speed change. There were taken into account the changes of the range of 0.1 
knots. In figures 7 and 8, it is shown the probability of occurrence of an increase and decrease in 
speed (by x-knots) on the route between individual gates. Function y, described by a 4" order 
polynomial, means that probability, and R” means the adjustment coefficient. 
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Fig. 7. Probability of occurrence of an increase in momentary speed (by x-knots) for the analysed group of vessels 
operating at the route between gates GD-HE and GA-HE (in the period of 01-30.06.2006) 
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Fig. 7. Probability of occurrence of a decrease in momentary speed (by x-knots) for the analysed group of vessels 
operating at the route between gates GD-HE and GA-HE (in the period of 01-30.06.2006) 


From the analysis of the charts shown in fig. 7 and fig. 8, it results that the changes in the vessels 
speed on the route from gate GD to gate HE and from gate GA to gate HE are slight, and the 
probability of occurrence of the changes in speed by more than one knot is very low (below 0.12). 


4. Conclusion 

In the paper, there were presented formulation and estimation of parameters of models of 
movement of the vessels operating in the Gdansk Bay region in a particular time period. The 
analysis was possible with the use of statistical data obtained via Automatic Identification System 
(AIS), which appeared to be a tool that can be effectively used in such purposes. The proposed 
method is possible to be applied in the process of determination of an approximated level of the 
toxic compounds emission in marine engines exhausts in given se region at particular time. 
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Abstract 


The area of the Gulf of Gdansk, just like port towns and coastal regions, is vulnerable to pollution contained in 
the exhaust gases of industrial plants, power plants, vehicles and harmful compounds contained in exhausts of vessels. 
In order to determine the share of vessels in polluting atmospheric air and to counteract the harmful effects of toxic 
compounds in marine engine exhausts it is necessary to know the values of those compounds’ emission from particular 
vessels, which is possible by knowing the movement parameters of vessels, the values of concentrations of particular 
compounds for these parameters and the atmospheric conditions in the region of their occurrence. 

The object of balancing emission of pollution contained in the exhausts from engines propelling vessels are the 
processes of global emission, averaged in a sufficiently long time period, which is determined first of all by the effi- 
ciency of averaging the variable conditions of the objects’ operation. Due to the impossibility of making emission 
measurements on vessels in the region of research, it is necessary to work out suitable models of estimating emission 
of toxic compounds, based on proper measurement tests simulating navigational conditions in the Gulf of Gdansk as 
well as mathematical models. 

The work presents conditions related to modelling the emission of exhaust compounds from vessels engines in the 
region of the Gulf of Gdansk, with particular consideration to estimating the values of toxic compounds concentration. 


Keywords: emission, dispersion, modelling, marine engines, exhausts 


1. Introduction 


The problem of air pollution in ports and approaches to ports is important inasmuch as the 
ports are usually close to or on the territories of large cities, and their restricted area causes large 
concentration of vessels in a small area. Widely conceived operational conditions are not insignifi- 
cant either, among them being the way of engine operation, the frequency of occurrence and char- 
acter of non-stationary states, transitory processes characterised by considerably larger emission of 
toxic compounds than during sailing in open areas with constant engine load. The way of opera- 
tion is also affected by external conditions, that is the effect of marine environment on the engine’s 
work. Certainly the kinds of fuel and lubricating oil affect exhaust toxicity in no smaller degree. 


The research currently carried out concerning atmosphere pollution caused by the emission of 
harmful compounds from traction engines (airplane and car engines) [1-3] provides very large con- 
tribution to the development of modelling of imission of harmful compounds emitted from diesel 
engines, but due to the difference in both topographic and hydrometeorologic conditions and the 
specificity of vessel operation, they cannot be applied for estimating imission in coastal areas [4]. 
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Modelling emission of particular exhaust components from exhaust systems of marine en- 
gines during movement in a particular area constitutes one of the key problems of balancing and 
making input data for dispersion models. 

The object of balancing emission of pollution contained in the exhausts from engines pro- 
pelling vessels are the processes of global emission, averaged in a sufficiently long time period, 
which is determined first of all by the efficiency of averaging the variable conditions of the ob- 
jects’ operation. 

The factors determining global emission of substances contained in marine engines ex- 
hausts have been classified and described in [4,5]. 

The process of modelling emission and dispersion of toxic compounds (ZT) in marine en- 
gine exhausts (which are the elements of modelling the imission of these compounds), is very 
complex and requires the knowledge of five basic parameter groups [5,6]: 

e vessel parameters: length, breadth, draft, technical condition of propulsion system, kind of 
propulsion (including kind and number of engines), kind and number of propellers, etc.; 

e vessel movement parameters: speed and course; 

e external conditions: wind direction and force, air and water temperature, atmospheric pres- 
sure, air humidity, state of the sea; 

e number of vessels, with consideration to category; 

e topographic parameters of the area: terrain kind and profile (water, land). 

Models of emission from land means of transport created in Europe, like HBEFA, 

COPERT, DVG and DRIVE-MODEM [7,8] endeavour first of all to take into account the 

largest possible number of parameters affecting emission, yet with such a large number of 

factors and the complex description of phenomena determining the emission process, sim- 
plifying assumptions cannot altogether be avoided. 


2. Modelling emission of exhausts from marine engines 

In work [9] there has been considered the trajectory of vessel movement as the realisation of a 
two-dimensional stochastic process {S(f) = (X(1),Y(t)): 120}, assuming that the process is of mul- 
tidimensional distribution of continuous type and continuous realisations. The realisation of such a 
process is a two-dimensional trajectory dependent on time {s(4) =(x(4),v(@)): teT}. 
The equation describing the mass of exhausts emitted can be presented as follows: 


B 
M=[f(s(0)| v(0)| dt, (1) 


where: 
v(t) length of vessel velocity vector. 


Assuming that in the time intervals [f,,,7,], i=1L...,N the speed of the vessel is constant 
Ivey, i=l N , (2) 
and the intensity of exhaust emission is constant 
SOIE = Vs tef[t,,,t,), i=1..N, (3) 


the lower estimation of emitted exhaust mass is obtained 
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M=YM., , (4) 


where: 
M,=y,v,At,, i=L12,.,N. 


The mass of emitted exhausts in particular area A in time interval [£;.,, t;] is the sum of masses 
emitted by all vessels in this time interval in the area. If wo, k =1,...,K denotes the mass of 


exhausts emitted by the k-th vessel, then the total mass of exhausts emitted in area A in time inter- 
val [£;.,, ti] is a random variable 


K 
Wx = > (5) 


k=1 


Random variable W, as the sum of independent random variables with normal distribution 
has normal distribution with expected value 


K K K N 
E(W,) = LW => E(AM®O)+M®]= XM. + cP 7 As] 
k=l k=l k=l i=l ; (6) 


Standard deviation of these random variables is equal to 


o(W,) = 2 PAM ©] = JE AP, (7) 


i=l 


In this model the number of vessels K has been assumed to be constant. In reality the number 


of vessels varies randomly in time, so it is a stochastic process [10]. 

Work [11,12] presents a method of approximate identification of vessel movement model, and 
for estimating the parameters and characteristics of exhaust emission from vessels sailing in the 
Gulf of Gdansk region it is necessary to know the number of vessels navigating in the region ana- 
lysed, their distribution with regard to kind, size, speed, power and kind of main propulsion en- 
gines etc. 


3. Modelling load of vessel main propulsion engines in order to estimate the value of toxic 
compounds emission in the exhausts 


The amount of emitted harmful compounds in marine diesel engines depends on values de- 
scribing the state of engine work like torque Mo, engine speed n, thermal state of the engine J, 
technical state of the engine Z (parameters of load exchange system, state of TPC system, techni- 
cal state and proper regulation of injection apparatus), condition of environment G (e.g. ambient 
temperature, pressure, air humidity) and variable sailing resistance of vessel O (vessel resistance in 
shallow water, vessel resistance while moving in a canal, air resistance and waving effect). It can 
thus be written down that the emission of the n-th harmful compound in exhausts e, will have the 
form: 


en f Mo, n, J,Z,G,O) 2 
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To calculate the vessel’s momentary power P* there have been used methods of calculating the 
power necessary to sail with speed v, based on calculations of resistance of sea-going vessels. 
There are a number of methods of calculating resistance of sea-going commercial vessels, among 
the most widespread being: Papmiel’s, Ayre’s, Lap’s, Taylor’s and Kabaczynski’s. The first two 
permit the determination of total (towing power) resistance, the three other only residual resis- 
tance. 

All methods (excluding Papmiel’s and Ayre’s) have vessel resistance comparable to resistance 
determined from model experiments, therefore additions have to be applied for them similar to 
those for resistance determined from model experiments. Only Papmiel’s method, which in the 
range of smaller speeds gives results considerably higher than those obtained by other methods, 
supposedly takes account of those additions. 

In connection with the amount of information available concerning the vessels’ dimensions and 
momentary speed v* obtained from AIS system [13,14], two methods were accepted for calculat- 
ing the value of momentary power P*: Papmiel’s and the Admiralty’s. 

In Papmiel’s method the power necessary for the vessel’s sailing with speed vs may be pre- 
sented by means of the equation: 


Vx v? v? 
Po caw ce Ą ję 2 
where: 
V — volume of the vessel’s underwater part [m°], 
L — vessel 's length on the waterline [m], 
vs - vessel’s speed [k], 
x — coefficient dependent on the number of shafts, taking account of protruding parts, 
A - length allowance coefficient calculated from the formula A = 0.7 + 0.3 L/100, 


y - hull fineness coefficient. 


In the case of Admiralty’s method the equation describing the power necessary for the vessel to 
sail with speed v, is significantly simplified and has the following form: 


„ph 


ÓW ak (9) 


where: 
P, — power required (towing power bound with hull resistance), 
Co — Admiralty’s coefficient for towed power, 
vs — vessel speed [k], 
D — vessel’s displacement [m°]. 


On the basis of formulae (8) and (9) there were calculated the values of momentary power P*, 
whose distributions have been presented in Figs | and 2. 

On the basis of results analysed it can be stated that the differences in calculated values of 
momentary power determined by Papmiel's and Admiralty's methods are significant. Therefore, 
particular attention should be paid to the means of determining momentary power, as successively 
the value of toxic emission for particular vessels is calculated on its basis, which entails further 
simplifications and errors. 
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Fig. 1 Distribution of momentary power P* [MW] calculated by Papmiel’s method for vessel group tested [9] 
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Fig. 2. Distribution of momentary power P* [MW] calculated by Admiralty 's method for vessel group tested [9] 


For a more accurate assessment of the emission value of particular toxic compounds there were 
worked out tests for examining the concentration of toxic compounds in the exhausts of main en- 
gines of vessels sailing in the Gulf of Gdańsk. For this purpose, out of the group of vessels sailing 
in the examined area in a period embracing 11 months 63 vessels were picked out with the highest 
incidence index, which gave 3543 observations [15]. 

On the basis of momentary power values of vessels tested P* = P,/P, 4 and 10-phase tests 
of exhaust toxicity were prepared. Seeing that most of the considered vessels sail according to 
typical propeller characteristic, modified test E3 acc. to ISO 8170 part 4 was accepted as initial 
test. The modification consisted in altering the values of weight coefficients, which in the original 
test characterise the whole range of the engine's operation, and in the case considered are to pre- 
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sent conditions of engine work on the Gulf of Gdansk fairways. Particular values of four-phase 
tests have been presented in Table 1[16]. 


In Fig. 3 there have been compared the emission values of nitrogen oxide calculated by means 
of E3 test according to ISO 8170-4 and prepared 4-phase test A (required power Po determined by 
Admiralty method) and 4-phase test P (required power Po determined by Papmiel method). 


The difference between results of NOx emission results obtained on the basis of E3 tests ac- 
cording to ISO8178-4 and 4-phase tests A and P seem to be insignificant (from 0.1 to 0.47 
g/(kW-h)). Yet if the typical powers of vessels sailing in the Gulf of Gdansk are taken into consid- 
eration, e.g. ferry “Stena Nornica” sailing from Gdynia (114 observations within 1 year, engines 
power 39600kW), then such slight differences in tests give NOx emission intensity respectively 
from 3.96 to 18,61 kg/h. In relation to NOx emission values from vehicles these are significant 
values and should not be ignored. 


Tab.1. Proposition of 4-phase exhaust toxicity tests of vessels sailing in the Gulf of Gdansk acc. to propeller charac- 
teristic taking consideration of the method of engine power determination 


Engine speed [% na] 100 | 91 | 80 63 
Load [% P,] 100 | 75 | 50 | 25 
Weight coefficient P” 0.20 | 0.38] 0.40 | 0.02 
Weight coefficient AJ 0.22 |0.20 | 0.41 | 0.17 
Weight coefficient of E3 test 

acc. to ISO 8178-4 0.20 | 0.50 | 0.15 | 0.15 


D weight coefficient P - required power Po determined by Papmiel’s method, 
? weight coefficient A - required power Po determined by Admiralty’s method. 
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Fig.3. Comparison of NO, emission values [g/(kW-h)] calculated by E3 test acc. to ISO 8170-4, 4-phase test A (re- 
quired power Py determined by Admiralty 's method) and 4-phase test P (required power Py determined by Papmiel’s 
method) 
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Recapitulation 


The modelling of emission of harmful compounds is a very important problem and a complex 
one at the same time. The scientific research work conducted at present devoted to pollution 
spreading concerns pollution in atmospheric air of stationary origin (power and industrial plants), 
motorisation and recently also airplanes. Current studies which concern first of all motorisation, 
because of, among other things, the size of marine engines, cannot find application for modelling 
emission of toxic compounds from marine engines, as the model structure depends not only on its 
destination but largely also on the amount and quality of input data. 

The possibility of obtaining data from AIS system like ship’s name, length, breadth and type, 
universal time bound with the vessel’s passing through the ,,gate”, course and velocity over the 
ground (COG, VTG) and the ship’s draft permit the creation of innovative models describing ves- 
sel traffic in the area researched and the emission of harmful compounds in exhausts both for one 
vessel and the whole researched area. 


The method of calculating power required for vessel propulsion has significant influence on the 
final results of determining values of toxic compounds emission in the exhausts of vessels. During 
analysis of data obtained from AIS system and simulated loads of main engines it was found that 
vessels mostly sail in the Gulf of Gdansk with a load of 0.4 + 0.6 P, and 0.9 + 1.0 Pa. Yet due to 
the lack of hull resistance characteristics of vessels considered it seems pointless to make tests in 
the function of vessel velocity. 

It should be added here that apart from problems that motorisation experts are coping with 
when modelling the emission of toxic compounds, in the case of vessels among the parameters 
disturbing the accurate determination of emission of particular compounds (due to lack of informa- 
tion or its variability) there should be additionally counted the engine’s technical condition, fuel 
apparatus in particular, and atmospheric conditions (particularly wind direction and force). 
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Summary 


In the real ship’s conditions, the methodology of taking an operational decision is based on the analysis of simple 
diagnostic data. The quality of a correct multi-criteria analysis of such data is the more dependant on the SDG 
diagnostic system the lower the professional competences of the engine operator. The complexity of the problem 
requires application of quick methods that analyse the information in a multifaceted way. The description of the 
current states of the engine operation can be presented in the form of a topographic map showing the current values 
of the operation parameters of the system, as scaled and correlated points on a plane. A change in the values of the 
operation parameters is imaged as a change in reciprocal location of points on the plane. On the basis of a 
topographic description, a classification of the engine operation state can be carried out by means of a neural 
classifier. Such a method of the engine operation classification will allow for obtaining a very high capability of the 
diagnostic system for adaptation and analysis of untypical states of the object and an easily construed visualisation of 
results. 


Key words: artificial intelligence, diagnostic of a ship Diesel engine, map of diagnostic parameters, diagnostic 
system. 


1. Introduction 


Different breakdowns of rates of a piston combustion engine operation can be encountered in the 
literature. 

In the operational practice, a marine engineer operates an engine that can be diagnosed with the 
methods available on a ship. The diagnostic system (SD) existing on the ship consists of the 
diagnosed system — the engine (SDN) adapted in a specified way for fitting of a diagnostic system 
(SDG) to it, consisting of sensors. The accuracy and reliability of the diagnosis depends on many 
factors, including the accuracy of selection of measurement points, the type of measured 
parameters, the measurement accuracy, the correct analysis of diagnostic data and the quality of 
the drawn conclusions — as good as the knowledge and the experience of the engine operator and 
the knowledge of an engineer who designs a diagnostic system. 
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An essential problem is the multi-criteria analysis of the measurement data and imaging of results 
of such analysis. The competences of operators are very different and they are often not able to 
correctly read and use the results prepared and given by the SDG, which influences the level of 
rationalism of their operational decisions. In the real ship’s conditions, the methodology of taking 
an operational decision is based on the analysis of simple diagnostic data that can include 
parameters of power media (most often the values of temperature and pressure) or accessible 
process parameters, e.g. the charge compression process in the cylinder or the combustion process. 
The quality of a correct, multi-criteria analysis of such data is the more dependant on the SDG 
diagnostic system the lower the professional competences of the engine operator. 


The complexity of the problem requires application of quick methods that analyse the information 
in a multifaceted way. Identification in the shortest time of processes taking place in the engine is 
very essential. Therefore the description of current states of the engine operation can be presented 
in the form of a topographic map of states that images the current values of the operation 
parameters of the system as scaled and correlated points on a plane. On the basis of the 
topographic description, a classification of the engine operation condition can be carried out by 
means of a neural classifier. Such a method of the engine operation classification will allow for 
obtaining a very high capability of the diagnostic system for adaptation and for analysis of 
untypical states of the object and an easily construed visualisation of results. 


2. Criteria structures taken for identification of the engine state 
Example criteria structures that can be taken for identification of the engine state include: 


- construction structure (different materials of specific physical and chemical properties, 
adequately fitted tribological systems, etc.), 

- organisational structure (co-operating subassemblies and assemblies of devices and the engine, 
operating within a precisely specified time, adjusted and controlled), 

- process structure (combustion process, cooling process, heating process, wear process (friction, 
oxidation, ageing, etc.). 


There is a set of couplings between the structures that can be determined in an experimental way 
but they cannot be predicted in a computational way with too many non-measurable unknowns and 
the necessity to apply many simplifications. 

Changes in the structures, in their serial dependences differ between themselves with their 
characters. 

Changes taking place in the structure of interdependencies of mechanical parts are visible and 
immediately detectable and identifiable (vibrations, noise, pressure loss of media, tightness of 
systems, e.g. hydraulic, compression pressure systems, etc.). 

Changes in the process structure are often detected with a long delay, are detectable but hardly 
identifiable, often with wrongly determined genesis and predictions (e.g. high temperature of 
exhaust gases — it may mean too low injection pressure, too large openings of an injector — but the 
reason can be different than the injector, e.g. too low fuel viscosity (too low temperature for a 
particular sort of fuel), wrong injection time, wrong ignition time, wrong time of valve opening, 
wrong process of head and cylinder liner cooling, leaky piston and cylinder system. 


The parameters of the engine operation determining the technical condition of the engine, closely 


attributed to the adequate structure and the area of the engine can be topologically imaged on the 
“map of parameters” with anomalies, possibly explicitly geometrically or colouristically deformed 
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and clearly suggesting the operational decision to the operator as an alternative to any graphs and 
fields of operation non-identifiable and non-comprehensible to the operator (with different levels 
of technical knowledge). 


Fig.l Dependence of the combustion process on the criteria structures of the engine: 
KS — combustion chamber: Tsp — combustion temperature, Psp — combustion pressure, tsp — combustion time 
D — supercharging: Tpd — supercharging air temperature, Ppd — supercharging air pressure, TSPR — turbocharger, 
UTT — liner-piston system: Ps — compression pressure, Ts — compression temperature, Tspl — temperature of exhaust 
gases, SzTł — piston speed, SzczUk — system tightness, WPtł — dimension of piston rings, WRtt — dimension of piston 
grooves, WTc— dimension of cylinder liner, WchT — cooling water for liner, WP — fuel injection: Tp —fuel 
temperature, Pp — fuel pressure, HSS — hydraulic controlling system, MSS — mechanical controlling system, STUW — 
technical condition of injection system, WK — crankshaft: Ne — shaft power, SUTr — lubrication of tribological 
systems: Tos — temperature of lubricating oil, Pos — pressure of lubricating oil 


3. Maps of diagnostic parameters of an engine 


The Fig. 2 shows the idea of creating maps of diagnostic parameters of an engine. A circle is the 
simplest geometric figure, on the basis of which the initial model has been based. It is possible to 
create a ‘set of maps’ divided into ‘packages of maps’, the quantity of which is limited only by the 
capacities of electronic equipment. Due to editorial limitations, the idea has been presented on a 
single map. 

The ‘map’ should, in a simple way, signal the trends of changes in particular structures. Centric 
circles set the rings (three outer ones separated from the inner ones by a dotted line) defining the 
values of parameters: the white ring means normal state, the broken ring (outside the white one) 
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means excess of parameters above the upper limit of normal values, the dotted ring (inside the 
white one) means excess of parameters below the lower limit of normal values. 

In order to facilitate understanding of the idea, the image area (image on the screen) has been 
divided into four zones. 


G AB | [A] 4, 


Fig. 2. Idea of creating of a map of diagnostic parameters of an engine: 


- field “I-1” excess of the upper limit of the normal parameter relating to the dimension or the 
operation of an engine element or subassembly 


- field “I-2” excess of the lower limit of the normal parameter relating to the dimension or the 
operation of an engine element or subassembly 


- field "II-b" excess of the lower limit of the normal parameter relating to the dimension or the 
operation of a device or the engine system 


z - field "II-y" excess of the upper limit of the normal parameter relating to the dimension or the 
WU: operation of a device or the engine system 


| | - field “IV-B” - a field determining the state of full efficiency of an engine 
- field “IV-C” - a field determining the state of partial efficiency of an engine 


SZER - field “IV-D”- a field determining the state of inefficiency of an engine 


La 


256 


[| - “outer circles” - a field determining normal values of parameters (green colour) 


-“outer circles” - a field determining excess of the lower value of the normal parameter (yellow 
colour) 


- “outer circles” - a field determining excess of the upper value of the normal parameter (red colour) 


> In the zone marked with “T”, excess of limits of parameters concerning particular systems or 
devices has been shown (depending on the class of the map detailing): 

- inthe field “1” visualisation of excess of the upper limit of the normal parameter has been 
shown, 
in the field “2” visualisation of excess of the lower limit of the normal parameter has been 
shown, 


> In the zone marked with “MI” and “III”, excess of limits of parameters concerning 
subassemblies or the whole engine has been shown 

After processing by the logic system, the influence of the parameter change on a subassembly or 
the whole engine can be shown: 

- in the field “II” visualisation of excess of the lower limit of the normal parameter of the 
engine operation resulting from change in one or in several parameters of the engine 
elements operation has been shown (e.g. damage to the injection system in one system), 

- in the field “IIT” visualisation of excess of the upper limit of the normal parameter of the 
engine operation has been shown, 


> In the zone marked with “IV” visualisation of the condition of the whole engine has been 
shown (“inner” rings limited from outside with a dotted line): 
- in the field “B”, the white arrow shows the field (white) of normal operation of the whole 
engine marked with green colour, 
- in the field “C”, the grey arrow shows the field (marked with a thin wavy line) of limited 
operation marked with yellow colour, 
- in the field “D”, the dark grey arrow shows the field (marked with a dense wavy line) of 
unsuitability for operation marked with red colour, 


In practice, such approach to the analysis of states of the engine operation can be accomplished 
with use of artificial neural networks. The advantages of such implementation methods are their 
high flexibility, versatility and possibility to use them for analysis of the object without the 
necessity of having the mathematical description of the object. The neural networks have got the 
ability for parallel processing of information and recognition of topographic maps reflecting the 
current state of an object under diagnosis. 


For recognition of the object’s condition one can use two types of networks: learning networks and 
self-learning networks. 

The learning networks are most often used for well defined problems where we have got the 
knowledge what the network’s response will be. For the learning networks it is possible to create a 
neural classifier of states (double-state or triple-state). However, the structure of such classifier is 
complex and the learning process is quite laborious. Input of new elements into the tested system 
requires changing of the neural network architecture and its re-teaching. 


An alternative solution is the application of the Kohonen’s neural network. Such networks have 
got the feature of topographic imaging of complex shapes, thanks to which they learn any imaging. 
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Another feature of such network is the process of self-learning consisting in competition between 
neurons. The neural network responds to the input signals. Obtaining of the final order of the 
system of neurons corresponds to the diagnostic maps of an object under testing. Input of new 
elements into the system does not require significant changes, and the neural network adapts itself 
to new tasks. A detector of the hybrid diagnostic system of a ship engine can make the neural 
classifier of the object’s conditions. The neural networks analyse signals coming from the 
diagnosed object. On the basis of the base, the expert’s system takes decisions by making use of 
the information about the current map of the object’s state gained through the neural network. 


4. Summary 


Modern ship engines need rational operation aiming at both increasing their reliability and their 
efficiency. Therefore development of diagnosing systems that make use of the most modern 
technological achievements enabling processing of data is evident. Depending on the type of 
diagnostic information and the expected results of operation of a diagnosing system, different 
mutually supplementing methods of data processing are used. 

Therefore, application of neural networks as a supplement of statistical and probabilistic methods 
used at different stages of the engine control and diagnostic processes allows for shortening the 
time for assessment of technical condition of elements, subassemblies as well as the whole engine. 
The idea of maps of diagnostic parameters presented in the paper is another trial of visualisation of 
the results of the SDG operation that allows a mechanic for determination of the condition of the 
engine in a short time as well as for monitoring the trends in changes of diagnostic parameters. At 
the same time, it may contribute to selection of an unambiguous operational decision of high 
degree of reliability worked out by the SDG to the mechanic. 
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Abstract 

Proper selection of materials for cylinder liner and piston rings, correct design of collaborating surfaces as 
well as the use of lubricating oils of better and better properties assure long and reliable operation of piston 
rings in small and medium size engines. In the case of high power engines, especially marine ones failures 
caused by incorrect collaboration of rings and liner still happen. Supply and correct distribution of lubricating 
oil over cylinder surface could be one of the causes of this phenomenon. 

Regulations introduced by classification societies make impossible research on piston ring — cylinder liner 
set of running engines. Mathematical models or simulation test benches are the possible way of carrying out 
such tests. 

This paper presents some test rigs and analytical models revealing their advantages and shortcomings. 


Keywords: combustion engine, friction losses, measuring methods 


1. Introduction 


Full information on phenomena accompanying the collaboration of engine kinematic pairs 
is available from tests carried out on running engines but such tests very often are impossible. 
For example, tests on collaboration of marine engine kinematic pairs such as piston rings and 
cylinder liner require an installation of certain type of sensors (pressure, temperature or 
distance, for instance) that means a need for changes in construction of collaborating parts. 
Such modifications on marine engines are forbidden by the classification societies because 
they can considerably affect the ship’s safety. 
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Tests on material models of piston—cylinder assembly parts could be a substitution for an 
investigation on phenomena accompanying the collaboration of these elements conducted on 
a real engine. Such models allow to reconstruct the operational conditions of frictional 
associations , geometry of mating parts contact, kinematics and dynamics of their loading as 
well as phenomena occurring by the process of lubrication. An additional advantage of such 
tests is a relatively low cost, short time needed for carrying out the investigation and no 
chance for engine failure. 

Another way to evaluate the phenomena present on engine subassemblies is a construction 
and running the mathematical models of these subassemblies. It is noteworthy that a 
mathematical model is merely an approximation and should be validated in a course of tests 
carried out on a real engine or at least on its material model. Successive chapters of this paper 
will present selected constructions of test benches where investigation on piston ring — 
cylinder liner collaboration could be performed and the most significant information about 
mathematical models will be provided. 


2. Model test benches for tests on piston ring — cylinder liner collaboration 


Tests benches where investigation on collaboration of piston cylinder group elements 
could be carried out can be classified into those entirely designed and constructed for this 
purpose and other ones constructed on the basis of a typical engine where certain parts and 
subassemblies have been modified. 

The model test benches of the first type offer wider possibilities for investigation, because 
the constructional details such as the range of operational parameters adjustment, location of 
sensors, selection of measurement devices could be foreseen at the very beginning stage of 
design. On the other hand, they are much more expensive than those using parts of a typical 
combustion engine. 

Taking into consideration methods employed for investigation one can distinguish test 
stands using optical and electrical methods. 

Test stands presented in Figs 1 and 2 are the examples for the use of optical methods in 
tests on evaluation of ring-liner collaboration quality. 
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Fig. 1. Model tests stand for optical measurement of Fig. 2. Schematic of the piston-cylinder 
oil film parameters: I — high pressure pipe assembly model, projector and observer 
supplying air to the actuator, 2 — lubricating oil (camera) positions for reflection method: 
supply pipe, 3 — joint, 4 — slider leading beam, 1 — glass cylinder, 2 — dark chamber, 
5 — air actuator, 6 — slider, 7 — glass plate, 3 — diaphragm, 4 — primary filter, 
8 — semitransparent plate, 9 — camera, 5 — impulse xenon lamp, 6 — camera box, 


10 — coherent light source (laser) [2] 260 7 — secondary filter, 8 — camera [2] 


The test stand presented in Fig. 1 has been constructed on the basis of a typical IC engine 
equipped with systems for research on oil film formed between models of ring (6) and 
cylinder (7). The model of ring fixed to a piston moves reciprocally relative to a steady glass 
plate lubricated with oil and ring pressure against the plate is exerted by a pneumatic 
servomotor (5). Through a semitransparent plate(8) and transparent one (7) a laser beam (10) 
comes to the oil layer (oil film) which makes that as a result of interaction between reflected 
and falling waves an arrangement of interference strips characteristic for a momentary oil 
layer thickness is being recorded by the camera (9). There is a possibility to reproduce the 
changes in oil film thickness throughout the entire cycle of engine operation thanks to the 
further analysis of recorded images. 

The test stand presented in Fig. 2 offers a quite different method of oil film evaluation. 
Cylinder liner made of cast iron has been substituted by a glass model (1). The cylinder 
surface is being lubricated with oil mist outflowing from slide bearings and additionally with 
oil sprayed to precisely specified regions of liner by a set of special jets (omitted in Fig. 2) 
which allows for a continuous adjustment of oil dose. Thanks to strong xenon lamps it is 
possible to observe and record (camera 8) the oil layer over the cylinder surface and to 
estimate certain oil parameters including oil layer thickness and oil film extent, in particular. 

Design of the stand presented in Fig. 3 is quite different from a typical model of engine 
piston-cylinder assembly. A flat slat (5) representing a piston ring slides over a plate (1) 
covered with lubricating oil. Differently than in real engine the oil layer is not renewable, i.e. 
fresh oil does not come to the plate in consecutive strokes. The slat is fixed in a groove cut in 
a flat model of piston (clamp 7) hanged on transducer (8). It is pressed on against the moving 
plate with a flat spring or hydraulic actuator controlled by a special circuit. It is possible to fix 
in the clamp a series of slats of profiles corresponding to the face profile of a real piston ring. 
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Fig. 3. Schematic of test stand: 1 — plate, 2 — actuator, 3 — control valve, 4 — position resistor, 5 — slat, 6 — flat 

spring, 7 — piston model, 8 — force transducer, 9 — support, 10 — slat pressure system, 11 — measuring&control 

system, 12 — hydraulic oil reservoir, 13 — oil pump, 14 — pressure accumulator, 15 — analogue control system, 
16 — digital control system [5] 


The hydraulic oil from oil reservoir (14) is being pumped to the control valve (3) by the 
oil pump (13). According to the position of control valve the oil is being pumped to the left or 
right side of the actuator plunger (2) which results in the movement of cart with the plate 
simulating cylinder liner surface. 

Either an analogue (15) or digital (16) control system can be used. The first one allows for 
a reciprocating movement according to gradually selected parameters of speed and amplitude. 
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When applying the digital control system it is possible to use a computer program which 
allows to realize an arbitrary motion of the plate including a reciprocating movement, typical 
for combustion engines. 

Due to a sophisticated measuring system (11) which is a part of the simulation stand it is 
possible to determine oil film essential parameters like thickness, pressure and friction forces 
connected with the motion of ring model. Also other sensors monitoring plate motion and 
recording its position, speed and acceleration are installed on the test stand. 

Taking into consideration features of the stand including possibility of tests on rings of big 
size and different profile geometry it can be assumed provisionally that this test stand allow to 
conduct tests on oil film parameters of the engine of high power. Next chapter will validate 
this assumption. 


3. Is it possible to test a ring-liner collaboration on a test stand? 


Both marine and generator engines belong to the group of high power engines (power 
higher than 100 kW/cyl) as well as the bigger railway engines which have the power indices 
close to the border value [3]. Among those mentioned the biggest are two stroke low speed 
marine engines. Their popularity results from the highest efficiency (nowadays beyond 50%) 
and possibility of consumption so called heavy fuel of relatively low price. 

The stroke to cylinder diameter ratio is the characteristic parameter for contemporary 
marine diesels of high power where it reaches the value of 4 (this secures continuity of oil 
film over cylinder surface thanks to the high value of the piston mean speed at low rotational 


speed). For engines of higher rotational speed the value of this parameter is far lower (see 
Table 1). 


Tab.1. Vital technical data of exemplary engines of high power [1, 3, 4] 


ine locomotive generator marine 
Technical data 12LDA28 | 214D40 L23 L21 L35MC | K80MC- 

D cylinder diameter [mm] 280 230 230 210 350 800 
S piston stroke [mm] 360 300 300 310 1050 2592 

k ratio (S/D) 1.28 1.30 1.30 1.47 3.0 4.00 
Rotational speed [rpm] 750 750 720 900 210-178 93-70 
Power per cylinder [kW/cyl] 94.7 712 130 190 650 3640 
A ratio 0.25 0.26 0.25 0.23 0.42 0.44 

Piston mean speed [m/s] 9.0 Ted 7.2 9.3 7.35 8.10 


Full reconstruction of engine kinematics on the model stand would mean the same piston 
displacement with similar speed. It means that it is impossible to fulfill conditions listed in 
Table 1 on the presented test stand [5,6]. Boundary value of speed that could be performed on 
test stand is about 12 rad/s which corresponds to the piston mean speed of about 3.8 m/s 
(when maximum is 4 m/s), i.e. the value several times lower than that on real marine engine 
(see Fig. 4). 

Course of speed vs. crank angle presented in Fig. 4 show only the maximum value of 
speed and its position expressed as corresponding angle. Presentation of speed fluctuations at 
selected points of its displacement seems to be a better form of such illustration (see Fig. 5). 
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Fig. 4. Courses of speed vs. 

crank angle for: 1 — model test 

stand 2 — the LDA28 engine, 

3 — the D40 engine, 4 — the L23 

engine, 5 — the L21 engine, 

6 — the L35MC engine 7 — the 
K80MC engine 


o 45 90 135 p[] 180 


Fig. 5. Courses of speed vs. 

displacement for: I — model test 

stand 2 — the LDA28 engine, 

3 — the D40 engine, 4 — the L23 

engine, 5 — the L21 engine, 

6 — the L35MC engine 7 — the 
K80MC engine 


0,3 0,6 0,9 s[m] 1.2 


According to the test stand technical data the range of plate stroke is several centimeters to 
less than one meter (0.8 m) which allows for reproduction of ring displacement of most of 
locomotive and generator engines but could not represent the long stroke marine engine (of 
far longer stroke). 

Digital control of plate movement makes possible any kind of displacement between both 
dead center with a constant speed of 4 m/s (horizontal line no 1 in Fig. 5). Though the ring 
movement with constant speed differs from that on a real engine it can be employed for tests 
because it facilitate the evaluation of the effect of ring selected design parameters on the 
formation of oil film. 

As it comes from the above considerations full reproduction of piston and rings movement 
(taking into account speed and range) would not be possible on the presented test stand, which 
does not mean that it is completely aimless. For the engines of stroke longer than the plate 
range another kind of tests could be performed, namely over the sections shorter than the 
stroke. On the other hand, for engines of shorter stroke tests could concern speeds lower than 
the nominal one. Thanks to the fact that both the oil layer thickness and the friction force are 
approximately proportional to the speed of ring movement, the value of these parameters can 
be estimated this way for a real engine. 

A large number of input values for which tests could be performed makes that 
investigation could be time consuming and expensive and the results obtained could be of low 
value. The reasonable solution could be development of stand mathematical model (as a 
computer program). Such creation and further application (after validation) could improve the 
effectiveness of research and carry them out within a range possible to be performed on the 
test stand. 
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4. Mathematical model of test stand 


Description of the test stand presented in chapter 2 shows that this solution does not allow 
for a full reproduction of phenomena encountered on the engines of high power (listed in 
Table 1). Beside differences resulting form different kinematics of the stand following 
shortcomings should be noted: omitting the effect of ambient pressure, low operational 
temperature (ambient temperature), ring sliding over unreconstructed oil layer, constant value 
of ring pressure against the plate. 


360 


== 
FILLET 3 


Parameters 


Calculation 
of oil film parameters 
and friction losses 


180 p [°] 180 p [2] 


Fig. 6 Block schematic diagram of test stand mathematical model [5] 


Fig. 6 shows the most important modules and mutual connections of the model. Within 
initial 180 grades of crank angle the slat encounters the fresh layer of oil on the plate which 
secures fully flooded contact with plate (at the inlet). Along with the increase in plate speed 
the pressure in oil film increases which eventually causes the rise in oil layer thickness. After 
the maximum oil layer thickness came along in the region of maximum speed, a slow 
decrease appears. Thanks to so called squeeze effect this pressure drop does not reach the 
value of zero at both dead centers. At consecutive strokes (180 to 360 degrees and beyond) 
the slat slides over the oil layer left after previous stroke. The layer thickness slowly decreases 
because the oil surplus is being swept towards the plate turning points. As a result the slat face 
area covered with oil also decreases. 

Computations carried out according to the described model give courses of oil layer 
thickness (in front, behind and under the slat), pressure distribution in oil film as well as 
courses of friction force. Results can be presented in tables or graphs and can be further 
processed. 


Fig. 7. A view of computer 

program chart — noticeable 

fragments of oil film parameter 

computations performed on 
the test stand 


0.000 100.00 12.00 100.00 12.00 


0.011 100.00 11.87 26.71 0.30 a a 0.00 0.00 
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Exemplary picture of simulation program chart with courses of minimum oil film 
thickness and the thickness of oil layer left behind the slat is presented in Fig. 7. 

Fig. 8 shows a comparison of minimal oil film thickness measurement results recorded on 
test stand (a) with results of model computations (b) carried out for input data corresponding 
to the stand operational parameters. As it can be seen maximum and minimum values of oil 
film thickness are the same, also the drop in thickness along the consecutive strokes is very 
similar. There are numerous discontinuities and waves on the course of recorded film 
thickness which result from plate shape imperfections not taken into consideration in the 
computational program. 


a) 


b) 


0 180 360 540 720 900 fi [°] 


Fig. 8. Courses of oil film minimum thickness recorded on model stand (a) and for mathematical model (b) 


5. Summary and conclusions 


Summarizing analyses presented in the previous part of this paper one can note that the 
presented test stand allows: 

— to evaluate phenomena accompanying the collaboration of piston ring and cylinder liner 
in a presence of lubricating oil, including oil film thickness and friction force, for a single 
ring (slat) or the complete set of rings, 

— to reconstruct partially the piston movement for long stroke marine engines 

- to reconstruct phenomena accompanying the piston movement for locomotive and 
generator engines, however at speeds lower than nominal ones. 

On the other hand, use of computational program allows to execute tests within the range 
unreachable for the test stand. 

At present efforts on development of the model stand are being carried out towards 
possibility of refreshing the oil layer over cylinder surface in a way similar to that applied on 
marine engines and perform the test at higher temperature of operation. 
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Abstract 


Durability is one of the most important features of components (elements and subsystems) and technical technical 
objectss as a whole. By durability we understand the capability of an item to maintain its useful properties under given 
circumstances of work and maintenance until a limit state is achieved. Borderstates are understood as states when an 
object for technical reasons, economic reasons or functional reasons (innovative) and/or for other essential reasons is 
no longer fit for use. The measure of durability of technical objects is a period of use after which a objects loses its 
serviceability according to its design. The moment a technical object loses its required properties is called the top 
border of durability. However, it is worth noting, that the moment starting from which an item is no longer fit for use 
is not always a distinct one. Therefore, a multi-aspect selection of criteria determinating the upper limit of operation 
time is the crucial issue. 


Keywords: durability, usability, maintainability 
KINDS OF DURABILITY 


A loss of required usable properties characteristics can result inter alia from technical reasons, 
economic reasons or perished innovativeness as it is present hereinafter. 

The economic durability is identified with retention usable properties until the moment 
when further maintaining of requested technical availability is found to be uneconomic. This is 
dependent on both economic and technical factors. An important measure of durability in light of 
the economic aspect is the unit price of maintaining an objects availability and the profitability of 
its use. 

The innovative durability is identified with the preservation of usable features until a new 
generation of competitive and non-compatible technical objectss appears, when at the same time a 
possible modernisation of an item is uneconomic. An important feature in terms of innovativeness 
is a prospective rise in profits due to the employment of a new generation technical objects. The 
necessity to terminate the use of a technical objects can be result from its incompatibility with the 
technospheric environment. 

The safety durability is identified with retention already existing usable properties until 
hazards occur. A need to terminate the operation of an object can result from non-fulfilment of 
objects safety requirements in the human- object-environment system . 

It is obvious that the technical durability -the physical durability to be exact - is what 
constitutes a reference base for all considerations on a lost worthiness of a given technical objects. 
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The technical durability is identified with retaining a required value of usable features until the 
possibility to maintain usable features to ensure a required level of the operational avaiability by 
means of repairs and technical services is irretrievably lost. 

The technical durability expires due to the depletion of technical expendability since there is no 
technical possiblility of upgradeability of an object . 

The measure of life (durability) of a technical objects is the service life after the termination of 
which a technical objects loses its capability of use according to its designed application. 

For human- object -environment systems one uses the term of liveliness, which applies to a 
multiple renovation (including a biological renovation) and sustaining of functions (including life 
functions) over a defined period of use under determined conditions. For the above systems also 
the notion of survival under sudden overloads (e.g. during an attack, sudden weather changes, 
hurricanes, tsunami, volcanic erosions) etc. is used. Due to this, an object loses irrevocably its 
capability of further functioning. 


LOSS OF USABLE PROPERTIES 


The technical durability is considered usually at the macroscopic level and the physical 
durability at the microscopic level. 

In case of the macroscopic level (considering constrains between objects elements and 
subassemblies which create the structure of a technical objects in the energetic meaning). On a 
microscopic level one considers interatomic constrains, which occur in material structures. The 
material of design elements is considered a system with a structure of constrains, superimposed 
onto microscopic material particles. 

Considering the durability of a design material one takes into account the laws of the physics of 
failure. Therefore the physical durability is expressed as a limit period resulting from the loss of 
the hardiness of its components against the effect of enforcing factors, in particular of working 
factors. 

The art of engineering consists, inter alia, in finding the relations between defects in the 
substance of a design material and the failure of the technical objects as a whole. 

Three destructive phenomena play the basic role: fatigue, metal corrosion and wear due to 
friction. The first process results in fractures and breaks of shaped pieces. The result of the two 
other phenomena, is a loss of the surface material layer, called wear. The results of these 
phenomena decide, inter alia, whether the usefulness of an objects can be technically restored. 

Exceeding the permissible values of usable features (states of characteristics) related to a 
correct functioning of an object, causes a variety of unwanted phenomena and events. 

Such damages may leave an object in an unwanted technical state, in which this item does not 
fulfill its required usable functions; loses the capability of repairing it and/or it does not satisfy 
safety requirements.. 

According to a thorough analysis the durability is define with two notions - strength and 
resistance (hardiness) 

The strength is a property dfine the capability of a design material subjected to allowable loads 
to maintain required usable properties. It is understood as a capability of maintaining - 
maintainability - of the internal cohesiveness of a material until permissible working loads are 
exceeded. 

One distinguishes an immediate strength and fatigue strength, i.e. a capability of a repeated 
taking on of usable loads. In the latter case, a destructive process, called a fatigue occurs and 
results in brittle ruptures, slip ruptures and split ruptures as well as in plastic and elastic 
deformations. 

The probability of occurrence of this kind of damage is a basis to determine the upper limit of 
the use-time of design elements. 
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The strength is a property, therefore a physical value, expressed in units of measure, which 
consequently can have particular numerical values. 

The resistance is a feature which define the capability of a structural material and associations 
of technical objects elements to maintain the required element characteristics and associations 
when subject to impact of external and internal forcing factors.. 

The requisite cohesiveness of a system is maintained thanks to the internal energy of cohesion 
binding material parts of a system. 

There are three basic factors, which contribute to maintaining the system’s resistance; these 
are: proper strength of an item, properly selected load acting on an item and a controlled speed of 
deterioration of objects characteristics. 

The resistance can be regarded as a non-measurable feature (a behaviour manifestation). E.g. 
technical objects can have a resistance or not. This resistance can be described only approximately 
as very low, low, high or very high etc 


PRESERVATION OF OBJECTS USABILITY 


Considering technical technical objects as consisting of assemblies, association of kinetic and 
static elements, one considers them as macro systems of a structure of constrains superimposed on 
cooperating design elements. 

As a result of destructive processes within systems of associated design elements natural and 
imparted usable constrains disappear as well as parasite constrains develop. This is the reason 
why components lose their worthiness and consequently rendering useless of the objects as a 
whole. One distinguishes two extreme states of technical objects. 

The first extreme state being the one in which an item loses its capability of processing and 
transferring energy (due to evolved additional parasitic constrains and/or disappearance of useful 
constrains). 

The second extreme technical state is a state when an item loses its capability of performing 
working moves and obtainment of a required speed of relocation — movement (due to a change of 
purposeful constrains to restricted constrains and due to an elastic deformation of constrains and 
kinematical pairs under the impact of forces and temperature). 

According to a thorough analysis the durability is described with a set of two notions 
maintainability and supportability. 

The maintainability — a feature of a machine and technical equipment, define its capability of 
performing maintenance as well as service and repair steps within a given time. 

These steps are: repairs to restore a requisite condition (state), a technical service to maintain a 
required condition and a technical inspection to determine the technical condition of a technical 
object. 

The tractability of maintenance consists in features, which enable easy services and repair of a 
technical objects as well as keeping a technical objects in a required condition by way of an 
appropriate strategy of maintenance. 

The strategy of maintaining of a demanded technical operational stand-by is a rational time 
schedule of inspections of the technical objects condition and preventing steps or upgrading steps. 
This strategy consists in a determination of such work periods of technical objects, which satisfy to 
the best advantage technical, organisational and cost limitations. 

The measure of the maintainability can be e.g. such a ‘physical value’ like performed technical 
steps or a property (a behaviour manifestation), e.g. a technical objects shows the capability of 
maintenance or repair or it does not .It can be also definite - not clearly - as a very low, low, high 
or very high etc. capability. 

The supportability — a feature expressing a serviceability and usefulness of technical 
equipment and means to maintain a machine in a required technical state. 
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As a result of very advanced technical progress, the use of machines is more and more 
dependent on professional maintenance. Performing maintenance and technical checks calls for 
modern inspection and measuring devices, availability of service with spare parts and appropriate 
consumables. Logistic service guarantees safe and effective use of technical objectss. 

A measure of the supportability can be e.g.such a physical quantity as a required diagnosing 
time, spare parts delivery time etc. or e.g. a feature of behaviour like accessibility or non- 
accessibility of services [0,1]. It can also be described approximately as a very low, low, high or 
very high professionalism of maintenance. 


CONCLUSIONS 


Therefore in order to estimate durability the following aspects have to be determined: 

e a technical aspect: Jack of technical means of restoring a required technical condition due 
to an overly dispersed energy of internal structural cohesiveness of an object 

e an economic aspect: non-profitability of further use due to an excessive rise in unit 
operating costs 

e aspects of innovation: innovative manufacture of technical objects of an increased output, 

e new properties and better compatibility as a result of new requirements for use, which are a 
consequence of a new generation of solutions. 
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Abstract 


The great number of factors generate situations that real dimensions of each machine elements are different from 
nominal ones. In the case of rolling bearings elements mentioned disagreements are small — their values are rank of 
micrometers. However, such small differences can generate significant changes of operational features of bearings. 
The analysis of influence of ball dimensional deviation on load distribution in ball bearings is presented in this paper. 
The probability of occurrence of ball diameters of the range of acceptable values is statistically determined. Random 
distribution of this important dimensions results in other than nominal level of internal stresses (in the contact area of 
balls and raceways) and therefore true values of some operational features of bearings are other — the most often 
smaller — than calculated ones. In presented analysis sequence of occurrence of different ball diameters is also taken 
into account. 


Key words: ball bearing, dimensional deviation, operational features of bearing, random character of balls diameter 
1. Preface 


A lot of factors which determine each manufacturing process cause situation which has random 
character. In this connection real dimensions of every machine elements, including rolling 
bearings, are encumbered with machining dimensional deviations. Therefore, real dimension can 
differ from nominal in some value, however, not greater than boundary value for specific accuracy 
class of machining. Consequently, in case of rolling bearings the essential factor is dimensional 
selection of their elements because it determines inner load distribution in contact area of rolling 
elements and raceways. The more uniform distribution - the achievement of calculated durability 
of bearings close to the true value is more probable. 

The statistical analysis of probability of existence in the bearings the balls’ diameters of 
specified values is presented below. The results of analysis are necessary for exact analytical 
designation of fatigue life of the rolling bearings, e.g. Applied in combustion engines. 

For standard rolling bearings analytical procedures are well-known and experimentally 
verified, instead for special bearings — this kind is analyzed in this paper — appears necessity of 
description new ones. 


2. Load distribution inside bearing 


The rolling bearings contain rolling elements (e.g.: balls, rollers, needles) of different diameters 
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(of course in the range of acceptable deviations for specific bearings’ accuracy class), just for that 
reason real load distribution, and thus also contact pressure between rolling elements and 


raceways, will be distorted in comparison to analytical — Fig. 1. 
a) b) 


E=1 Ye =180° 


P max cos a 


Fig. 1.Theoretical load distribution in angular ball bearing: a) distribution of forces in contact area of balls 
and raceways, b) distribution of internal load of bearing circuit for determined load conditions 


The values of internal load existing in the rolling bearings and stresses in the contact area of 
balls and raceways one can calculate by well-known relationships, e.g. [2, 5]. For inner forces P, it 
is the following equation: 


3 
P, =P 1- l- 0087) 0) 
i 26 
in which: 
P 
max — —— > (2 a) 
J ,(é)z sina 
or 
P 
max = = . > (2b) 
J,(6)zsina 
where: 
Pmax - maximum internal force, 
e - load distribution factor, 
y - position angle, 
P,, -load components, respectively: axial or radial, 
Jxy -load integrals, respectively: axial or radial, 


zZ - number of the balls, 
- contact angle. 
For calculation of the maximum contact stress value Omax the following relationship is used [5]: 


ELA M | G3) 


a*-b* 


where: 
a*, b* - dimensionless axis of contact elipse, respectively: semimajor 
or semiminor, 
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Sp -curvature sum. 

In consideration of the fact that loads values are function among others of rolling bearings 
dimensions (especially contact area geometric features), calculated values of inner loads existing 
in the bearing, assigned to specific angles, one should admit as approximate values. 

As a result of changed distribution, attitude towards nominal load distribution, and thereby 
other values of contact pressure, bearings’ fatigue life can be different (most often smaller) 
because amplitude values of stresses in contact area of rolling elements and bearings’ rings also 
are different (usually greater than nominal). 

In order to settle quantitative determination of that influence, below the analysis is carried out. 
To this aim theory of probability is applied [1]. Analysis is conducted on the example of angular 
ball bearing applied in the front wheel hub of a bicycle. 


3. Analysis of balls dimensions 


The accuracy of rolling bearings is determined conjointly by geometrical constructional 
features of following bearings’ elements: 

e rolling elements (balls, rollers, needle rollers, barrels, cones). 
e inner ring, 
e outer ring, 

In manufacturing consideration, conveniently is to assume that bearing element diameter (in 
analyzed example — ball) is redundant dimension. Its value is assorted to specific, measured values 
of raceways diameters of rolling bearings — on the outer and inner rings. Thanks to it, the proper 
value of clearance is obtained in the bearing. 

In the all manufacturing processes surface geometric structure (SGS), with random and 
determined components, is generated on machined elements [7]. Therefore, true dimension is 
stochastic quantity — different from nominal dimension assumed by designer. For the aim of balls 
selection, in order to limit range of these differences, in standard [6] are created nine accuracy 
classes. In each class, values of boundary deviations of: diameter, geometrical tolerance and 
surfaces’ roughness are determined. On account of diameters’ spread in the inspection lot, in 
quoted standard the division of individual accuracy classes to dimensional groups and subgroups is 
suggested. 

The rolling bearings are composed of finite, strictly determined number of rolling elements 
(balls) z. Therefore, the probability of existing of balls with determined diameters should be 
limited and refer to the lot number z. The distribution of diameter deviations in the lot of balls is 
presented in Fig. 2. 

Individual quantities in the mentioned figure are defined in standard [6]. The way of their 
values delimitation is also described in this reference. In below presented considerations the most 
fundamental of them is Vp, — scatter of diameters in the lot. It is the difference between the 
greatest DwmL max and the smallest Dwmi min average balls’ diameter in the lot. 


scatter of diameters in inspection lot Vou. 


BSK 
minimum diameter Dym min maximum diameter Dym max 


NM —— 0 —— = 


lots 


average diameter Dyn -| k _| scatter of observed diameters Vpws 


Fig. 2. The distribution of bearings’ balls dimensional deviations in inspection lot 
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Particular attention should be paid to the fact that average ball diameter Dwm is arithmetic 
average of the greatest and the smallest individual balls’ diameters. 
4. Randomness of balls’ diameters - statistical distribution 


Factor X, refers to true ball diameter, have Gaussian distribution with average value m and 
standard deviation s. This fact may be written as follows: 


X~N(m, s), (4) 
In the bearing exist the balls with the diameters in specific accuracy class, therefore new 


random variable Y is introduced. Its features limit the set above mentioned what can be note as 
follows: 


X= (5) 
0 for |X -m|>a 


This variables set has Gaussian distribution, with two-sided section, limited by assumed 
deviations values a of average balls diameters. This type of distribution is shown in Fig. 3. 


f(x) 


m-a m mta X 


Fig. 3. The graph of probability density function of random variable x 


The probability density function of analyzed factor X has alternate form: 


(x-m)? 


= for |x-m|<a 


fzo) = > (6) 


0 for Ix—m|>a 


where: 
c — constant, calculate from equation: 


sS 
in which 


(x), for x = a/s — distribution function of standard Gaussian distribution, 
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which can be described by relationship: 
2 
REZ 
= 2 
P(x) = Tr fe du | (8) 


and 
u — working variable. 

The distribution function of random variable X , in individual function intervals, assumes 
values in accordance with dependences mentioned below: 


0 if x<m-a 
F(x)= =") d 2) if m-a<x<m+a, (9) 
c sS sS 
1 if x>m+a 
Let Xin X2n, ... Xnn are ordinal statistics from random test (X, X .., Xn) arising from 


distribution with distribution function F(x). The value of average r-th ordinal statistics X,., one can 
assign from relationship: 


n- 1 m+a 

BX, =a) PAGNE -rG ar, (0) 
m-a 

In the case of analyzed real rolling bearing, parameters’ values in equation (10) are following: 

e n=z=9, 

e r=1,2,...,9. 

The statistic assumed in presented analysis: 


W, = Ama = Kia ad 1) 


is called the range of the sample and it is difference between maximum and minimum average 
diameter in the sample — see Fig. 2. For example, in 40" accuracy class of the bearings, value of this 
quantity, according to standard [6], is equal +16 um. Average value of the range of the sample, for n 
elements in the sample, can be calculated from relationship: 


EW, = ju =F"(x)=[1= F) Vie, (12) 

Analytical solution of considered problem makes possible to illustrate the dependence of EW, 
on standard deviation s. This way can be determined which part of the range (m-a, m+a) is 
covered by ordinal statistics X; 9, X29, ..., X99. 

In this aim dependence of average value of range of sample EW, on standard deviation is 
evaluated for the whole lot. For better interpretation of computation's results, values of ratio 
EW,/2a are given on the Y-axis — Fig 4. 

With the help of above described procedure, values of balls diameters are calculated, in 
assumption that machined balls are subjected to Gaussian distribution with two-sided section. In 
the calculations the true constructional features of tested bearings are taken as input data. 
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EMS 


Fig. 4. Dependence of EW, versus standard deviations (value and unit of s — adequate to basic size) 


As results of calculation the following values of balls diameters dy were obtained: 4.754, 4.757, 
4.760, 4.762, 4.762, 4.762, 4.764, 4.767, 4.770 mm. They composed input set of quantities in 
bearings’ fatigue life calculations. 


5. Probabilistic analysis of bearings fatigue life 


As earlier was found, differences between nominal and true diameters of balls in the bearings, 
which are contained inside the range of particular sample (in considered example the lot of analyzed 
size 1s z = 9 and it is equivalent to balls number in the bearing), cause that inner loads distribution, 
shown in Fig. 1, is only theoretical. The true value of amplitude, as the result of such changing 
load, will be also other than theoretical. Thus, there are good reasons to state that analytical fatigue 
life will be different from determined on experimental way. 

The sequence of the balls of differentiated diameters occuring in the bearing is also random 
factor. Differences of balls’ diameters cause differentiation of contact stress - see equations (2) and 
(3), therefore, order of balls sequence — Fig. 5, generates changes of stress amplitude — quantity 
which decides about fatigue life of bearings. 
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Amplitude 
of stresses, ©; 


Number of cycle with individual amplitudes, n; 


Amplitude 
of stresses, ©; 


Number of cycle with individual amplitudes, n; 


Amplitude 
of stresses, 6; 


Number of cycle with individual amplitudes. n; 


Fig. 5. The examples of possible cycles of load in the contact ball-raceway caused by nominal dimension deviations of 
ball diameters: seq.1) decreasing sequence, seq.2) ascending sequence,seq.3) mixed sequence; 
n — nominal level of stresses’ amplitude 


The above figures show that balls’ sequence existing in the bearing directly determines run of 
load cycle, therefore, it can have significant meaning for fatigue life of analyzed rolling bearings 
[8, 9]. 


6. Closure 


In the operational processes of all machines one observed that durability of correctly serviced 
bearings’ pairs is different (usually greater) than determined on computational way. This is caused 
by analytical methods imperfection. That is why evolution of computational methods is observed, 
e.g. [3, 4, 10]. It appears in including increasing number of the factors, which determine durability 
of bearings. 

Application of the above presented analysis in constructional practice should increase 
calculations accuracy and thereby, to induce better closing of analytical and practical results. This 
way the mathematical models, describing observed changes will be more correlated (congruent) to 
real run of researched phenomena. 
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Abstract 


It is well established that one of the major factors of accident prevention on board is the perfect theoretical and 
practical knowledge possessed by engine room officers while operating engines and auxiliary equipment. 

This paper describes the latest developments in 3D computer simulation applications, designated for the 
familiarization with marine machinery, specially taking into account the emergency operating procedures training. 

The experiences in 3D computer simulation application, the benefits and advantages of use of computer 
simulation in educational process of engine room officers in the Gdynia Maritime University are equally presented in 
this paper. 

This paper describes an example of application of new 3D simulation techniques in engine room simulator based 
on modern - computer controlled engine room with low speed main engine MAN Diesel LMC type, applied on a 
container ship, where trainees have possibilities to develop operational skills, update their know-how and refresh 
emergency procedures. 

This installation plays a vital role in the safety on board. Thus, it is of extreme importance that trainees acquire 
the capacity to react in a prompt and effective manner to emergency situations. This new 3D simulation technique 
specially emphasizes the relation between simulation and realism of machinery operation. The simulator described in 
the paper provides for a new approach to navigation through the different system’s elements, allowing for an easy and 
quick access to basic engine room operation (valve opening/closing, setting position of switches, push-buttons etc.). 
This has been possible due the application of state-of-the art 3D visualisation with zoom techniques. 

The basic tasks for computer simulation in maritime education are equally described, taking into consideration 
the new methods and emergency procedures training. 

This papers’ conclusion is that the use of 3D computer simulation in maritime education results in increased 
emergency preparedness and in consequence, leads to hazard mitigation and reduces the risk of human error in the 
operation and maintenance of marine equipment. 


Keywords: 3D computer simulation, marine engine room simulators. 


1 Introduction 


Marine engine simulators allow for operation of emergency situations that are not permissible 
under normal exploitation conditions due to safety limitations. Simulator’s software includes also 
assessment features that enable objective review of trainees acquired capacities. During simulator 
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exercise the instructor is able to apply various exercise set-ups (initial conditions) and scenarios 
that include different fault finding tasks. 

For this reason, engine room simulators are more and more used in maritime academies as a 
valuable asset for educational process [1]. The application of engine room simulators is also 
recommended by STCW 95 IMO Convention [2]. 

It is worthwhile mentioning that marine engine room simulators have also some basic 
disadvantages. Namely, they include lots of simplifications, abbreviations and schematic 
presentation of machinery systems as a result of the fact that they are presented only in 2D 
visualization. Hence, the trainee with perfect knowledge of simulator operation can experience 
serious problems with real ship power plant operation, because the graphical presentation and 
operating procedures of the simulator are distinct from the reality. 

For this reason, manufacturers of engine room simulators begin to apply 3D graphical system’s 
layout presentation in order to provide a machinery configuration as close as possible to reality. 
The main problem in creation of 3D simulators is to provide for proper navigation through the 
system’s elements [4,5,6]. Engine room is a complex, multi level and complicated set of sub- 
systems, equipment and machinery and this is a new challenge for entities creating such kind of 
simulators. 

It is also necessary to allow for an easy and quick access to basic engine room operation (valve 
opening/closing, setting position of switches, push-buttons etc.). It is possible to achieve this 
feature by applying zoom techniques for selected elements of the system. Users of 3D simulators 
should also be able to observe the system’s elements from pre-select specific parts of the engine 
room. 

Based on the author’s experiences with the application of different types of simulators, a better 
solution consists in navigation by mouse cursor and zooming facilities. 

The application of new 3D simulation techniques in marine engineering education shall be 
analysed on the examples of full mission, hardware type engine room simulator with low speed 
main engine. 

A new technique of navigation through the system’s elements has been applied in this virtual 
reality simulator, providing for a solution of the main problem in creation of 3D visualisation. The 
latest development includes also a combination of 3D and 2D diagram presentation, which enables 
to follow how a certain device really functions and gives a complete picture of its structure. The 
presented solutions have improved considerably the level of simulator fidelity in relation to real 
machinery. In consequence, it was possible to eliminate the disadvantages of the engine room 
simulator with typical 2D presentation consisting in a schematic and simplified presentation of 
machinery systems. 

The application of virtual simulation in teaching the operation of complex marine machinery 
leads to a better understanding of the functioning principles of both the equipment and the systems 
in comparison with traditional educational methods. As a result, trainees are far better prepared to 
deal with real life operation of machinery, thus increasing in a considerable manner the standards 
of safety of ship operation. 

On the other hand, software, due to its features enables the trainee to repeat in an unlimited 
number of times the required operations, thus to achieve the necessary preparedness level. 

The engine room simulator based on the medium speed engine room simulator is one of the 
first simulators which use hardware type of consoles combined with 3D visualization. 

This simulator specially enhances the operational procedures related to emergency situations, 
like electrical black-down, emergency manual operation of the main engine with propulsion 
system as well as auxiliary machinery in case of remote control failure. As it has been said before, 
these procedures may not be trained in real life conditions due to safety constraints. From didactic 
point of view the best solution is to combine hardware version of engine room simulator with 3D 
visualization. Such combination improves in a considerable manner the safe operation of marine 
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engine room, as the crew members have previously been trained in relation to various fault 
scenarios. 


2 LER3D Low speed engine room simulator’s description 


The application of new 3D simulation techniques in marine engineering education shall be 
analyzed on the examples of software type engine room simulator with low speed main engine. 
The LER3D virtual reality simulator provides for a new approach to navigation through the 
different system’s elements, allowing for an easy and quick access to basic engine room operation 
(valve opening/closing, setting position of switches, push-buttons etc.). This has been possible due 
the application of state-of-the art 3D visualisation with zoom techniques. The latest development 
includes also a combination of 3D and 2D diagram presentation, which enables to follow how a 
certain device really functions and gives a complete picture of its structure. The presented 
solutions have improved considerably the level of simulator fidelity in relation to real machinery. 
In consequence, it was possible to eliminate the disadvantages of the engine room simulator with 
typical 2D presentation consisting in a schematic and simplified presentation of machinery 
systems. 

This simulator is designated for training students of maritime academies as well as for different 
types of marine vocational training centres. The simulator has universal features and may be used 
both for training merchant and navy fleet crew. 

The main purpose of the LER3D simulator is the practical preparation of the trainee for engine 
room operation, and more particularly: 

- familiarization with the basic engine room installation (compressed air system, fresh and 
sea water cooling system, lubricating, fuel oil system etc.), specially taking into account 
training on the base of modern, computer controlled engine room; 

- acknowledgment with main engine and auxiliary equipment exploitation procedures; 

- propulsion system maneuvering ; 

The software also generates the main engine room’s sound. 
The simulator has been developed in compliance with: 

- STCW Code: Section A-1/12 and Section B-1/12. 

- ISM Code: Section 6 and Section 8. 

A general view of the LER3D simulator engine room (main engine ) has been presented on fig.1. 


Fig. 1 Low Speed Engine Room Simulator LER3D — General view of main engine 
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It is well known that one of the major factors of accident prevention on board is the perfect 
theoretical and practical knowledge possessed by engine room officers while operating engines 
and auxiliary equipment. The basic role of Low Speed Engine Room Simulator LER3D is the 
familiarization with different operational modes, required for achieving a high level of emergency 
preparedness [6]. 

The basic role of this simulator is the familiarization with different operational modes. 

This simulator allows not only for training under normal operation conditions but also 

for emergency operation procedure training. The more familiar the trainee is with the 
equipment the faster and more effective are his reactions to a state of emergency [3]. 

On fig. 2 main engine’s emergency control local station is presented. 


Fig. 2 Low Speed Engine Room Simulator LER3D - Main engine emergency local control station 


On of important element for training emergency operating procedure in LER3D simulator is 
emergency power plant. The basic role of Emergency Power Plant is to supply electric current, in 
case of failure of main diesel generators. 

General view of 3D Emergency Power Plant simulator is shown on fig.3. 
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Fig.3 Emergency Power Plant — 3D simulator’s general view 


Emergency Power Plant diesel engine view with control panel zoom selection is presented on fig. 
4. 


Fig.4 Emergency Power Plant simulator - Diesel engine with control panel 


Diesel engine starting procedure can be effectuated in two ways: 

- from control panel with application of electrical starter (24 V available) 

- by emergency hydraulic starter (24 V not available) — fig. 5 
In case of diesel engine emergency operation, the trainee ought to create proper pressure in 
hydraulic bottle and by manual lever start the engine. Emergency hydraulic engine’s starter is 
shown of fig. 5. 
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Fig.5 Emergency hydraulic starter 


One of possibilities to improve realism and fidelity of operation procedures training is application 
of touch screen. On fig. 6 use of touch screen is presented for hydraulic starter emergency 
operation. 


Fig. 6 LER3D Engine Room Simulator — Application of touch screen for hydraulic starter emergency operation 


Emergency operation procedures are emphasizes in hydraulic systems operation. 
On fig. 7 steering gear installation (rotary vane type) is presented. Trainee, by pressing push- 
buttons on solenoid control valve can change rudder position in emergency ( manual) way. 
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Fig.7 Steering gear installation - Solenoid control valve with emergency push-buttons 


Another example of engine room installation where it this possible to train emergency situation is 
fire fighting system. An example of S-type fuel oil separator fire simulation is shown on fig. 8 
LER3D Low Speed Engine Room Simulator in relation to fire fighting system consist of: 

- CO2 system 

- Water mist system (fig. 9 — 10) 
Described engine room simulator allows for very detailed operating procedure training like CO2 or 
water mist system activation. 


Fig. 8 S-type fuel oil separator fire simulation 
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Fig. 9 Water mist system — activation panel on engine control room 


Fig. 10 Water mist activation in Diesel generators section 


LER3D virtual reality engine room simulator can support and facilitate preparation of trainee for 
use of full mission, hardware type of engine room simulator. 

This software type simulator is first of all intended for individual students’ training ( one student — 
on PC). Typical configuration of PC classroom is 8 till 12 stations + instructor station. En 
examples of PC classroom with two screen solution on each station and PC projector is shown on 
fig. 11. 
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Fig.l] LER3D Engine Room Simulator — PC classroom 


3. Conclusion 


The application of engine room simulators with 3D visualization, in maritime education leads to 
a better understanding of the marine machinery and results also in increased emergency 
preparedness and in consequence, leads to hazard mitigation and reduces the risk of human error in 
the operation and maintenance of marine equipment. 

The application of simulator gives the operator of equipment the opportunity to get acquainted 
with the system’s operation, before starting real life usage. For this reason, more and more often, 
manufacturers of marine equipment supply their products accompanied with software type of 
simulators. 

As it results from the above, the latest developments in simulation techniques, including 3D 
presentation enhances the above mentioned benefits, as it brings the engine room simulators closer 
to reality. In consequence, the gap between operating marine machinery in simulation conditions 
and in real life is decreased. 

As it has been mentioned above, proper navigation through the 3D simulator’s elements is the key 
point in order to achieve the didactic purposes. The new concepts of view selection, zooming 
features of elements and operation by cursor and mouse clicking, as applied in LER3D low speed 
engine room simulator, described in this paper, appears to be very effective and easily adaptable 
by trainees in practice. 

Due to the specificity of operating marine equipment in real life conditions, the didactic goals in 
marine education are directly linked with achieving preparedness for emergency situations. Such 
preparedness may only be achieved if the trainee is familiar with both the equipment and its 
operating modes, including emergency situations. 

In the near future, this type of 3D solutions should be applied more and more often in engine room 
simulators design. The presented simulators are related to marine machinery, but the concept of 
composition and navigation through the system’s elements can be easily applied for the purposes 
of any type of technical equipment and shall contribute in a similar manner to hazard prevention. 
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To summarize, the application of engine room simulators with 3D visualization is a valuable 
hazard prevention tool as it reduces the level of human error in the operation and maintenance of 
marine equipment. 
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Abstract 


The results of investigations connected with the implementation of the procedures of monitoring the technical 
condition of machine engines and their investigation for the chosen arrangements of mechanical vehicles were 
introduced. The examples of the rule of inference of diagnostic were passed for the gathering of diagnostic 
parameters and the opinion of the state. 


Keywords: machine state recognition, procedure algorithmization, conclusion rules 
1. Introduction 


The usage in exploitation process of machines technical state evaluation methods, being the 
basis of the automation of state recognition, requires the determination of diagnostic parameters 
set, diagnostic tests assignment and method optimization of diagnosis, genesis, forecasting. The 
realization of these tasks depends on many factors related to the level of machine complexity, the 
quality of exploitation process and the course of wear process. For this purpose it is necessary to 
identify the methods of machine state recognition, create procedures possible to use in the machine 
state recognition process, as well as determine the rules of concluding being an element of the on- 
board software of machine state recognition system [2, 3, 5, 6, 10, 11]. 

Machine state recognition is a process which ought to enable: 

a) real-time estimation of the machine’s technical state on the basis of diagnostic researches 
results, through the control of state and the location of damages in case of the machine’s 
disability; 

b) anticipation of the machine’s state in the future on the basis of incomplete history of 
diagnostic researches results, which enables time estimation of the machine’s flawless 
operation or the value of work done by the machine in the future; 

c) anticipation of the machine’s state in the past on the basis of incomplete history of 
diagnostic researches results, which enables the estimation of the machine’s state or the 
value of work done by the machine in the past. 


2. Problem characteristics 


The problem of state evaluation, prognosis and genesis of the machine’s state is equally 
important at the stage of construction preparation, production and exploitation, and the main 
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problems appearing during the solution of the machines’ state recognition task are [4, 6, 7, 8, 10]: 

a) formalizing the aim of diagnosis, forecasting and genesis of the machine’s state; 

b) the description of the machine’s state changes in the exploitation time; 

c) building a diagnostic model — the relation between state characteristics and diagnostic 

parameters; 

d) the solution of state diagnosis task; 

e) the solution of state forecasting task; 

f) the solution of state genesis task. 

The machine’s state W(©,) at the moment in time ©, can be characterized with the use of the 
parameter values set {y;(©); j=1,....m}. The machine at the moment ©, (machine’s state diagnosis 
task) is in the state of aptitude S°, when satisfied is the condition: 


WO = W° <> Y GL...) [yah < ty; @v)} < {Yie}, (1) 
where: tyja). {yj,e} — sets of upper and bottom boundary values of diagnostic parameters. 


Respectively, it is possible to formulate the condition of aptitude at the moment O, 
(forecasting task) [9,10,14]: 


W(@n )= W° © V G=l....m) [fyja) < ty; (Oomi )} < fyjsh]. (2) 


whilst the elements of the set {y(@pi1 )} are unknown, hence the necessity of foreseeing them in 
the assumed time range tı. The element t; means the time range for which the forecasting process 
is realized (the element c; is also called advance or “prognosis time horizon”). In this aspect, the 
time evaluation of machine’s passage from the aptitude state into the inaptitude state is determined 
by the results of diagnostic parameters prognoses {yj(@piti)} signalizing the exceeding of 
boundary values. 

Similarly we can formulate the condition of aptitude at the moment Qy.» (machine’s state 
genesis task): 


W(On-2 ) = wey (j=l,....m) [tyja) < ty; (@v-2)} < {Yg} 1» (3) 


Whilst some elements of the set {yj(©p-.2)} are unknown, hence the necessity of foreseeing them in 
the assumed time range t2. The element t2 means the time range for which the genesis process is 
realized (the element t2 is also called advance or “genesis time horizon”). In this aspect, the 
estimation of the machine’s state or the value of work done by it in the past are determined by the 
results of diagnostic parameters values genesis ty;(Ob-x2). 

The main problems appearing at the solution of such defined tasks are: the selection of “the 
best” diagnostic parameters describing the actual state and their change in exploitation time, 
diagnostic test determination, the determination of forecasted value of diagnostic parameter for 
prognosis horizon Ti, yjp(Qpt+ti) with the help of “the best” prognosis method, and the 
determination of the next diagnosis or operation ©, the determination of genesis value of 
diagnostic parameter for genesis horizon Ta, yjp(©p-t2) with the help of “the best” Genesis method, 
and the estimation of the machine’s state, work done by it in the past, and the determination of the 
cause of an existing damage. 

The above used notion “the best” is connected with the acceptance of appropriate criteria and 
considering these problems in the categories of optimal solution search, whilst taking into 
consideration many estimation criteria, it is necessary to consider these problems in the categories 
of polioptimal solution [1, 10]. 
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3. The algorithm of machine state recognition 


The scheme of determining procedures for machine state recognition is shown in the Fig.1, 
whilst the algorithm of the use of the appointed optimal diagnostic parameters set, optimal 
diagnostic test, optimal genesis method and optimal method is shown in the Fig.2. 

Summing up the considered in this chapter problems concerning theoretical bases of the state 

recognition process methodology, it is necessary to state that: 

1. The process of machine state recognition includes the following kinds of diagnostic research: 
state evaluation, genesis and prognosis. 


Data 
j sl), (A, 0 
EG 
Y parameter determination | w» 
algorithm 


Optimal test 
determination algorithm 
Optimal operation date 
determination algorithm 
Past machine state 
estimation algorithm 


Perfonnmg recessy 
operational actions 


Fig.2. Scheme of machine state recognition realization 
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2. In the process of machine state recognition, an important element is the determination of the 
diagnostic parameters set: 

a) for the evaluation of the machine’s state considering the criterion of machine states 
differentiation; 

b) in the process of prognosis and genesis with the use of criteria of: correlation between the 
diagnostic parameter value and the state and machine’s exploitation time, and the 
information capacity of the diagnostic parameter. 

3. In the procedure of machine state evaluation, analyzed was the possibility to use the following 
methods for examining the relations between diagnostic parameters and the machine’s states: 

a) mean test examination; 

b) examining the distance of trust ranges of the diagnostic parameter mean value; 

c) examining the changeability of the diagnostic parameter value. 

4. The presented estimations allow to state that it is possibile to create appropriate algorithms of 
the machines’ technical state recognition, and appropriate concluding rules vital for creating 
the software of the on-board machine technical state recognition system. 


4. On-board state recognition system 


Created rules of diagnostic models identification, established procedures of acquisition and 
processing of information on state changes of examined models, their implementation in 
information technology, as well as proposed state recognition procedures (diagnosis, genesis, 
prognosis) are the basis of on-board diagnostic systems construction. In this chapter we presented 
the assumptions for the construction of the on-board machine state recognition system, the rules 
for its design, and the procedure of the diagnostic system construction [10, 11]. 

The idea of the on-board state recognition system ought to be included in the organizational 
range of the new machine’s design, including the problems determining: 

a) functional features; 

b) construction features; 

c) exploitation and operation conditions; 

d) operational resource base potential; 

e) economic calculation. 

Taking into consideration the above assumptions, the creation rules of the on-board machine state 
recognition system should include: 

1. The project analysis — concerns the analysis of needs and possible problem solutions, the 
analysis of economic calculation concerning the construction of the system, as well as the 
analysis of functional, technical and economic requirements including the relation of 
manufacture and exploitation costs of the system to the acquired benefits in result of its use; 

2. System design — concerns the determination of the architecture of microprocessor, modules, 
interfaces and other components, as well as software in the aspect of the system useful 
features fulfillment, including: 

a) logical project — concerns logical aspects of system, processes and information flow 

organization, 

b) functional project — concerns the function description of system elements and their co- 

operation, 

c) construction project — concerns the structure description of system elements (e.g. processor, 

memory, communication, input and output, clock, power supply); 

3. The construction of system work simulator (computer) — its purpose is to ensure the initial 
evaluation of system work with simulated machine states; 

4. System implementation — the aim is to create and build a physical model performing the 
assumed functions of the system; 
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5. System quality assurance — concerns testing the programs and exploitation research; 

6. System documentation preparation — including the construction, requirements and limitations, 
functioning and operational procedures. 

The object of the machine state recognition system working involves: 

1. In the area of State Evaluation: 


b) 


c) 


concluding on the machine state based on the relations between the measured values of 
diagnostic parameters set and their nominal values on the basis of the created state control 
test; 

concluding on the machine state location based on the relations between the measured 
values of diagnostic parameters set and their nominal values on the basis of the created 
damage location test. 


In the module State Evaluation we experience: 


a) 
b) 
c) 


diagnostic matrix creation on the basis of input data; 
the possibility to edit the diagnostic matrix; 
saving to a text file. 


2. Inthe field State Genesis based on the concluding on a possible cause of the located disability 
state of the machine according to the algorithm: 


a) 
b) 
c) 


d) 


optimal diagnostic parameter (maximum weight value); 

optimal genesis method (minimum genesis mistake value); 

minimum value of the distance between the diagnostic parameter value with the range of 
genesis mistake and the boundary value of the diagnostic parameter din; 

correlation of the minimum distance value dmin with the machine’s disability states si, with 
the simultaneous testing of circumstances and exploitation conditions of their occurrence, 
as the potential cause of the located machine 's disability state. 


In the module State Genesis we experience: 


a) 
b) 


c) 
d) 


the determination of genesis value of diagnostic parameter and genesis mistake; 

the determination of the minimum distance between the genesis value of the parameter and 
its boundary value; 

testing the influence of exploitation factors (number of parameters, row size) over the state 
genesis; 

the visualization and record of approximation or interpolation function for chosen 
parameters of a chosen object in the form of a graph and a report (Fig. 3). 
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Fig. 3. Interpolation method with the third-degree compound function from the State Genesis module (estimating the 
genesis value of the parameter NO, (weight w,=0,9229) for the interpolation with the combined function, type 3) 


3. In the area State Prognosis on the estimation of the date of the next machine operation Qa 
(accepting the minimum date value Qa, , Od) on the basis of the analysis of values of dates 
generated by two prognosis methods for the next operation date ©4ı and Op according to the 
algorithm: 

a) optimal diagnostic parameter (maximum weight value); 

b) optimal prognosis method (minimum prognosis mistake value); 

c) minimum value of the next operation date Oy. 

In the module State Prognosis we experience: 

a) the determination of the prognosis value of the diagnostic parameter and prognosis 

mistake, 

b) the determination of the diagnosis and operation date of the machine; 

c) testing the influence of exploitation factors (number of parameters, time row size, 

prognosis horizon value) over the state prognosis; 

d) the visualization and record of analyzed prognosis models for chosen parameters of the 

chosen object (Fig. 4). 

For the realization of the above system functions it is necessary to use the object programming 
procedures. Then the basic module of the machine state recognition system is the knowledge base 
<OBJECT, ATTRIBUTE (state characteristic) — VALUE>. The object defined in the state 
recognition system are the systems and units of the machine. 

The information on the machine state have hierarchical structure, where general information take 

the highest level, e.g. for the State Evaluation (machine state control), and lower levels are 

predestined for more detailed information (system damage location). 
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Fig. 4. Brown-Mayer model from the State Forecast module (Estimating the prognoses value of the parameter py. 
(weight w,=0,858) and date O,(t=AO, t=2AO, t=3 AO) for Brown—Mayer model, type 1 (a=0,2)) 


An ideal state recognition system is an on-board diagnostic system performing state control, 
damage location, state genesis and state prognosis functions. In this case the cost of the machine 
rises, however the machine’s exploitation effectiveness grows as all the state recognition functions 
are performed. Such solution of the state recognition system can be appropriate for critical 
machines or other special machines. 

A less costly solution is an on-board state recognition system which performs only state 
control functions. In this case the external state recognition system can forecast the state or locate 
the damages of the object and perform the genesis of the disability state. It can be a universal 
system used for diagnostic researches of different machines. 


5. Dedicated conclusion rules 


The analysis of research results (for example) of machine state prognosis methodology allows 
to formulate dedicated conclusion rules of type “IF — THEN” or “IF — THEN — ELSE” in the area 
of: 

a) diagnostic parameters optimization; 

b) state prognosis. 

In case of the combustion engine UTD-20, the generated rules have form: 

a) for diagnostic parameters set Y°: 

- if wy; 2 0,05 then yj € Y°, 
- orifwij= Wijmax then yj € Y°; 
b) for state prognosis: 


295 


- if Wij = Wijmax and if wyj = 0,9 then yj €Y° and the set Y° is a single-element set, 
y’ =y°! 
- if Wij = Wijmax and if wyj < 0,9 then yj €Y° and the set Y° is not a single-element set, 


- if the prognosis mistake of Holt method (with appropriate values of the parameters a, B) 
for the set Y°! < prognosis mistake of the Brown—Mayer method (with an appropriate 
value of the parameter a) for the set Y°, then the method of value prognosis of the set 
Y°! is the Holt method (with appropriate values of the parameters a, B), otherwise the 
prognosis method of the value Y°" is the Brown—Mayer method (with an appropriate 
value of the parameter a), 

- ifthe value of the next examination date of the engine UTD-20 Qa (X) < value of the 
next examination date of the engine @a (Y°'), then the method to estimate the next 
examination date of the engine is the method of levelling the prognosis mistake value, 
otherwise it is the prognosis method of diagnostic parameter boundary value, 

- if the prognosis mistakes for methods: Holt (with appropriate values of the parameters 
a, B) or Brown—Mayer (with an appropriate value of the parameter a) for diagnostic 
parameters of the set Y°° take minimum values, then prognosis methods of values of 
appropriate diagnostic parameters of the set Y°° are the above methods, 

- if the value of the next examination date of the engine UTD-20 Qa (Y°°) < value of the 
next examination date of the engine Od (Y°°) then the method to estimate the next 
examination date of the engine (for the considered diagnostic parameter) is the method 
of levelling the prognosis mistake value, otherwise it is the prognosis method of 
diagnostic parameter boundary value, 

- ifthe value of the next examination date of the engine UTD-20 Oa is determined for Y°’, 
then this values is the weighed value of the value Ogy. 

The presented conclusion rules in range of machine state prognosis, after performing 
appropriate verification researches, could be the basis for dedicated software of a machine state 
recognition system in an on-line mode (for an on-board system) and off-line (for a stationary 
system). 


6. Conclusion 


From the analysis of action requirements, and the configuration of the on-board machine state 

recognition system, it results that the architecture of the system should assure: 

a) the system configuration in the range of beforehand determined needs, including inserting 
the appropriate number of diagnostic parameters, their boundary values, and nominal 
values of diagnostic parameters, machine state, machine operational time; 

b) measuring and recording diagnostic signals measured values according to determined 
conditions (the beginning and end of the measurement, which values, and when they are 
recorded); 

c) diagnostic concluding on the basis of relation analysis between the standard values and the 
measured values on the basis of diagnostic conclusion rules analysis; 

d) machine state visualization, including exploitation decisions generation (able, unable, 
damage location, other). 

From the above it results that it becomes necessary to create a data base in which, apart from 
boundary and nominal values sets, and the set of diagnostic parameters values recorded during 
exploitation, there are diagnostic conclusion rules. 

The presented conclusion rules can be the base for programming the machine state recognition 
system in an online mode (for an on-board system) and offline mode (for a stationary system). 
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The presented algorithm of conclusion rules generation unambiguously identifies the machine’s 
set, system or the machine, in the aspect of their state recognition, which should enable the 
creation of a dedicated software (for machines’ sets and systems) of the state recognition system. 
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Abstract 

In the paper has been made an attempt to replace the conventional system of fuel dose control and injection 
passing angle with electronic control system, which has been realized this manner, that selected hydraulic 
accumulator, injector C-R type and fuel dose controller has been attached the marine engine. 

It has been assumed, that by controlling the current impulse of controller it is possible to model whichever size of 
fuel dose, injected to the combustion chamber and to control the injection lasting time. 

Fuel which feeds the engine must flow through differently formed/shaped channels of particular component 
elements of fuel system. In relation to this, there has been made an attempt of presenting the method of modelling fuel 
flow through the chosen sections of fuel system with different shapes of cross-section. 


Keywords: damping by the parallel surfaces, flow through the apertures, energetic loss 
1. Introduction 


In the construction of model describing the course of complex phenomena occurring in the 
process of fuel injection in bunker systems it is essential to consider many specific problems 
connected with the flow of this liquid through channels with complex shapes. 

In this paper there has been made an attempt of describing the fuel flow through the apertures, 
pipes with section different from the circular, also the manner of calculating the hydraulic 
resistance of two flat surfaces has been presented. 

The significant problem at estimating the operation of hydraulic appliances, which are parts of 
gear feeding the vessel engines by fuel under high pressure, is the proper estimation of tightness of 
moving connections of parts cooperating with each other. 

In the mobile joints between two cooperative parts (e.g. between piston and cylinder in the 
piston pump, between rotor and housing in the gear pump, between needle and atomizer cylinder) 
the narrow apertures always arise what is conditioned by defined constructional backlash. As far as 
the basic criterion for flows in the pipes is obtaining the lowest resistances, the aim is to obtain the 
highest resistance in the flows through apertures that guarantees possibly lowest flow delivery. 


2. Flows through the apertures between needle and atomizer cylinder 


Considering the very low size of one of cross dimensions of aperture, flow through the aperture 
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is usually a laminar motion [1,2,3]. It was experimentally affirmed that Reynolds’ critical ordinal 
for apertures (related to the minimal size of aperture), when the transition from laminar flow to 
turbulent flow occurs, ranges 

Re= YË x 600+ 1000, 

v 

where: 
Vg — mean velocity in the aperture, 
s — minimal size of aperture, 
v — kinematic viscosity coefficient. 


2.1. Frictional flow in the flat aperture at zero pressure gradient 
Let us take under consideration the flow in aperture between two parallel flat plates with 


infinite span and one of them moves to other that is immobile at constant velocity vę on the thin 
layer of liquid with thickness s (Fig.2.1.). 


Fig.1. Flow in the flat aperture caused by the motion of upperplate 


Pressure in the inlet section of such formed aperture is the same as in the outlet section (dp/dx 
= 0). Because of adhesion forces, the elements of liquid directly adjoining to plates remain in the 
release condition toward them. Whereas the indirect thin layers of liquid will translocate to each 
other, velocity gradient will be different from zero (dv/dy = 0) in normal direction to the pan 
motion. Assumption of infinite span of plates allows considering the flow in aperture as a flat 
flow. On the assumption that temperature in the layer of liquid is constant also its viscosity has not 
changed. 

The infinitesimal element of prism form with section dx -dy and unit spam dz = 1 is taken under 
consideration in this layer of liquid. From the condition of force equilibrium on the direction of 
axis x we obtain: 

p:dy:I- p-dy-1+(z+dt):dx-1-T:dx-1=0, 
and further  dr=0—Tr=C where C means constant. 
Distribution of shear stress in considered flow in the aperture is homogenous. 

Introducing Newton's formula of the dependence of shear stress on velocity gradient in normal 
direction we obtain 


T= pee C; and after merger v= m Y= Crs 


dy 
Constant C i C; are defined on the basis of marginal conditions. For the approved system 
connected with the immobile plate is 


tor y=0>v=0 V=SV=y so C,=O and C= u% 
S 


Substituting defined values to (2.1) we obtain the relation on velocity distribution 
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v= m” , so the velocity distribution in the aperture is linear distribution. Flow delivery 
S 


through the section F=s-b is equal to Q = j v: b: dy= i s:b 
0 


Vo 
2 

Previously introduced relations allow calculating the speed of shaft rotating at constant angular 
velocity around the axle which overlaps with the pan axle at presence of very small aperture in 
comparison to diameter of (Fig. 2.2) s/D<<1,where: D — diameter of shaft; s — size of radial 
aperture. 


It is obvious that mean flow velocity in the aperture Vir = 


At low relative aperture, the curvature of layer of liquid could 
be omitted and treated as flow in the flat aperture with linear 
distribution of velocity and homogenous distribution of shear stress. 
Defining the shaft length as /, moment of friction would be 


Meni 21 oh ; 
2 S 


where: vo=@-D/2 — speed of shaft, so: 
3 
Fig. 2.2. Laminar flow in the M= pon 1 2, 
ring aperture caused by the ; 4 ; ; 
speed of shaft It should be emphasized that the laminar flow in aperture 
- tween 


the shaft and co-axial pan occurs for the Reynolds’ critical ordinals Re< 30: > : 
S 


2.2. Flow in the flat aperture under the influence of pressure differences 


The second flat flow in the aperture, which has broad practical application, is the flow between 
parallel immobile plates under the influence of pressure difference. Such flows occur in the flat 
and cylindrical distributor, in the piston pump and hydraulic engines in piston-cylindrical opening 
and in the valves and other hydraulic sets. Considering very small nominal dimensions of 
apertures, flow velocities in the apertures have not reached such high values, to cause the turbulent 
flow in them. In principle, they will be the laminar flows. 

The laminar stationary liquid flow in the flat aperture between two parallel plates is considered, 
omitting the phenomena, which occur on the edges of aperture and thermal processes. 

Receiving reference standard (axis x agrees with the flow direction and places/situates in the 
longitudinal centre of aperture, axis y in perpendicular direction to the walls of aperture) was 
isolated the elementary prism with 2y x/ x I dimensions (Fig. 2.3) [1,2]. 


Fig. 2.3. Scheme of flow in the flat aperture caused by the pressure difference 


From the condition of force equilibrium influencing on this elementary prism we obtain 
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D2: y- I= py 2 y1-2-t-1-1=0 
where: 
pı and pz — pressure in sections I — I and II — II, 
T— shear stress on the lower and upper wall (considering flat flow the stress on the side walls gives 
resultant equal to zero). 


After substitution and reduction we obtain 


d 
(p,- p) y= 2-1, 


dy 
: = A A 
defining Da , dv= eR ydy, 
I I 2:u-1 
: . Ap 2 
Integral of this equation vV=-—— y+C, 
2.u:1 
From the marginal conditions yH=t 5 ; v=0, 
was determined constant C of above equation 
ae ae 
2.u-1 4° 


The formula for velocity distribution in the aperture will take the final form 


A —— 
2.u:1 |4 
Velocity distribution in the aperture is parabolic. 
Velocity reaches the maximal value in the longitudinal centre of aperture 


Ap > 
Vinax PW > 
Mean velocity of the aperture 
2 _ Ap P 
Vir > 73" Vows 12.1 S > 
Flow delivery through the aperture with unit breadth 
_ 4p; 
ESA 
and for the aperture with breadth, b>>s 
Ap 3 
Q= bua 


Pressure decease at velocity vsr 


Ap 


2 


12: 1-Q_ 12-V- Ly Vig 
bs Bs 

Flow loss through the flat aperture amounts to 

1201 vy 
gs 

or in analogous form to the formula for flow loss in the pipelines 


1 j 


j y 
AN EN 
D, 2:g 
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where: 
D, = 2:s — hydraulic diameter, 
A = 96/Re — flow loss coefficient for the aperture, 
Re = (vy--Dy)/v— Reynolds’ critical ordinal. 
Power lost at the flow through aperture is equal to work which is necessary to translocate the 
liquid through aperture at pressure difference equal to pressure loss 
s (Ap) b 
2-u 1l 
Above considerations and introduced formulas are correct on the assumption of constant 
viscosity in aperture. In reality/as a matter of fact the viscosity of liquid depends on temperature 
and pressure, which change along the aperture. For practical equations with sufficient accuracy is 
accepting the mean arithmetical value of kinematic viscosity coefficient 


N,= Q: Ap= 


> 


_Vi ty 
Vir 7° 

where: v; and v2 define values of kinematic viscosity coefficient which matches temperature in the 
inlet/outlet of aperture. 

The most common flow in the hydraulic appliances is the flow occurring between two parallel 
walls, which is caused by the pressure difference at simultaneous parallel translocation of walls. 

The examples could be: the flow between rotating toothed—wheel rim and gear casing and the 
flow between cylinder and atomizer needle at their mutual relative motion in the engine atomizers. 

If it is considered the laminar flow between two flat parallel walls and one of them moves with 
constant velocity vo at pressure gradient different from zero (dp/dx)=(Ap/l)z0, it could be regarded 
as total of two flows: fractious flow with linear velocity distribution and flow caused by the 
pressure difference with parabolic velocity distribution. 

When wall moves according to the direction of flow presented in Fig 2.4., velocity distribution 
of flow in the aperture to the co-ordinate system in the longitudinal centre plate is analytically 


defined by the formula: 
y 1| Apes |, 49 
v= = y'|—+—]+ -| Z ; 
Ky) „| A TA 


Fractious flow 


Pressurized flow 


p. p. 


Fig. 2.4. Velocity distribution in the flat aperture at pressure difference and upper plate moving according to the 
direction of pressurized flow 


If the wall moves in the opposite direction to the flow caused by the differential pressure, 
velocity in the cross-section of aperture will be 


Vi A) ZAPEZA 
v= = |>+—|+ | Z ; 
W) „| Ą | F 
Velocity distribution in this case presents Fig. 2.5. 
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Knowing the velocity distribution, the value of mean velocity and the flow delivery is easy to 
define 


Y Fractious flow 


Pressurized flow Resultant flow 


P.>P,; 
dp/dx<0 


Fig.2,5. Velocity distributionin the total pressurized flow in the aperture at motion of upper plate opposite to the 
direction of pressurized flow 


_ Ap-s +% 
Vir = a ee bj 
l2:u:l 2 
sign ,, + ” — when plate moves in direction of lower pressure, 
sign ,, — ” — when the direction of plate is opposite. 


The relation defines flow delivery through the aperture with span b 


Q= Ap s |, -b, 
(2-6-7 2 


(sign ,, +” like at defining mean velocity). 
2.3. Hydraulic loss in the pipes with section different from circular. Hydraulic radius 


The loss in the rectiaxial pipes with circular cross has been previously considered [1,2,3]. 

There should be put a question to what degree the presented considerations and formulas could 
be correct in solving the practical tasks in case of pipes with section different from the circular and 
pipes with circular section. 

At calculating the hydraulic loss, the significant source of energy loss is shear stress appeared 
on the pipe walls, which was caused by liquid viscosity. Loss caused by the liquid friction will be 
the higher, the higher is circumference of washed by the liquid in relation to pipe section. At the 
same delivery and the same section, the flow resistance in the pipe with rectangular section will be 
higher than flow resistance and it in turn will be higher than the resistance in pipes with circular 
section. It is because circle has the lowest relation of circumference to section among all plane 
figures. Pipes with circular sections are the most beneficial taking under friction loss. So it is not 
difficult to draw to conclusion that measure, which is characteristic of values of loss occurring in 
pipes with different sections, is the relation of liquid section F, to wetted circumference /,,, i.e. to 
the circumference where liquid meets with walls. 


s 


This relation is called the hydraulic diameter r,= 


Fig. 2.6 presents examples of hydraulic diameters for cases, which are most common in 
practice. 
For the flat apertures, where the width of aperture a, is significantly higher in comparison to 


high s, the values of hydraulic radius n= 


304 


Although there can not be the dynamic similarity of flows in the pipes with circular and non- 
circular-section, still on the basis of experiments it could be assumed that the character of variation 
of coefficient 2 in pipes with different sections will be analogous at the same Reynolds’ critical 
ordinals. 

1 2 
f° d _F, aa _ ab 

Th Th Th Th 2-(a+6)’ 


z:d 4 Le 4a 4 

Especially at lack of experimental data for not very accurate calculations of flow loss in the 
pipes with sections different from the circular could be used formulas given for circular flow 
without fear of making glaring mistakes. 


ODE 


Fig. 2.6. The examples of hydraulic radiuses for the most common pipe sections 


Because hydraulic diameter for the circle-section 7, = d/4, formula on calculating the hydraulic 
I Vir 
4r, 2-8 
where loss coefficient A is defined for the particular zones on the basis of formulas for sections 
with circular-section, where Reynolds’ critical ordinal is following 
dą sy, 
Re = Th" Vir , 
v 
2.4. Hydraulic resistance of two parallel round surfaces immerge in liquid and approaching 


loss for pipes with any section will has following form h=A 


In the present point the motion of two parallel plates immersed in liquid and approaching to 
each other will be considered, what concerns the sections of motion of atomizer needle in extreme 
positions [4]. 

Two parallel flat round plates with diameter R are located one above other in low distance from 
each other; space between them is filled with liquid. By pushing out the liquid the plates approach 
to each other at constant velocity u. Task lies in calculating the resistance, which plates meet. [4]. 

To calculate above task the installed set of parallel round plates was oriented to cylindrical 
coordinates with the beginning of set in the middle of bottom plate (which is assumed as 
immobile). The motion of liquid is axially symmetrical and because the liquid layer is very thin, in 
fact it is the motion along the radius (v„<<v,), where — << — ; 

z 


So the motion equation has form: 


Ev @P _o 1 ó(ry) Ów 


=—, a =0 b 
í ZZ Oo 74 (3 r ór 74 ©) 
with following marginal conditions: 
v-= v,= 0, when z=0; v-=0, v=—u,when z=h; p= p, when r=R 


where 7 — dynamic viscosity, h is the distance between plates and pressure p, is the exterior 
pressure. From the equations (a) we obtain: 
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1 4 7 
Merging equation (b) after dz we obtain u=—- > [T vy dz= eS (r L) 
r 0 


from which p=p,t m s (R = P), 
Total resistance force influencing on the plate is equalto F= di - - z R : 
. A 


4. Conclusions 


Simulation calculations for the chosen sections of feeding system has been presented on the 
basis of theoretical and empiric relations without verification and tests on the real feeding 
installation in fuel system C — R. type. 

As it results from the initial observations of measure results on the stand, the above-proposed 
relations are the output basis for detailed analysis of phenomena, which occur in the selected 
sections of fuel gear. 

To verify the approved model, the experimental tests should be conducted on the test stand, 
which is replica of feeding system1SB and the obtained result should be used to adjust the 
approved calculation model. 
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Abstract 


The paper deals with the problem of influence of changes variable stator vanes axial compressor settings of gas 
turbine engine on work parameters of compressor and engine. Incorrect operation of change setting system 
of variable vanes could make unstable work of compressor and engine. This paper presents theoretical analysis 
of situation described above and results of own research done on real engine. On the base of results of experiment there 
were determined mathematical equations determining relationships of changes of particular engine work parameters 
in function of variable inlet guide stator vanes setting angle. There are presented results of the solution 
of mathematical equations, which describe the changes of engine work parameters values too. 


Keywords: gas turbine, axial compressor, variable stator vanes, modelling 


Parameters, abbreviations and subscripts 


on - air stream outlet angle with stator vanes, 

QKW - setting angle of variable stator vanes, 

Bu fy  - air stream inlet and outlet angles in rotor vanes, 
Cla - axial component of air stream absolute speed on rotor blades inlet, 
Cim C2,  - circumferential components of air stream absolute velocity on the inlet and outlet rotor blades, 
co - combustor, 

HPC - high pressure compressor, 

HPT - high pressure turbine, 

UK - compressor efficiency, 

i - air stream inlet angle on rotor blades, 

LPC - low pressure compressor, 

LPT - low pressure turbine, 

ŁK - variable stator vane, 

m - air mass flow, 

n - compresssor rotor speed, 

0 - angular velocity, 

Pfuel - fuel pressure, 

Boo - nominal engine power, 

PT - power turbine, 

T's - compression ratio, 

u - circumferential speed, 
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Wj, Wa -air stream relative speed on inlet and outlet rotor blades, 


VIGV -variable inlet guide vanes, 

Ac, Aw, - air stream whirl in the rotor, 

VIGV -variable inlet guide vanes, 

z - number of inlet guide stator vanes, 


1. Introduction and purpose of research 


When in the compressor construction is assembled system of setting change of variable stator 
vanes its task is to make optimal cooperation engine units during the permanent improvement 
of compressor characteristic. Perturbations in the operation of this system could cause changes 
in work of compressor and engine similarly as in the case of changes caused by changes 
of rotational speed or polluted interblades ducts of compressor. 

Compressor stage unitary work on radius is defined on the base of equation of angular 
momentum and it has form 


lst = OT: (Coy — Chuy) = U-Ac, = UAW, (1) 


where: 
©0 — angular velocity, u — tangential velocity, r — rotor radius, 
Clw C2u — circumferential components of air stream absolute velocity on the inlet and outlet 
rotor blades on radius 7, 
Ac,, Aw, — air stream whirl in the rotor. 


That work is constant on whole depth of rotor blade. The sum of works is the unitary work 
of stage [2]. Involved change of variable stator vanes angle setting at a constant level rotational 
velocity (constant w) caused change of air stream inlet angle in rotor vane /, (Fig. 1). It caused 
change of axial component of air stream absolute velocity on inlet cza what is equivalent with 
the change of air mass flow m and change of air stream whirl Aw, in rotor. It influences 
on efficiency and work of stage. 

The purpose of investigations, which was carried out on real engine was determination of 
influence of incorrect operation of axial compressor inlet guide variable stator vanes control 
system of gas turbine engine on compressor and engine work parameters. 

Compressor characteristic is a relationship between compression ratio z’s, compressor 
efficiency 7's and air flow mass m and compressor rotational velocity n. It makes possible 
to determine the best condition of compressor and another engine units mating. The characteristic 
is used to select optimal conditions of air flow regulation and assessment of operational factors on 
compressor parameters. 

Therefore compressor should be so controlled in operational range of rotational velocity that 
the compressor and engine mating line has a stock of stable work. The main rule of compressor 
control during the change of their rotational velocity or flow intensity is to keep up the stream inlet 
angles i value near zero. One of the most popular ways of axial compressor control is changing 
their flow duct geometry by application of inlet guide stator vanes or variable stator vanes 
of several first compressor stages [2]. 

This solution makes it possible to change of air stream inlet angle on rotor blades of 
compressor stages by change of stator vanes setting angles during the change of compressor 
rotational velocity. Fig. 1 illustrates the rule of regulation of variable stator vanes. 

For average values of operational range of compressor rotor speed is situation in Fig. 1b — 
speed values and directions with subscript number 1. In this situation is intermediate angle setting 
of stator vanes. Air stream inlet angle on rotor blades do not cause disturbance of stream flow by 
interblades ducts. For lower values of compressor rotor speed and in consequence lower values 
of absolute axial component velocity cja, it is necessary to reduce the stream outlet angle 
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of variable stator vanes a, (Fig. la). The angle reduction range should allow keeping the same 
value of stream inlet angle on rotor blades. Analogical situation takes place during the work 
of compressor with higher rotational speed. For higher rotational speed absolute axial component 
speed cj,” increases. In this situation for keeping stable work of compressor and in consequence 
constant value of stream inlet angle on rotor blades, it is necessary to increase the stream outlet 
angle of variable stator vanes — Fig. lc. Application in gas turbine engine construction of control 
system of flow ducts geometry has a bearing on a run of unstable processes. 


a) i b) c) d) 


Fig. 1. Essence of control of compressor s axial stage by changing the setting angle of stator vanes ring at 
changeable air flow velocity; a) decreased axial velocity, b) analytical axial velocity, 
c) increased axial velocity, d) schema of flow round of axial compressor rotor blades during constant rotor speed 
and constant air stream inlet angles; k — variable stator vanes ring, w — rotor vanes ring 


2. Object of research 


The object of research is type DR 77 marine gas turbine engine, which is a part of power 
transmission system of a warship. It is three-shaft engine with can-ring-type combustor chamber 
and reversible power turbine (Fig. 2). 
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Fig. 2. Block diagram of DR77 gas turbine engine 
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In compressor construction configuration of this engine there are used inlet guide stator vanes 
which make possibilities to change the setting angle incidance (change of compressor flow duct 
geometry) in dependance on engine load. This process is operated by control system which 
working medium is compressed air received from last stage of high pressure compressor. On Fig. 3 
is presented block diagram of flow control signal of variable stator vanes system. 


-= vicy | I 
I 
MOVING RING| I 
I 
ro ror rl 
ŁK|[£K] [EK] [eK] i 
1 2]|..||z I 
I 


Fig. 3. Block diagram of stator vanes change setting mechanism 
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Compressed air from the last stage of high pressure compressor is supplied to working space 
of control actuator by cleaning and cooling block. Compressed air exerts pressure on control 
actuator elements. It causes moving of control piston which is connected with moving ring. This 
ring moves on circumference of compressor body. Ring is connected with stator vanes by levers. 
When the ring is moving stator vanes realize rotational motion changing the air stream outlet angle 
a, (Fig. 1). In cleaning and cooling block are made holes. During research air stream was bleeded 
by the holes and less air was supplied to the actuator. It caused change of setting angle axw 
of variable stator vanes. In consequence of that change flow duct geometry was changed. 

The experiment was carried out on an engine load 0,5Piom with taking into consideration 
atmospheric conditions. For this load setting angle of variable vanes has value axw = - 4°. During 
change engine load in the whole range from idle to full load setting angle axw of variable vanes 
changes in range from -18° to +18°. During experiment a few parameters of engine work was 
measured and registered. It was made for three different setting angle @xw of variable vanes: 
A-akw = -4°, B- akw = -11°, C- ægw = -18°. Tab. 1 presents measured and registered parameters. 


Tab. 1. Parameters of engine DR77 work measured during research 


Parameter Measurement range Parameter name 
NLPC 0 — 20000 low pressure rotor speed 
Nypc 0 — 22000 high pressure rotor speed 
power turbine rotor speed 
REECE pressure on compressor inlet 


pressure on low pressure compressor outlet 
air pressure on high pressure compressor outlet 
fuel pressure before injectors 
-203 - 453 air temperature on compressor inlet 
273 - 1273 exhaust gases temperature on inlet power turbine 


3. Results of research 


Change angle vanes setting from position A to position C caused the increase of air flow 
resistance by stator vanes. In consequence of that subatmospheric pressure on the compressor inlet 
pı decreases. It causes pressure decrease in next parts of compressor and engine flow duct. In this 
way reduced air density flowing by compressor, for stable quantity of stream fule supllied to 
combustor, causes increase of compressors rotor speed. The most noticeable is increase of low 
pressure compressor rotor speed caused by directly influence on this compressor incorrectly 
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setting variable stator vanes. Gasodynamical connection between the low pressure compressor and 
the high pressure compressor absorbs disturbances work of low pressure compressor which are 
transferred on high pressure compressor. Therefore range of change high pressure compressor 
rotor speed is lower than low pressure compressor. In this experiment it is below 1% and it is in 
measuring error of sensor range. 

Change of subatmospheric pressure is above 5% undisturbed value of this parameter. Changes 
of low and high pressure compressor outlet presure are adequately above 1,3% and above 2,4% 
undisturbed value of angle setting aw = - 4°. 

Changes of pressure and air mass flow intensity values accompanied disturbed work 
of compressor, during constant fuel mass flow intensity in combustor, caused enrichment of fuel 
mixture. As a result of that, temperature combustor outlet gases increases. In experiment was 
confirmed the tendency changes of gases tempertaure values even though the range of those 
changes is in measuring error of sensor range. 

On the base of results of experiment there were determined the mathematical equations modelling 
the changes of particular engine work parameters in the function of variable inlet guide stator 
vanes setting angle Qky : 


n snc — 0,7449a Kw 2 F 2,6020kw F 9234,5 (2) 
N swc = 0,0204a kw 2 = 1,12240akw + 12598 (3) 
pi =-10°axw? - 10°axw + 0,0077 (4) 
p21 = 10" agw +0,0029a Kw + 2,9814 (5) 
pa = 2:10 akw” + 0,01430kw + 8,1771 (6) 
Tp =0,0204axw * + 0,16330xw + 526,33 (7) 


Fig. 4 presents results of mathematical modelling of engine work parameters. Modelling was 
cary out an state engine load what was equivalent unchangable fuel mass flow. In this case range 
of change of variable inlet guide stator vanes setting angle aky was widen from -18° to +18". 
Researches in range @xw from -4° to +18° were not possilble to realize on real engine. It is caused 
by technical restrictions on the engine. 
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Fig. 4. Change of values of engine work parameters in function of variable inlet guide stator vanes setting angle 
gotten during mathematical simulation 
4. Conclusions 


Change of values of variable inlet guide stator vanes setting angle axw from -4° to +18° caused 
the increase of stream outlet angle of variable stator guide vanes a, (Fig. 1). It decreases air flow 
drag on low pressure compressor inlet that caused decrease of subatmospheric pressure. During 
keeping the constant engine load (constant fuel mass flow) absolute axial component velocity cja 
increases. It exerts an influence on air mass flow m increase. Simultaneously the absolute axial 
component velocity cja increase caused decrease of air stream whirl in rotor Aw,. In consequence 
of that low pressure compressor rotor speed increases (Fig. 4a). In connection with decrease of 
subatmospheric pressure the increase of air pressure on low pressure outlet compressor is caused 
(Fig. 4d). In spite of the slight decrease of high pressure compressor rotor speed the increase of air 
pressure on low pressure outlet compressor involves the increases of air pressure on high pressure 
outlet compressor (Fig. 4e). This slight decrease of high pressure compressor rotor speed caused 
increase of gases flow drag in the next gas turbine engine units for the combustor. The effect of 
above is a slight increase of exhaust gas temperature on power turbine inlet. 

Multi-shaft construction of gas turbine engine reduces effects of incorrectly setting of variable 
vanes. Therefore compressors of three-shaft gas turbine engine do not require variable stators 
vanes as many stages as compressor of two-shaft engine with the same achievements. 
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Preliminary research confirms the necessity of inspection the correct operation of variable 
stator vanes system control. It makes possibility of elimination this factor from group of factors 
informing about technical state of engine which are identified during the diagnostic inspections. 
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Abstract 


Contemporary empirical researches on the object, which is combustion engine, are proceeded basing on the 
theory of experiment. Available software applications to analyze the experimental data commonly use the multiple 
regression model, which enables studying effects and interactions between input values of the model and single output 
variable. Using multi-equational models gives free hand at analyzing measurement results because it enables analysis 
of effects and interaction of many output variables. In this paper author presents advantages of using the 
multidimensional regression model on example of researches conducted on engine test stand. 
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1. Introduction 


Contemporary empirical researches on the object, which is combustion engine, are conducted 
basing on the theory of experiment. Basic purpose of such proceeded scientific researches is to 
prove the relation between the input signals (introduced by researcher) and output signals 
(observed by him). The final purpose of statistic analysis of the results of measurements is to 
determinate functions of the test object and the empirical model of functional engine. Related to 
this task very wide equations using calculus of probability, theory of stochastic processes and 
differential calculus are very time-consuming and without computer technique and specialist 
software are practically impossible. In the process of solving problems concerning 
interrelationships and complementary approximation issues, statistic correlations, relevancy 
assessment and inaccuracy of measurements and also adequacy of functions of test object with 
mathematical issues and graphic defining special points are used available computer programs, 
among all package STATISTICA PL. Mentioned programs are based on the analysis of variation 
and they assure: 

- defining coefficients of function regression in the test object; 

- assessing principal effects and interactions; 

- defining correlations of input and output values; 

- assessing adequacy of function in the test object (empirical functional model); 

- defining mathematical dependences assumed by the operator, which results from the 

elaborated functional model. 
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It should be emphasized that the statistic computer analysis could relate to various models 
which do not concern interactions at various stages of complication, which are accepted in the 
description of model of input values. Simultaneously there is possibility of rejecting (disregarding) 
in the statistic analysis both at choice selected input values describing the object and various types 
of interactions. It means that choice of the appropriate (most adequate) model depends on the 
operator and his knowledge about the specialist theoretical basis of the researched issue. 

At assumed lower accuracy of projecting reality and in practice in case of defining the 
character of changes (trend) of output values there is possibility of significant simplification of 
approximating polynomials by concerning only input values and only these of interactions which 
are statistically significant. Also the stage of approximating polynomials decides about the stage of 
model complexity and complication of basic values. So it is justifiable pursuing to creating models 
with the possibly simple form and the most profitable linear models. It is supposed that with 
regard for making some mistakes it is better to describe researched issue with non-linear character 
in small linear sections than in single complex non-linear complete section. 

Commonly used software securing the experiment planning and its further analysis do not give 
free hand in analyzing collected material but they use prepared above presented schemes of 
analysis. So interfering in the program (program package) is impossible. Noticeable in the recent 
period development of social, medical and economic sciences caused rapid progress in using 
statistic methods which secure planning the experiment [3,9,10]. Especially econometria has great 
achievements in this field and new approach to the statistic analysis could be successfully used in 
technical researches [2]. Among all using the multi-equational models gives possibilities to study 
the correlations between input and output values and additionally concern feedbacks between 
output variables and gives possibility of their direct analysis. Such assumption in contrary to 
commonly used multiple regression is closer to real conditions even if considering the Diesel 
dilemma i.e. dependence between CO, HC and NOx. concentration. 

Below are presented results of researches on the fuel feeding system of the engine (injection 
system) using the double-value fractious plan and multi-equational model. 


2. Researches on the fuel feeding system of the engine using double-value fractious plan 


The object of research was fuel feeding system in the single-cylinder test engine 1SB installed 
in the Exploitation Laboratory of Shipping Power Stations in the Naval Academy [11]. 

To identify the influence of technical condition of engine on the energetic parameters of 
engine, there were defined sets of the input values (given parameters) and output values (observed 
parameters). 


1. Set of the input values X: 
- xı- rotational engine speed n [rev/min]; 
- X2—engine torque Ty [Nm]; 
- x3—leak of the cylinder-injection pump piston set Spy /, um]; 
-  x,— leak of the discharge valve of pulverizer needle S., [um]; 
-  xs5— leak of the skirt of pulverizer needle S;/, um]; 
- x6- leak of the needle cone in the pulverizer setting S, /, um]; 
- x7—erosive wear of the pulverizer nozzle Se [um]; 
- xa — coking of the pulverizer nozzle S; /, um] A 
-  x9—strain injector spring AP /MPa]. 


2. Set of the output values Y: 
- yı- fuel injection advance angle aww / OWK]; 
- y2—fuel injection angle a„/ OWK]; 
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- y3— injector opening pressure Powr [MPa]; 

- y4-maximal fuel injection (forcing) pressure Pytmax) [MPa]; 

- y5- speed of pressure accumulation in the cylinder (Ap/Aa); [MPa’ OWK]; 

- _y6— speed of pressure accumulation in the injection conduit (Ap/Aa) [MPa OWK]; 
- y7— fuel consumption per hour B /g/h] 

- yg— outlet exhaust temperature from the cylinder T,; /K/; 

- y— mean indicated pressure pmi [MPa]; 

- yzo— compression pressure during the fuel injection pe [MPa]; 

- yıı-highest compression pressure Peymax) [MPa]; 

- yız maximal combustion pressure Pmax [MPa]; 

- y;3— angle at the moment of maximal combustion pressure pmax [ OWK]; 

- yı4— carbon monoxide concentration in the outlet exhaust manifold Ccow [ppm]; 
- yı5— carbon monoxide concentration in the crankcase Cco [ppm]; 

- yı hydrocarbon concentration in the outlet exhaust manifold Cuca) [ppm]; 

- yz7. hydrocarbon concentration in the crankcase Cyc [ppm]; 

- yış- nitric oxide concentration in the outlet exhaust manifold Cyoxv [ppm]; 

-  yqo- nitric oxide concentration in the crankcase Cyoxs) [ppm]; 

- 29— oxygen concentration in the exhausts Co2 [%]; 

-  2,~air-excess coefficient 2. 


For this present paper analysis of the input values were limited to seven values (v;4 Vis, V16 
VI1Z V18 V20, V21). 

As a result of conducted analysis there double-value fractional plan with possibly highest 
resolution (R = III) and maximal incomplication of interaction of values describing the functional 
empirical model of fuel feeding system of engine was worked out to realize laboratory artificial 
tests. Presented in the beginning remarks of general nature and results of analyzed own researches 
justified linear model [7,11] assumed in the tests of fuel feeding system of engine. Therefore 
statistic analysis is limited to justification of chosen model, defining approximating polynomials 
and also characteristic and estimation using possible statistic measures. 

As a result of conducted analysis was accepted model considering double-factor interactions, 
in which the highest values of determination coefficient R? = I and total of rest MS = 0 
characterize all polynomials approximating output values. Values of these measures indicate that 
the model is according to the theory of experiment the accepted model is most adequate. 
Considerably lower values of determination coefficient and significant total of rest are 
characteristic of the model without interactions, and what is more these differences are depend on 
defined output values. 

Graphic confirmation of rightness of decision can be also exceptional diagrams in Fig. 1. 

Determined approximating polynomials allow to define any dependences between individual 
variables and also to calculate and estimate the influence of introduced (simulated) failures (wear) 
of the elements in fuel injection equipment on the work and toxicity indicators of engine [3,8]. It is 
assumed that it is possible to define relations (correlations) between the parameters of structure 
and exhaust toxicity indicators directly or indirectly using the engine operating indicators. It is 
assumed that this way it would be possible to select diagnostic parameters of defined elements or 
sets of engine fuel equipment among the exhaust compounds. 


3. Researches on the fuel feeding system of engine using the multi-equation models 
As it was mentioned above, commonly available models which were used to test analysis and 


based on an analysis of the multiple regression have not given possibility of research in the model 
of input variable (1) connections [2]. Fundamental feature of the models with correlative equations 
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is the fact that they allow the existence of feedback between the input variables, which is certainly 
a real issue. 

Below will be presented theoretic basis of the multi-equational models and its practical use on 
example of the presented in advance plan of experiment. 
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Fig. 1. Dependences between approximated and measured values in the model with interactions 
fuel consumption, b) concentration of carbon oxides; c) concentration of hydrocarbons; d) concentration of nitric 
oxides in the outlet exhaust manifold 
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Dependence between the input signals x,,x,,---,x,, and the output signals y,,y,,***,Vy 


could be described by means of the linear system of equations 
Yi "de + bzy: FA to + Bi Ju + Ay + yy FayX, ++ Ay Xy + 
Vo = Dy VY + bzy: + dy ++ biu Vy + dy FAX FAX; ++ Ay Xy + 67 
Ys = Dy Vi + Dy Vy + Dy Vg Ht Dau Vy + dzy + X, FAX ++ dzyXy +Ë; 


Ju = OW Fb FF bym ma + Ao + dyr + dy ++ AyyXy + Sar 
where: 
y,,¢=1,2,...,M@ - explained variables (output), 
x;,j=1.2,...,N, - explanatory variables (input), 
b, —is a coefficient appearing in i- this equation at j - this 
explained (output)variable, i, j =1,2,---,M 

a,,- is a coefficient appearing in i- this equation at j - this 

explanatory (input) variable, 7=1,2,---,N, j=0,1,---,N, 
&,- is unobserved random element in i -this equation. 

System of equations (1) can be described in the matrix form 


BY = AX+6 
where: 

1 -bp -bim ao dn" Aim 
B= — by, 1 s =p ! Ka dy Ana `° Am | 
— by — byo 1 ayo dm a yn 

Yı Xi Ši 

y= Va | Za Xz = > 

yM Xy Šk 


(1) 


(2) 


Identification of the configuration (1) means the issue of selecting coefficients of the system of 
equations (1) at values of input signals known from the measurements on the real object of input 


values Vivo Yayo mvs valz K.: 


therefore defined from the measurements X,,,%,,:*:,%y, V=1,2,-:-,K values of input variables 


X,,X,-..,X, Would be presented in the form of matrix of measurements values of input signals 


1 Xii ... AM 
~ |l % «+ X 
12 N2 
X=]... . ARE 
1 Xk `U Xk 


where, similarly to above presented plan of experiment the individual input values are placed in 


sequence in the columns: 

xo — indeterminate element, 

x, — rotational engine speed n, 

x2 — load of engine at torque Ty, 

x3 — wear of the cooperative surfaces of cylinder and fuel pump piston Sw, 
x4—loss of leakproofness of the pumping valve Szr 

x; — wear of the skirt of pulverizer needle S; 


319 


x6 — wear of the taper part of sealing needle in the pulverizer setting S,, 
x7 — wear of the pulverizer nozzle Se, 
xs — coking of the pulverizer nozzle S;, 
xo —power loss of strain of the pulverizer spring Ap. 
However, recorded values jj „Y,,,***',Vy, V=L2,:--,K output signals y,, y,,---, Yy were 


written in the form of matrix of measurement values of output signals 


Ni Voy = Jm 
r.|e a w Tio] 
Vix Vax a Ja 


where individual output values are placed in sequence in the columns: 
yı — air-excess coefficient A, 
y2 —carbon oxide concentration in the outlet manifold Ccok, 
y3 — carbon oxide concentration in the crankcase Ccos, 
y4 — hydrocarbon concentration in the outlet manifold CHck, 
ys — hydrocarbon concentration in the crankcase Cucs, 
Yo — nitric oxide concentration in the outlet manifold Cyoxk, 
y7 — oxygen concentration in the outlet manifold Coz. 
System of equations (1) can be written in the reduced form by multiplying equation (2) by 
matrix B~ opposite to matrix B assuming that its determinant is different from zero, then 
B'BY=(B'AX)+B'Ę, 


from this 
Y=B"'AX+B 6. 
Marking 
Ii=BA, WBE (3) 
following form of the reduced model is obtained: 
Y=IIX+q (4) 
where: 
Tio Tu "U Tiy nı 
II - Tg Ta TaN Lae No 
Tuo Tur U Tun Nu 
Reduced model in the form of system of equation is following: 
si = Tyg ETXA FF My ty HN 
Vo = Tyg F My FF My yXy +N (5) 
Ju = Tuot t+ yy Xy + Na 


Coefficients T,,7T,,,...,Ty, 1= 12,...,M from the above system of equation were selected 
such way to function 


k 

Jio Ty Ty) = DZ + np + Tp, +... + Ny Xy — vay » i=12,...M (6) 
v=l 

reaching the minimum, and problem of selecting best model from the range of equations (5) in 

sense of minimizing the quality identification indicators (6) was solved using the theorem of 

orthogonal projection [1], so the optimal coefficients 


te, i=1,2,...,M, j=0,,...,N 


ij? 
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reduced form of model (5) can be defined from the identity 
LU” =(XX)'XTY (7) 
where matrix of coefficients is following: 


0 0 0 
Tio Ax" Tuo 
0 0 
AT Ti Ta Tim 
II . 2 
0 0 0 
Tiy Tay U T 


The next step of analysis is defining the matrix of variation and co-variation of estimator II, 


examining the confidence range for individual coefficient of multiple regression, , which is a 


coefficient of matrix II, defining the coefficient of multiply correlation R [1,5,6]. 

All of those variables which possibly could have influence on the forming variable value 
Y should be considered while constructing the regression model. Not every of these variables are 
significant in the model. To verify which of input variables have not significant influence on the 
output variables Y the relevancy test should be use to every of obtained coefficients of model at 
individual variables. This test allows to verify the hypothesis that value of regression coefficient is 
zero. Only after rejecting such hypothesis we can claim that specific variable is significant in the 
linear regression model. Variables, at which regression coefficients are not significantly different 
from zero, should be removed from the model and construct model with lower amount of 
explanatory variables [2]. 

Final multi-equational model in which all coefficient are significant and it can be used in 
practice e.g. for diagnostic purposes we obtain only at the second or third stage, what is more at 
every stage coefficients of multi-equational regression of individual model are estimated, their 
statistic relevance and removes variables with regression coefficients insignificant from zero [2]. 

In considered case statistic ¢ has arrangement ¢ — Student at K-N-1=20-9-1=10 freedom stages 
and relevance indicator a = 0,1, and read off the schedule table t — Student's crucial value tę = 
3,169 [1]. After conducting the series of test it is seems that the criteria are too severe. It was 
experimentally certify that lower values t, significantly approach the model. 

After verifying the relevancy of its parameters and rejecting insignificant values as a result, it 
comes to the considerable simplification of models. The result of analysis is that following 
dependences occur [11]: 


Yı =f(%1,%2) 
J =f (%p,%7) 
V3 =f (X25 Xg) 


Ja = J0,X4,X5,x5) 
Js =J(22,X3,X4,X5,X5,X7) 
Vo =J(24,X5,X6,X7) 
Ja = f (X1; X2,X5) 
In accordance with (1) an equation y1 describing change of air-excess A has form: 
Vy =D Vy FAY: FbyVą +DisYs + byGYę + byy + dy tax FApX. 


After reducing (considering relevance of coefficients) the equation was accepted 
Yi = Op V2 + Oy V7 + dg FdpX,. 
Differences between the measured parameters and parameters obtained from the multi- 
equational model y,-j,, ve 1,20 are insignificant. It is confirmed by the courses of dispersion 


compared in the Fig. 2. The result of their analysis is that adjusting obtained model to values 
obtained in the experiment are significant. 
In case of other models the equations of output variable are following: 


321 


J2 = by) + by3 Vs + bąJą + b25 Y5 + b26 Y6 + b27 Y7 + a20 + 414 + 42% 

J3 = 539 Yo + B35 Y5 + b37 Y7 + 438 Xg 

Ją = b41 Ji + bys Ys + bąq Vq + 240 + 444 X4 + 445X5 + 246 X6 

Js = bsą y2 + bs3 Y3 + bsą V4 + bs7 ¥7 + aso 

Y6 = bo + 467 ¥7 + 260 + 464X4 + 265X5 + 467 X7 

Yı =bn + baja + 270 + 42% 
Analyzing equations obtained as a result of model equation researches it should be state that: 
in case of first equation i.e. dependence describing changes of air-excess coefficient /, it 
depends on the changes of load of torque Tą and along with them appear changes of carbon 
oxide concentration in the outlet manifold and changes of nitric oxide concentration. Relation 
between changes of CO and NO; is accurate, because in both case amount of oxygen decides 
about emission of these compound. In case unsupercharged engine of Increase of load also 
causes the decrease of amount of fresh load in the cylinder. 
in case of carbon oxide concentration in the outlet manifold Ccox (y2) the equation combines 
dependence between carbon oxide concentration in the crankcase Ccos, hydrocarbon 
concentration in the outlet manifold Cyc, and in the crankcase Cyc;, nitric oxide concentration 
Cwoxk, OXygen concentration Co? and two input parameters i.e. rotational engine speed n i and 
load T4. Analyzing values of the coefficients the most significant are coefficient which are 
responsible for the amount of oxygen in the combustion chamber. The most insignificant is 
coefficient which forms value of NOx. 
carbon oxide concentration in the crankcase Ccos (v3) is described by the changes of carbon 
oxide concentration in the outlet manifold Cco, nitric oxide concentration Cyoxk and input 
value i.e.: coking pulverizer nozzle S;. The greatest influence has variable Cyosxx, and it results 
from the previous researches, because the linear dependence between emission of Ccos and 
Cnoxk is noticeable. 
in case of concentration of hydrocarbons in the outlet manifold Cyc, (74) the equation 
combines dependence between air-excess coefficient A, concentration of hydrocarbon in the 
crankcase Cpcs, concentration of nitric oxides Cyo,, and input values, i.e. loss of 
leakproofness in the pumping valve S; wear of the skirt of pulverizer needle S; and wear of 
taper sealing part of pulverizer needle in the injector setting S,. The most significant is air- 
excess coefficient A. Also significant is variable Cyoxx though it has opposite direction to 2. At 
the similar level of impact remains input variable S; i.e. wear of the skirt of pulverizer needle. 
concentration of hydrocarbons in the crankcase Cyc, (ys) is described by the output variables, 
first of all concentration of hydrocarbons in the outlet manifold Cyc, which as it was known 
from the previous researches changes in proportion to concentration of Cyc,, at similar level 
of impact remain Ccox and Cyoxt., although they are inversely proportional to Cyoxs. 
changes of concentration of nitric oxides Cyox (vo) describes firt of all the dependence of air- 
excess coefficient A, i.e. factor which is directly responsible for forming the nitric oxides. At 
the lower level is variable Coz. Changes of concentration of nitric oxides Cyo,, are also 
described by input variable S, i.e. wear of the skirt of pulverizer needle. 
changes of oxygen concentration in the outlet manifold Coz (v7) are described, similarly to the 
changes of A from the first equation, first of all air-excess coefficient, from he first equation 
first of all air-excess coefficient, oppositely correlated nitric oxide concentration in the outlet 
manifold Ccox. Similarly correlated is also input variable which comes from engine load by 
torque Tig. 
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* 
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Fig. 2. Dependences between the approximated and measured values in the multi-equational model 
a) concentration of carbon monoxide in the manifold; b) concentration of carbon monoxide in the crankcase; c) 
concentration of hydrocarbons in the manifold; d) concentration of hydrocarbon in the crankcase 
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4. Conclusions 


Presented description of the active experiment space by the multidimensional models gives 
great possibilities in analysis of measurement data and scientific conclusions. Furthermore, 
assuming that coefficients’ matrix II’ is orthogonal, there is a possibility of fulfilling reverse task, 
that is assessing, with complex relevance at known input variables, which describe work point i.e. 
engine rotational speed n and torque load Ty, the other input values. In the nearest future authors 
will work on this issue. 
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Abstract 


In the thesis some problems concerning fabrication of composites based on present knowledge and author’s 
own experience have been revealed. The methods of fabrication of the composites take into consideration following 
operations like size-reduction, agglomeration and compounding (mixing components) in solid (grainy) state to obtain 
the composite of determined properties assigned for processing operations like injection moulding, extrusion or press 
moulding. 


Keywords: grainy materials, rubber powder, mixing, composites 


1. Introduction 


Production of materials from polymer about required properties joins mostly with selection of 
suitable polymer materials (more and more often of composite materials — obtained from waste 
material “thermoplastic materials” and waste material “rubber”) also selection of proper methods 
of processing [1-5]. Processing of composite materials will demand complex interlocking in 
technological process of each composition process [4-6]. The principle problem from scientific 
point of view and utilitarian are prospecting of conditions connected with constitution required 
proprieties of composite materials also qualification of conditions processing these composites in 
aspect their uses in practice [7-9]. 


2. Aim of this paper 


Aim of this work is present idea permiting on qualification of conditions constitution of 
composite materials in technological process, aim of receipt of products from materials about 
given properties and definite qualities. One accepted argument, that exists one profitable harvest of 
composition process constituting properties needs will be this injection, extrusion or pressing) of 
composite materials (produced from waste material foil” and rubbers) and that are special 
manners of their realization (for example characteristic conditions of process agglomerating-plant 
of foil, of crumbling for rubbers and mixing) which proper selection will be decide about qualities 
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of products received from these materials. Basing on the literature study, preliminary research in 
industry [10,11] and a long-term international co-operation [6,12-14] and especially with the 
Technical University of Chemnitz/Germany, composition processes were developed (see Table 1). 
These are: 

- Process of crumbling of input — materials (connected from them such problems how 
obtainment reproducible of fraction, estimation of size and of shape and of external surface 
of elementary grains), 

- Process of mixing of input — materials (in solid state constants - about figure of grains) and 
of additional components (often also about grainy figure) and estimation of degree level of 
mixed composite materials with methods direct or indirect, 

- Process of processing principle obtained of composite materials (in dependences from 
needs will be this injection, extrusion or pressing). 

In Table 1 (column 1) one seized in synthetic manner these composition process, which in 
essential manner can influence on proprieties of useful composite materials. From table results, 
that condition rational elaborations in wanted solutions is nearer recognition replaced of 
composition process, which qualitatively would qualify influence select factors constructional- 
technological (appointed in column 3 embracing for example kind, dimensions of material 
crumbled and mixed), described by functions of object of researches, on physical sizes qualifying 
effectivity suitable composition process (appointed in table 1 in column 2, as measure). Will 
demand this so uses and leadership proper researches methods. 


Tab.1. Summary composition of composition process constituting properties of composite materials 
[4,6,7,9,13,14] 


=" Technological- 
Measure qualifying ; 
po: SE ae constructional factors 
Name of composition process effectivity of composition : SĘ 
essential for composition 
process 
process 
1 2 3 

A Agglomerating-plant OBI orgue Mo; M,, Fe = ffa, a2 a 
force of cut: F, Pe sie hs ? 

: torque Mo, Mo, Fe = fai, a........... an) 

IB Crumbling OB II farce of cut Ë; 

C Mixing OB III power B nme tss. ih =, b) 


level of mixed Mm 


D Process of processing 
(pressing, injection, 
extrusion) 


pressure pp, Po bp = NOn 62 evasinsi Cn) 
OB IV time t, of pressing 


Tensile strength Rm, 


IE Verification of researches | OB V level of mixed. Mu, 


F Other OB N cece een e ence nen eeeeaees YP. @O = f (x1, OE AN. Xn) 


3. Experimental part 
To reach founded aim of work one passed research laboratory basing on program of 


researches, using following materials 
- from recycling (agglomerate) polythene (PE-LD) about grain class equal 1.0 mm: 
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crumbled waste material rubber (about nearing proprieties to practical Tyres of vehicles of 
rubber type SBR) about grain classes given below. 

On suitable position research [12, 14], one crumbled waste material of rubber, with aim of 
obtainment of suitable grain classes. In result passed strainer analysises in (accordance to PN- 
71/C-04501) one obtained rubber powder about following grain classes: 0.2; 0.4; 0.6 and 0.8 mm. 
On Fig. 1 one showed of for example distribution curve grain for different classes grain rubbers. 
btained fractions of rubber powder one added (in masses) according to material from recycling in 


following quantities: 5, 10, 15, 20 and 25%. 
a) 


FRITSCH PARTICLE SIZER ANALYSETTE ZZ 


506 
zize (microns) 


Prhekorper Nr 3 


b) 


PRITSCH PARTICLE SIZER ANALYSETTE 22 


size (microns) 


Prbekorper Nr 4 
Fig. 1. Example- distribution curve for grain classes 0.4 (and) and 0.6 mm( b) 
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Process mixings of components in state solid and estimation of quality of process mixing one 
realized according to provided for procedures in description of position research [13, 14], with 
utilization of computer programmer COREL PHOTO-PAIN. Received samples at random 
according to literature [9], one subjected to estimation of degree confusions( with immediate 
method) showed proprieties of mixture nearing two component to homogeneous [10,11]. From 
replaced received at random samples composite one executed with method of injection and of 
pressing pressure — according to C (tab. 1) samples about standardized shapes and dimensions to 
researches strength with aim of verification of mechanical propriety according to D (tab. 1). 
Results of researches of influence suitable grain classes (0.2; 0.4; 0.6; 0.8) and contents of rubber 
powder V[%] in recycling PE-LD on endurance on extension one represented graphically on Fig. 2 
and 3. 


T particel size 
5 pin mm 
= 17 
E 
a 
S 
o 
= 
2 
D 
7) 
|- 
2 

0 5 10 15 20 25 

value of rubber V [%] 


Fig. 2. Influence contents of rubber powder V[%] and sizes of grains( p = 0.2; 0.4; 0.6 and 0.8 mm) on endurance on 
extension Rm [MPa] — for pressing 


particel size 
pin mm 


tensil strength Rm [MPa] 


0 5 10 15 20 25 
value of rubber V [%] 


Fig. 3. Influence contents of rubber powder V[%] and sizes of grains(Q = 0.2; 0.4; 0.6 and 0.8 mm) on endurance on 
extension Rm [MPa] — for injections 
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Analysing nearer courses of each function one can notice, that essential influence on research 
feature in section from 0 to 25% has manner of preparation of samples, of melting’s of crumbling 
and content of powder in recycling. prepared Samples by pressing show considerably higher 
values than injected. However fall of strength is not similar- identical for all of grain classes. 
Interesting is that least fall one can notice for composite in all section of content powder rubber for 
grain classes equal 0.2 mm. 

This will create possibility of minimalizing fall of strength by using of rubber powder about 
dimensions below 0.1 mm. Use in future researches and uses of rubber powder on level sub- 
micro(10°°) or most profitably on level nana(10”) will permit probably to obtain composite about 
unparalleled to this times proprieties chemical-physics [2,6-10]. As goes for fall of tensile strength 
for remaining grain classes (composite with participation of powder about dimensions 0.4; 0.6 and 
0.8 mm) this is more considerable to greater in content. On all observed microstructure, in this also 
for example shown on Fig. 4 one can notice comparatively little adhesion among warp (PE-LD) 
and with interpolations (grains of rubber) - refers this especially microstructure shown on fig. 4c 
and of fig. 4d. 


rubber powder about class 0.4 mm, b- 20% and 0.4 mm, c — 10% and 0.8, d— 20% and 0.8 mm 


This observation will demand of however further researches especially in range of measurement of 
adhesion on border: polymer material — grain of powder. Research should be so realized with 
regard of wider section of classes grain and greater contents of rubber powder 


4. Conclusions 


There were no procedures to analyze the fabrication of composites made of thermoplastic film 
wastes and rubber wastes in connection with properties of components in solid (grainy) state 
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before compounding and press moulding by now. Because of specific basic operations the tests 
were carried art on special test stands using phenomenological methods. The results have enabled 
determination of the most advantageous directions for designing technological lines within the 
range of utilization of selected plastics and rubber wastes. The results might also be used for 
designing other constructional-technological solutions which would concern fabrication of 
composites based on other secondary polymeric materials. 
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Abstract 
The article presents methods of estimating the availability of vehicles supported by numerical examples. All the cases 
concern the vehicles used In air base logistic system. Markov chains and processes hale been used to work out the 
mathematical model of their usage. 


Key words: availability, usage, service 


1. Introduction 


The mathematical model has been built on the basis of analyzing states collection for 
vehicles used in air base logistic system. Each of the technical objects can be at certain time ¢ in 
one of extinguished states, which make numerable (finite) set of states. The usage is understood as 
the movement of the vehicle on the extinguished states collection. 

The model of vehicle usage as a random process X(t) with finite set of states. X(1)= S; means 
that in time ¢ the analyzed vehicle is in a state S;. The realization of the process is meant as a 
sequence of extinguished states and their durations. The sequence, their durations and frequency 
are dependent first of all on work organization, the types of vehicles, the structure of subsystems 
collaborating in usage process [7]. 


2. The description of usage states for vehicles 


Used vehicles can be in different states whose number is limited. Three models of 
emergency vehicles in air bases A, B and C of different states collection have been analyzed. 
Model A (Fig.1) does not include the system of vehicle restoration meant as the replacement of 
used vehicle with the new one (in working order). 
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S, — performing a task 

S2- being on duty (waiting for a task) 
S3— everyday service 

S4- renovating 

Fig.1. Transition graph of vehicle usage process — model „A ” Ss -1-st level service 

S6- 2-nd level service 

S7— periodic service 

Ss - vehicle supply 

So - diagnostics 


Unknown theoretical probabilities of transition are estimated by empirical probabilities in a 
year time — transition frequencies for a vehicles as: 


0, =n,/n, (1) 
where: 


n, = Yin, - the number of transitions from S; state; 
j 
nj— the numbers of transitions from state i to j in a period of time. 


Stochastic process X(t) in a continuous time is ergodic, if at least one positive boundary 
probability of finding a vehicle in state S; for t>co exist, which is called ergodic probability p;: 


too 


p, =limP(X(t e<t;t+At>)=S,); p, 20; DP; =1 (2) 
J 


where: 
p;— boundary probabilities; 
P(X (t e<t;t+At>)=S ,)- the probability of being in state j in time interval <¢, t+ At> fora 


vehicle. 


Boundary (ergodic) probabilities fulfil standardizing conditions which means that at least 
one of them is positive. In relation to Markov processes it was proved [1, 2, 3] that if boundary 


probabilities exist, they can be calculated from boundary matrix in n steps M, =M}? . In other 


words the linear equation or equivalent matrix equation must be solved, i.e. coming from 
continuous time f to discrete time n being of number of further experiment of observing the vehicle 
in time Af: 


AP, =limp,(n) => pp; © Mi |p, |=[p,], przy Yip, =1 (3) 
i j 
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where: 
M A - transposed transition matrix My; 


[pi] — vector of ergodic probabilities; 
pi — probability of transition from state i to state j. 


Giving consideration to a condition: M” -[p, |=[p, | and Xp, =1, for A model 
j 


boundary probabilities must be calculated by doing the following systems of equations shown in 
this matrix form: 


0 0 p, 0 0 0 0 p, 0 Pi Pi 

0 0 0 0 0 0 0 Py Po || Po b, 

0 p, 0 0 0 0 0 p, 0 P, P, 

0 0 0 0 0 0 0 Pu Por | | Pa P4 

0 0 0 0 0 0 O ps 0 [p| ps |=| Ps (4) 
0 0 0 0 0 0 0 Pe 0 Ps Ps 

0 0 0 0 0 0 0 Ps P» || Pr P7 
ps 0 0 0 0 0 0 0 pos | | Ps Ps 

| Pio 0 0 Py Ps Pos Pr Peo 0 | LPs] LPo | 


dp) =! (5) 


Model B of a vehicle maintenance description additionally takes into account the 
replacement state S79 (Fig.2), which is the reflection state with the same probability transitions S79 
(P910 = P109) and zero p1010=0. 


S, — performing a task 
S» — being on duty (waiting for a task) 


S; — everyday service 
l S4- renovating 

Ss -l-st level service 

S6- 2-nd level service 

S; — periodic service 


Ss - vehicle supply 
So— diagnostics 
Sio- the vehicle replacement 


Fig. 2. The graph of vehicle transitions with replacement 
possibility Sio 


Also for B model the following systems of equations for ergodic probabilities p; have been 
estimated: 
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0 0 0 0 0 0 0 ps Po 0 ||P Pı 
0 0 p, 0 © 0 © 0 0 olip P, 
p, 0 0 0.0 0 0 0 0 0 Ps P; 
0 0 0 0 0 0 0 0 py 0 |] p, P4 
0 0 0 0 0 0 0 0 py 0 |) Ps |_] Ps (6) 
0 0 0 0 0 0 0 0 p, 0 ||» Ps 
0 0 0 0 0 0 0 0 p 0 ||p, P: 
Pa Pa Ps Ps Ps Po Po 0 Po 0 || Ps Ps 
0 Pa 0 Poy 0 0 0 Pog 0 Pow Po Po 
6.6 0 6 6 0 © © pw 0 | pul [Po] 
10 
DP, =! (7) 


jal 


In model C undetermined (universal) state S,,, (Fig. 3) has been taken into account. It 
separates the states from S7 to S; and has zero probability of return (p7;,1;=0; pi1j20; pi1i20 for 
ij#l1). S,, state is meant as time loses of a vehicle being of each of the remaining states resulting 
from organizational reasons (technical etc.). in a real time S7; is parallel, but in the graph it has a 


terraced character. f 
Sı — performing a task 


S> — being on duty (waiting for a task) 
S; — everyday service 

S4- renovation 

S; -1-st level service 

S6- 2-nd level service 

S7 — periodic service 

Ss - vehicle supply 

S9— diagnostics 

Sio- the vehicle replacement 

S,1- undetermined (universal) state 


Fig.3. The graph of vehicle transitions. The vehicle renovation and replacement SA well as finite times of 
transition between states are covered 


Ergodic probabilities p; can be calculated in accordance with (3). Suitable systems of 
equations for C model have been written in matrix form: 
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0 0 0 0 0 0 0 Pig Pig 0 Pin Pi 
0 0 Px» 0 0 0 0 0 0 0 Pou || Po 
Px 0 0 0 0 0 0 0 0 0 Pau P3 
0 0 0 0 0 0 0 0 P% 0 Pan Pa 
0 0 0 0 0 0 0 0 P59 0 Psat Ps 
0 0 0 0 0 0 0 0 Poo 0 Pou || Ps |=0 (8) 
0 0 0 0 0 0 0 0 Pr 0 Pu || Po 
Pu Puz Psa Psa  Pss Pe Py 0 Px 0 Pau Ps 
0 Po 0 Pos 0 0 0 Pos 0 Pojo Pon Po 
0 0 0 0 0 0 0 0. Pros 0 Pros | | Pro 
[Pui Puna Pus Pus Pus Pus Pur Pus Pus Puno 0 |[Pu| 


DP; =! (9) 
j=l 
After having solved the system of equations (8) i (9), taking into account the empiric 
probabilities for presented models A, B and C in discrete time the following results have been 
achieved (Fig. 4): 


$1 S2 53 S4 S5 S6 S7 S8 s9 S10 S11 
Fig. 4. Boundary probabilities of Markov chain for A, B and C models 


model C 
model B 
model A 


Ergodic probabilities pj(t=n) for S; state estimated from Markovs chains are the 
probabilities of input for each state. For example probability of duty p3=0,23 means that the limit 
of duty state number is 23% of all vehicles states number in a considered period of time. It does 
not mean that the vehicle was on duty in 23% in average in 7,. However probabilities p;(n) for 
Markov chains are ergodic and applicable to state collection process but not a physical time. 
Therefore they cannot be interpreted in a real time. Only after the standardizing of the transition 
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matrix, which means taking into account the duration time for each states (transitions from discrete 
times to real one by means of transition intensity A, ) achieved results will show the real situation 


of vehicle usage. 
3. A method to determine functional availability of vehicles 


To estimate ergodic probabilities for continuous time in compliance of Markov processes 
matrix equations are fulfilled [2]: 


(A")[pj] = 0 (10) 


and 
be p,=l (11) 
j 
where: 
A=[A,j] — intensity matrix with diagonal elements -A;j and Aj. 


For A model the ergodic probabilities satisfy the equation system with its standardization 
condition (12): 


-A,, 0 Ay 0 0 0 0 Ay 0 Pi 
0 -4, 0 0 0 0 0 Ap Aa || Do 
0 Ay hg 0 0 0 0 Ag 0 B: 
0 0 0 -A,, 0 0 0 Az. Ai P, 

0 0 0 0 -4, 0 0 Ag; 0 |-|p; |=0 
0 0 0 0 0 -A, 0 Az 0 Ps 
0 0 0 0 0 0 -A;n 49 0 P; 
Ais 0 0 0 0 0 O —Ag Ag Ps 
L Ar 0 0 Quy Ås o În fy Ass | | Po 

9 
> p,=1 (12) 


For B model the ergodic probabilities satisfy the equation system with its standardization 
condition (13): 
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= Aia 0 Asi 0 0 0 0 As 0 0 Pi 
0 -—%, 0 0 0 0 0 Ag Ap 0 PD» 
0 Ay -Åy 0 0 0 0 Ag; 0 0 D3 
0 0 0 -A,, 0 0 0 Aa Aa 0 Ps 
0 0 0 0 -A,; 0 0 Ags 0 0 Ps|_o 
0 0 0 0 0 -Aw 0 Ags 0 0 Pe 
0 0 0 0 0 0 -A,, 4% 0 0 P> 
Arg 0 0 0 0 0 0 Aggy Ag 0 Ps 
Any 0 0 Ag Asy Ao Ang Ag Ao Aiow Po 
L 0 0 0 0 0 0 0 0 Agi — ioio | | Pio 
10 
Xp =! (13) 


|-4, 0 Aai 0 0 0 0 Ag, 0 0 Ay | | pi | 
0 =4, 0 0 0 0 0 Ap Ap 0 0 Po 
0 Ay -Ay 0 0 0 0 As 0 0 A, || Ps 
0 0 0 — dą 0 0 Az Noa 0 Aan P4 
0 0 0 0 — Ass 0 0 Ags 0 0 Asi Ps 
0 0 0 0 0 — Keg 0 zę 0 0 Ae p, |=0 
0 0 0 0 0 0 Ay Ag 0 0 Any Py 
Arg 0 0 0 0 0 0 —Agg Avg 0 Agi Ps 
Aro 0 0 Ag Aso As Arg Ags Avy Aros Ao Po 
0 0 0 0 0 0 0 0 Asio Aiow Mon Pio 
L Aa 0 Ans Ana Ans Ane dar Ans mo Ano 74un | LPu 
ul 
pa (14) 


jal 


After having solved the above-mentioned system of equations the boundary probabilities of 
the Markov process are obtained p;(t) for a vehicle (if they exist). If we put the transitions 
intensities to a system of equations (14) and are not able to solve it that means that the whole 
process is non-ergodic (not periodic, not regular). One thing we can do is simulation, which shows 
the irregularity of the process. Moreover it allows to estimate the level of non-ergodicity. Another 
way [4,5,7] is to calculate the availability according to the relationships (15): 
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KAC EM (15) 


i=l i=l 
where: 


i=l 


> T, —total sum of medium times of staying a vehicle in availability; 


> T, — total sum of medium times of staying a vehicle in unavailability . 
i=l 
The medium time of staying a vehicle in a considered state collection could be estimated by 


taking into consideration historical events[6]. Empirical results for the vehicles used in logistic air 
base system have been presented in tables 1-3. 


Tab. 1.Collection of data about medium number of entrance, medium transition probabilities, medium time of staying 
(hrs) and intensities of (l/year) S,-S,, states for a vehicle — C model 


State Sı S2 S3 S4 S5 So S7 Ss So Sio Su 
n. 214 161 330 8 4 2 2 358 | 45,656 2 416 
p;(n) 0,138 | 0,104 | 0,213 | 0,005 | 0,002 | 0,001 | 0,001 | 0,232 | 0,029 | 0,001 | 0,269 
L. 869,3 | 5656,3 | 13,75 252 22 16 240 | 0,2265 11 240 | 1439,43 
A. 10,077 | 1,5487 | 637,09 | 34,762 | 398,18 | 547,50 | 36,50 |38675,5 | 796,4 | 36,5 6,086 
Tab. 2.Collection of data about medium number of entrance, medium transition probabilities, medium time of staying 


(hrs) and intensities of (1/year) S,-S10 states for a vehicle — B model 


State Si S2 S3 S4 S5 S6 S7 Ss So Sio 
i, 173 161 165 4 2 1 l 179 23,65 1 
p,(n) 0,243 0,226 0,232 0,005 0,002 0,001 0,001 0,251 0,033 | 0,001 
i, 1159,3 | 5656,3 22 504 144 288 240 4,376 22 720 
A 7,556 1,548 398,2 17,38 60,83 30,42 36,50 2001,8 398,2 12,2 


Tab. 3.Collection of data about medium number of entrance, medium transition probabilities, medium time of staying 
(hrs) and intensities of (I/year) S,-Sy states for a vehicle — A model 


State Sı S2 S3 S4 Ss S6 S7 Ss So 
n, 173 161 165 4 2 1 1 179 22 
P; (n) 0,244 0,227 0,233 0,005 0,002 0,001 0,001 0,252 0,031 
i, 1159,33 6376,3 22 504 144 288 240 4,376 22 
A, 7,5561 1,3738 398,18 17,38 60,83 30,42 36,50 2001,83 398,18 
Where: 


n,— the average transition number from S; in 7, period; 
@,— frequency S; state in a considered collection state 7 — s (s = 9; 10 or 11); 


t,— the average duration time for S; state before entrance to all S; states in Ty period; 
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Â, - intensity of S; state; 


Mathematical dependencies for parameters mentioned above have been presented in [5] 
thesis. 


S1 S2 S3 S4 S5 S6 S7 S8 S9 S1 S 
Si(t) 


A model 
B model 
C model 


Rys. 5. Probabilities p;(t) of staying in availability for vehicles 


Fig. 5 describes probabilities staying in availability p;(ż) for vehicles in a considered S; states, 
which are calculated as: 


A 
p, =—, while availability has been calculated according to (15) relationship; 


To 
where: 
T; =t, +1, 


i=3 
4. Conclusions 


Estimated annual average sum of tenses £, +f, for vehicles population in availability is from 
6526 hrs (realistic C model) to 7535,63 hrs (optimistic A model — table 3), it means that the 
availability is from 0,745 to 0,86 in Tọ period. 

The availability of vehicle population is high, but there is also time margin in a year ty = 
1439,43 hrs (model C), resulting mainly from time losses during replacement S79 and renovating S4 
of the vehicles (tab. 1-3). It enables to improve the availability up to 16%. Therefore model C with 
Sı; state is justify and realistic. Withdrawal from usage decrease the number of vehicle population 
on average 7% in a year, but the replacement time is usually longer than one month. Szo state 
cannot be omitted because it significantly changes the duration time balance for all states of the 
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vehicles. Models B i C with S;, state reflect the simulation of renovating and replacement of the 
population in long periods of time. 

All presented models are stable and their non-ergodicity in real time is only according the 
simulation 0,2-1% for availability state after 3 years. The reason for their non-ergodicity can be 
vast (broad) spectrum of duration times for considered states collection. 

However Markov chains for A, B and C models are ergodic, but boundary probabilities p; 
concern collection states of the process and they cannot be interpreted in real (physical) time. 

It has been illustrated by the proportion of the estimated real time of refuelling max. ps(t)= 
4,376 hr (tables 2 and 3) to misinterpreted ergodic time ps * I year ~ 2200 hours for discrete time 
estimated from ps(n). 
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Abstract 


This work presents the main descriptors of the diagnostic system of machine exploitation. This problem contains: the 
measurements of technical state symptoms, the determination of their boundary values, and the frequency of diagnostics. 
This problem and the tasks in the system of machine exploitation supported with computer techniques constitutes about the 
rank of the discussed problem. 


Keywords: machine diagnostics, symptoms of the state, boundary value, frequency of diagnostics 


1. Introduction 

The knowledge of the machine’s technical state results from the need of making rational decisions 
about the “quality”, and further treatment of the machine. It can be a decision on further use, on 
undertaking appropriate preventive interventions, or on the introduction of construction, technology or 
exploitation changes. 

Discussed in this work are chosen problems of machine state diagnosis, emphasizing the problems 
of the new strategy of machine exploitation, including state symptoms, determination of the boundary 
value of measured symptoms, and determining deadlines of consecutive terms of diagnosis. 

The introduction of diagnostic systems makes it possible to improve the organization and 
management of machine usage in industrial institutions with the help of computer technology. 


2. Main problems in machine diagnostics 

The growing level of complexity of machines and criticality of their function within safety and 
economics, force constructors and users of these objects to know their current technical state, and to 
use them considering prognosis. This is possible, if at construction stage equipment and diagnostic 
procedures are integrated with the object. 


Generation of vibration signals in the description of machine condition changes. 

Evaluation of the machine dynamic state with the use of generated by them physical processes 
requires the association of functional parameters of the evaluated object with the set of measures and 
the opinions of output processes. 

While functioning of the machines, in the consequence of the existence of a number of external 
factors (environmental forcing, from different machines) and internal (aging, wastes, co-operation of 
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elements), there are disorders of balance state in the machine, which spread in a springy medium — the 
material of which the machine is built. 

The disorder has dynamic character and sustain the conditions of balance between the state of 
among the condition inertia, elasticity, the suppression and extortion. 

The disorder spreads from the source in the form of waves in the way dependent from physical 
properties, and the configuration borders, dimensions and shape of the machine. This in result causes 
wave energy scattering, their bending, reflection and mutual superposition. The existence of sources 
and spreading of disorders causes the occurrence of vibrations of machine elements and the 
surrounding environment. These processes are the basis of building a model of signal generation 
determining the manner of constructing, functioning and state changes of the object. 

The sequence of assumptions leading to the signal generation model can be presented with the use 
of a cybernetic model, as in the Fig. 1. 


Of forcing Characteristic Dynamic  y(t.8.r) 
internal signal setup Agnostic 

. pi(t.8,r) hit, 6,r) signal 

T- periodic 

transformation x (t.0) 
Of forcing kinematics 
outside 
m bguick interactions, — interactions 


( opetational > 


Fig.1 Model of diagnostic signal generation of a machine 


The presented way of interpreting the signal y(© ,r) is, in a general case of machines of periodic 
working, true, but not always as simple as in the Fig. 2. 


drive 


T - periodic 
transformation 
kinematics 


pi(t,0,r) y=Zh*gisdri 


Fig.2. Transformation of the characteristic signal g; into the output signal y 
as the model of signal generation in machines [2] 


The example of such formulation of the problem is the main transmission of the bridge of a 
vehicle whose generation model is presented in the Fig. 3. 

The received output signal in any place of the transmission casing is a weighed sum of answers to 
all elementary events Un(t, „r), occurring always in the same sequence in separate dynamic partial 
systems of impulse function of transition hn(t, © ,r). These effects, after going through certain 
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dynamic systems, sum up and undergo an additional deformation on the transmission corpse, whilst 
the change of place of receiving the “r” signal is also connected with the transmittance change. 

The notion n(t, ©) describes a random effect occurring because of the existence of dynamic micro- 
phenomena, such as friction, roughness, etc. The output signal of any receiving point can be 
approximately expressed with the equation: 


yx(8,r) = a (k) hj (t,0,r) * [uj (t,0,r) + n (t,6,r)] (1) 


where: the impulse transition function h(*) corrupts also the properties of the corpse, a(k) gives 
different weighs of summing connected with the receiving point ”7“*. 
Main problems of machine diagnostics include: 
- acquiring and processing diagnostic information; 
- building of models and diagnostic reports; 


trunk of thel 
gear 


Sequence of 
Orginal events 


possible places 
partial 


dynamic setup | 


of the reception 
of the signal 


| ot | 


attacking roller | X1 (t,0) 


- bearings 


attacking roller- 
wheel plate 


"| entrance of tht tool 2 
ino Tesses econstruction of 


the differencing 
mech. - bearings 
hs (1,0) 


axle shaft - Xn (t.0) | | 4 
differencing meclt 
Un (1,0) in (t, 


Fig.3. Model of the diagnostic signal generation of the toothed transmission 


u: (1,0) 


- the diagnostic inference and boundary value; 
- the classification of machine state; 
- the expected time of next diagnosis; 
- presenting the decision information. 

Measuring system for the aims of the modern diagnostics of machines consists of two basic 
elements: 
- the equipment in which distinguished are the modules: conditioning subsystem and processing of 
signals, the subsystem of processing signals of the phase gauge, the subsystem of the industrial 
computer, the subsystem of power supply, 
- the software, consisting of the modules: operating system (VxWorks), the software of the modules of 
signal processing and analysis, software assuring communication between the layers of the system, 
software for storage and archivization of measuring data, software managing the work of the system 
(configuration of the system, testing the system, the initialization of measuring sessions). 
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The introduced structure of the measuring system uses the newest solutions, both in hardware and 
software. The applied solutions easy extension of the system and including it for any diagnostic 
system. 

The problems troubling the practice of applying methods of diagnosis (Fig. 4): 

1. The time of forming the diagnostic symptom. 


„Sg (standards, DTR, assigning) 


Ir, 
generating o predicting of the symptom 


ti t t thst 0 


Fig.4. System diagnostic actions 


2. The change the boundary symptom - preventive system. 

3. Complex state evaluation: the measurements of symptoms, reference to the boundary value, 
prognosis of the state, estimation of next diagnosing, the genesis of the cause of the measured 
symptom changes. 


3. Diagnostic system of machine exploitation 

Contemporary machines determined by such features as: functionality, reliability, availability, 
safety, mobility and flexibility of operation. Formation and maintenance of these features is possible 
with methods of technical diagnostics which enables: 

- diagnostic construction and production of new machines and maintenance of machines in the 
condition of functional ability. 

The use of the machines is characterized by: 

- sets of randomly changing times of the correct work; 

- randomly changing moments of the beginning and changing lengths of times of task duration; 
- intensive work of people and machines in randomly changing period of exploitation time; 

- the influence of randomly changing conditions of exploitation; 

- various kinds of tasks executed in short periods of time. 

The needs and conditions of market economy justify the necessity to introduce modern authorized 
strategy of the machine production and exploitation. In the proposition of this strategy, we do not lose 
the so-far accomplishments of the newest strategy of exploitation according to the state, but it is 
creatively modernized. The proposed strategy of exploitation - ASEM — indicates the creator and 
responsible for the product by name. The manufacturer interested in the quality and later the sale is 
responsible for the product from the very idea, through construction, production and exploitation, until 
utilization after the liquidation of the object. 

The same manufacturer constructs and produces their products based on the newest achievements 
of technical thought, they protect their product with their own service during exploitation, and also 
they provide objects with diagnostic (preferably automatic) means. 
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The effectiveness of solutions in applying this strategy requires the improvement of: machine 
diagnostic models, the methods of diagnosing and prognosis of the machine state, economic, exact 
and reliable diagnostic equipment, the principles of the diagnostic vulnerability formation, the 
algorithms of conducting the maintenance of the machines in the ability state, methods of evaluating 
the efficiency of diagnosis and machine exploitation system. The specified questions embrace the 
whole problem and unambiguously establish the directions of the development of machine technical 
diagnostics. 

The users of machines are interested particularly in their task ability, for determining which the 
following are necessary: 

- determining the symptoms of ability state; 

- determining the boundary values of ability state symptoms, 
- determining the class of object ability, 

- determining diagnosis periodicity. 

The distinguished diagnostic tasks will be selectively discussed below, whilst their detailed 
description can be found in the author’s works [6,7]. 


4. Boundary values of state symptoms 

Task ability in the symptom notion is unambiguously described by the boundary value of the 
measured state symptoms. 

Exceeding the boundary value means that the machine enters the state of accelerated wear 
characterized by high probability of a sudden breakdown. 

Realized most often in the industrial practice passive and passive-active diagnostic experiments 
give state symptoms which are compared during concluding with boundary values available in many 
national, foreign, trade norms, or with data from own experience. When, however, there are no such 
norms for the examined machine, helpful can be the statistic description of a random exploitation 
process with the use of decay thickness or the frequency of the observed symptom occurrence. 

Estimating the symptom boundary value for a safe shutdown of a machine before failure can be 
realized with the use of statistical methods. 

The formula for determining Sgr minimizing the probability of failure, with the set, acceptable 
probability of redundant repair A can be written down in the form [1]: 


“| S 
P. {| — aS =4 
g (2) 
J X, 
where: Pg - the ability probability. 


According to Birger [1]: A= k (1-Pg), where: k - the coefficient of the store (k = 1-3 for usual 
damages, k = 3-10 for dangerous damages), Pg - availability the machines determined from the 
dependence: Pg=Nz / Nz + Nn, where: Nz - the number of fit machines, Nn - number of unfit 
machines. 

The row of simple transformations leads in the effect to the dependence: 


= F; 
A> =sto, T (3) 


The received estimation of symptom boundary value based on the mean value, dispersion and 
repair politics, gives good basis for simple determination of boundary values of examined state 
measures in the industrial practice. 
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5. Diagnosis periodicity 

The growth of intensity of the occurrence of damages along with wearing away the exploitation 
potential of the machine forces the need to optimize diagnosis periodicity. From the course damage 
intensity of the machine, it is concluded that that in the period of intensity growth it is necessary to 
increase the frequency of diagnosing. 

This helps to reduce expenditures on the exploitation of the machine (the decrease: wear 
intensiveness, waste of the fuel, spare parts, material use), and at the same time the costs and time 
consumption, and the machine turn-off time grow. 

The optimization of diagnostic periodicity comes down to two basic questions: how often to 
perform diagnosis?, in what range to perform the next diagnosis?. 

There are several possibilities to determine diagnosing periodicity (the method of symptom 
boundary values, the methods of the smallest sum of exploitation costs), whilst for their realization 
indispensable are numerous statistical data, often troublesome (in the sense of amount and reliability) 
to obtain. 

In this work, the question of diagnostic periodicity has been considered symptom notion, using a 
known symptom boundary value. Performing n-measurements of a chosen in a separate procedure 
signal measure (symptom), and determining on their basis a boundary value according to the 
dependence (2), there is a need to determine the date of the next diagnosing ta. The essence of the 
method presented in the works [6,7] show that the date of next diagnosing is estimated from the 
dependence (Fig. 5): 

1-P)(S.-S 
i= en) O, (4) 


m 


Fig.5. Periodicity of diagnosing in the symptom notion 


6. Management of exploitation system 

All of economic organizations have a certain system of management which fulfills its 
requirements within the realization of an accepted strategy. This is particularly essential for those 
which have a considerable influence on the course of production process (logistics, exploitation, tools 
and equipment), or control its fixed assets of considerable, from the point of view of the company, 
values (movement maintenance, repairs, check-ups) [4]. 

Subsystem functions: 

- it conducts the classification and record of all fixed assets, 
- it proposes basic technical-economic indexes, 
- it supervises the exploitation of fixed assets, 
- it analyses data from monitoring and makes decisions, 
- it infers the liquidation of fixed assets, 
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- it plans, supervises and realizes all kinds of examination, maintenance and repairs, 
- it establishes basic norms, records and accounts for performed works, 
- it plans the supply of spare parts and necessary materials for repairs, 
- it offers and motivates leasing, offers and motivates outsourcing, 
- it organizes storage of spare parts, their distribution and accounting for, 
- it plans investment tasks, organizes and realizes the purchase of machines and equipment, 
- it organizes the receipt of fixed assets, 
- it prepares the technologies of repairs. 
Analyzing the range of functions attributed to the system for realization, it is possible to determine 
what groups of data ought to enter it, as well as what data it generates. 
The model of exploitation system management was built on the basis analysis of two basic 
criteria, i.e. the flow and type of data, and functions realized by individual modules. The structure of 
exploitation system management, together with the flow of data, is shown in the fig. 6. 


MODULE FOR 
THE FOREGIN 
REALIZATION 5 
MODULE FOR THE 
REALIZATION OF TASKS 
PE4 


Fig.6. Model of exploitation system management 


Separate modules creating the structure of the system realize the following functions: 
1. The module of data processing is responsible for processing data sent to the system. Carriers and 
media of the transmission can be considerably diverse. 
2. The managerial module, to which data of various degrees of aggregation enter from the module of 
data processing. It can be stated that about 80% of the data is processed according to the assigned 
algorithms, creating a basic set for the need of SE management. 
3. The logistic module, which delivers indispensable materials, equipment, components and 
standardized machine elements for the needs of realized repairs; manages stock and analyses the level 
of stores, runs the record of distributed supplies, analyses their waste for individual orders, organizes 
and supervises transport of purchased technological equipment, co-operates with company’s logistics 
in the range of the economy of scrap-iron. 
4. The task realization module which realizes or supervises the realization of examinations and a 
bigger part of repairs. 
5. The module of the strange realization, registering the range of repair works ordered to a third party. 
6. The control module, checking the quality and range of realized works, outside and own. 
7. The module of technical base renovation, purchase of machines, renovation and repairs. 
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8. The personnel training module, motivation and training of workers. 
9. The accounting module, creates abstract summaries, as well as controls the needs for the realization 
of certain operations. 
10. The technical module realizing the functions: planning, constructional, technological, technical 
state evaluations of the possessed equipment, record and updating, the emission of records. 
PU - the exploitation subsystem which exploits machines and equipment. 
PP - Remaining subsystems. Relations of these subsystems with the exploitation management system 
are defined to a smaller degree. 

The construction of the model of the exploitation management system allows to identify the basic 
elements of its surroundings, as well as of the modules creating the subsystem itself. 


7. Summary 

The accomplishment of diagnostics in recent years using the achievements of many fields of 
science, allow treat it as a tool of formation and evaluation of machines, at all stages of their 
existence. 

Looking at present trends of machine development it should be recognized that currently the 
growth of their quality is contained mainly in the sphere of automation. Automatic acquisition of 
measurable features is becoming the only objective way of evaluating and forming the quality of 
machines. 

The range of investigations in the field of methodology of diagnostics includes such questions as: 
the source of diagnostic information, signals and diagnostic symptoms, the principles of detailed 
methods of diagnostics, modeling in diagnostics, diagnostic experiments, supporting diagnostics with 
modern computer technologies, diagnosing in systems of human engineering and social engineering, 
and the organizational and economic aspects of applying diagnostics. These question respectively 
apply to: source of information from the physical side and from the informative side, further bases of 
methods and investigative techniques, simulation and experimenting in diagnostics, and modern 
inference and visualization of worked out diagnostic-exploitation decisions. 
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Abstract 


In this work presented are chosen problems of machines’ technical state diagnosis with the use of identification and 
technical diagnostics methods. Relations between methods of dynamic state evaluation and methods of technical state 
evaluation were indicated. Example modal analysis results illustrate the complexity of projecting dynamic state 
researches into diagnostic researches of machine state evaluation. 
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1. INTRODUCTION 


Destruction processes of technical systems force the need to supervise changes of their 
technical state. It is possible with the use of technical diagnostics methods. 

Methods and means of modern technical diagnostics are a tool of machine state diagnosis, 
which is the basis of decisions made at each stage of their existence. 

Many previous works of the author [1,2,3,4,5] clearly indicate connections between machine 
dynamics and technical diagnostics, especially vibration diagnostics. The bases of identification, 
modeling and concluding fully convince towards the dominating role of vibrations in machine 
state identification [6,7,8,9]. 

Properly planned and realized experiment is the base to obtain diagnostically sensitive signals 
which processed will determine state diagnosis procedures. The processing includes: creation of 
numerous signal measures in time domain, frequencies and amplitudes, selection and reduction of 
the number of signal measures, creation and analysis of effectiveness of cause-and-effect models, 
as well as evaluation of the righteousness of made diagnostic decisions. 

The realization of these tasks is possible only through broad support from information 
technology, which in this work is presented in the form of SIBI programs. 

Practical applications of the presented ideas has been verified in researches on complex objects 
exploited in difficult climate conditions. 


2. DYNAMICS AND DIAGNOSTICS 
Into quality measures of machine’s technical state, i.e. its dynamics, included is the level of 


vibration amplitudes, as well of the machines as the lot, and also of relative vibrations of separate 
elements and parts. Overall vibrations of the machine can be perceived as an external symptom 
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while they are responsible for the level of interferences emitted into the environment. Relative 
vibrations of separate elements, however, influence mainly the state of internal forces in the 
machine, i.e. at its level of dynamic stress amplitudes. 

Identification can concern both the construction of models and the reconstruction of the 
examined model state, which leads straight to the problem of technical diagnostics. 

The process of diagnostic identification includes: 

* modeling (symptom or structural), 

* identification experiment (simulation and/or real), 

* estimation of diagnostic parameters (state features or symptoms), 
* diagnostic concluding. 

The specificity of diagnostic identification tasks is different from general identification in the 
way that it includes a number of additional elements enhancing this process. These are: 
- constructing models of diagnostic signals generation, 

- choosing features of object structure state and diagnostic symptoms, 
- modeling cause-and-effect relations, 

- evaluating the accuracy of choosing variables for the model, 

- determining boundary values of the symptoms, 

- classifying the states and determining diagnosis periodicity. 

Methods of identification can be divided concerning: the kind of identified model, the kind of 
experiment, identification criterion applied, as well as estimation procedure applied. In general 
these are: methods of analysis, time, frequency, correlation, regression, factor analysis, as well as 
iteration methods described in works of many authors [2,4,6,5,8]. 

For simple objects, a good tool to evaluate their changeable dynamic state are methods of 
simple identification which use amplitude-frequency spectrum. Searching rezonans frequency and 
amplitude value in this frequency with the use of tests (impulse, harmonic and random) are 
relatively well masterem in research techniques of our enterprises [2,5]. 

Another way of describing and analyzing the dynamic state of machines is a modal analysis 
used as a theoretical, experimental and exploitation method. It uses frequencies of own vibrations, 
values of suppression and forms of vibrations to describe the changing machine state, and it is 
used to improve the finished elements method. The presented procedures are based on the 
knowledge of the system model, and the conclusions drawn from the actions on the models depend 
on their quality. Depending on the aim of the performed dynamic analysis of the object, different 
requirements are set for the constructed models, and their evaluation is conducted with different 
experimental methods. 


3. DESCRIPTION OF OBJECT STATE CHANGES 


The dynamic state of the object can be, in the easiest case, described with a model of | degree 
of freedom — Fig.1. A conventional description of this model are known relations (1-4) indicating 
that vibrations well reflect the state of the machine. A description of this model can be achieved 
within m, k, c categories, or through a, v, x researches. 


Fig.l. Model of a system of 1 degree of freedom. m, k, c = vibration process (a, v, x) 
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Identification of his model (1) from the experimental side is the a, v, x measurements for 
different time moments, which reflects the changes of the object state and is widely applied in 
vibration diagnostics. The solution of the task in the m, k, c, categories, however, requires a 
number of solution conversion of the equation (1) for determining: 

Cy = 2M0 c, =4llmf 
k=m-o° k = 411’? mf? (5) 
Determining the value (5) requires realizing identification experiment from which the frequency f 
or frequency w can be determined. Here is useful the simple identification or modal analysis 
directly giving the values of own frequencies w from the stabilization diagram — Fig.2. 
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Fig.2. Stabilization diagram for œ determination 


The problem becomes more complicated for models of many degrees of freedom (more than 
3). Here also the problem of object state identification can be solved from the measurement side 
(a, v, x), Chile from the side of determining m, k, c own problem needs to be solved. 


2 
(K-o M)-q, =0 (6) 
Equation (6) presents a linear system of homogeneous algebraic equations: 
2 2 
(e, -0 M, i + lk, -0°m, y, pot (k, 70 M, y, =0 
2 2 2 
len -0 myk, tl, -0 mk, tot (by, 70 m, |, =0 
ae de dai teense C (7) 
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A solution for qź0exists when the main matrix determinant (K -0M)=0, i.e. 


det(K —0M )= 0. Solving the system of equations (7) own values can be determined, and from 


them the frequencies of own vibrations, indispensable for the object identification (A = @* = S; 
m 


4. IDENTIFICATION RESEARCHES SOFTWARE (SIBI) 


More and more frequently conducted identification researches of machine dynamic state, used 
for the evaluation of the state changes, fault development and location of the occurred state causes, 
were the basis for creating a specialized software system. It allows acquiring and processing 
measurement data, creating many measures of diagnostic signals, examining their diagnostic 
sensitivity, statistic processing and diagnostic concluding. The program was named Information 
System of Identification Researches (System Informatyczny Badań Identyfikacyjnych — SIBI). 

The structure of the program is a module construction which includes the following modules 


(Fig.3.): 
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Rys.4.1. Główne okno dialogowe programu SIBI 
Fig.3. Main dialog box of SIBI program 


A. READ UNV module which allows processing data from UNV format into XLS format. 

B. SYMPTOMS module which allows defining, determining and creating matrix of many 
measures of vibration processes. 

Modules A and B are part of software responsible for acquiring and processing vibration processes 

in order to acquire observation matrix of vibration estimators. 

C. OPTIMUM module uses the method of ideal point for individual evaluation of the sensitivity 
of measured symptoms of vibration processes. 

D. SVD module used for determining generalized damage measures, and for the evaluation of 
damage development. Using the SVD method allows a multidimensional description of the 
state of the examined object. 

E. INPUT/OUTPUT module used for the analysis of similarities between vibration processes, 
and for determining different exploitation measures of the examined object. 

F. OPTIMIZATION module used for creating models and data in genesis (with approximation 
and interpolation methods), diagnosis and prognosis of object states. 

G. NETWORK module using neuron nets for state classification on the basis of obtained results 
in the form of time rows. 

C, D, E, F, G modules are elements of 2 parts of the software allowing the performance of statistic 

concluding and cause-and-effect relations, as well as visualization of the obtained results. 
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5. CAUSE-AND-EFFECT MODELING 


Many state measures acquired in experiments requires the reduction of over measurement, 
which is possible with the use of OPTIMUM procedure (statistic evaluation of separate measures) 
or SVD (for the multidimensional approach). Optimized set of symptoms is the basis of 
constructing cause-and-effect, most often regressive, multidimensional models (Fig.4). 


y = -2,68061w1 - 0,54083 row] - 0,49318 x1 + 2,01273 w2 + 0,35480 row2 + 2,26940 x2 - 
0,02717 H(f) + 0,06833 H(f)L + 0,01696 g2xy - 92,00391 ARMS(t) + 12,99146 bkurt + 
239,69713 Cs - 200,58670 I - 44,37385 


Dependent variable ~~ Real data 

Model 

R2=0.83343 

F=9.173 _ Fstat=2.2258 
Pvalue=5.237e-006  a=0.05 
JBtest: OK 
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Fig.4. Regressive model determination 


The wellness of the model is evaluated with the help of the determination coefficient R”, and 
the number of component symptoms determines its accuracy — Fig.5. 


5 etymatorów 
13 etymatorów 


R° = 0,0857 


Wartości 
Wartości 


Estymatory Estymatory 


Fig.5. Number of measures v. accuracy of the model 
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6. CONCLUSIONS 


Presented in this work considerations concern the modeling of object dynamic state with the 
use use of description and researches within the range of identification, distinguishing modal 
analysis and ideas directly supporting different methods of forming machine dynamics. 

The knowledge of the dynamic state and structure of the system allows to describe its 
behavior, and allows creating prognosis models of the system behavior in the function of dynamic 
evolution time, based on the model of the technical state symptoms growth. Most often, however, 
there are no known equations describing behaviors of the system in the function of dynamic 
evolution time, which accounts for the need to apply new tools to examine the dynamic state. 
There is, therefore, the requirement to experimentally verify analytical technical models as the 
proper one is a model which is verified in practice. An experiment is, therefore, often only an 
inspiration for further researches leading to the optimization of the construction. 
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